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20.  The  scattering  of  acoustical:  waves  from  fluid-load  elastic  spheres  with 

10  <  ka  <  100  (where  k  is  the  wavenumber  of  the  acoustic  wave  in  the  liquid 
and  a  the  radius  of  the  sphere)  is  studied.  The  main  emphasis  is  on 
understanding  the  scattered  pressure  near  backscattering.  By  carrying  out  a 
Sommerf eld-Wataon  Transformation  (SWT)  it  is  shown  that  this  pressure  Includes 
contributions  from  specular  reflection,  transmitted  bulk  waves,  and  surface 
waves.  It  is  shown  experimentally  that  surface  wave  contributions  to  near 

i 

backscattering  display  a  Bessel  function  angular  dependence  characteristic 

of  weak  axial  focusing  along  the  backscattered  direction;  l.e.,  glory 
scattering.  This  focusing  is  modeled  by  alternate  methods  which  agree  and 
whieh> complement  each  other.  One  method  allows  a  simple  picture  of  the  glory 
phenomena  while  the  others  (which  includes  the  SWT)  are  more  mathematically 
rigorous.  The  absolute  value  of  individual  Rayleigh  and  whispering  gallery 
surface  wave  contributions  is  found  by  this  modeling„>to  be  equal  to  A  J^(kbf)  c 
-  Tn-thlaexpreSalon  \  ¥*  is  the  angle  relative  to^badescattering  and  b  and  A 
are  parameters  yhich  depend  on  the.  surface  wave  contribution  being  examined  and 
are  determined  by  the  models. The  theoretical  results  £ec — A — and — dr^are 
tested  experimentally  using  the  Rayleigh  surface  wave  contribution  to 
scattering  from  a  tungsten  carbide  sphere.  The  SWT  results  are  also  used  to 
find  the  far  field  backscattering  pressure  as  a  function  of  ka.  These  results 


in  turn"  are  used  to  better  understand  and  interpret  the  recently  developed 

.  0-""’ 

Resonance  Scattering  Theory  (RST)  as  applied  to  acoustical  scattering.  It  is 
shown  that  surface  wave  contributions  to  backscattering  may  be  summed  in  a  way 
analogous  to  the  frequency  response  of  a  Fabry-Perot  resonator.  /  ■  "  '  ) 
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ACOUSTICAL  SCATTERING  FROM  AN  ELASTIC  SPHERE  IH  HATER:  SURFACE 
HAVE  GLORT,  RESONANCES ,  AND  THE  SOtMERFELO-HATSOH 
TRANSFORMATION  FOR  AMPLITUDES 
Abstract 

by  Katin  Lm  Villi saa,  Ph.D. 

Wan h lug ton  Scats  UnlTarsity,  August  1985 

Chair:  Philip  L.  Mars ton 

Th«  scattering  of  acoustical  waves  from  fluid- load  elastic  spheres 
with  10  <  ka  <  100  (where  k  is  the  wavenumber  of  the  acoustic  wave  in  the 
liquid  and  a  the  radius  of  the  sphere)  is  studied.  The  main  emphasis  is  on 
understanding  the  scattered  pressure  near  backscattering.  By  carrying  out  a 
Sommerfeld-Watson  Transformation  (SWT)  it  is  shown  that  this  pressure  includes 
contributions  from  specular  reflection,  transmitted  bulk  waves,  and  surface 
waves.  It  is  shown  experimentally  that  surface  wave  contributions  to  near 
backscattering  display  a  JQ  Bessel  function  angular  dependence 
characteristic  of  weak  axial  focusing  along  the  backscattered  direction;  l.e., 
glory  scattering.  This  focusing  is  modeled  by  alternate  methods  which  agree 
and  which  complement  each  other.  One  method  allows  a  simple  picture  of  the 
glory  phenomena  while  the  others  (which  includes  the  SWT)  are  more 
mathematically  rigorous.  The  absolute  value  of  individual  Rayleigh  and 
whispering  gallery  surface  wave  contributions  is  found  by  this  modeling  to  be 
equal  to  A  J^(kby).  In  this  expression  y  is  the  angle  relative  to 
backscattering  and  b  and  A  are  parameters  which  depend  on  the  surface  wave 
contribution  being  examined  and  are  determined  by  the  models.  The  theoretical 
results  for  A  and  b  are  tested  experimentally  using  the  Rayleigh  surface 
wave  contribution  to  scattering  from  a  tungsten  carbide  sphere.  The  SWT 
results  are  also  used  to  find  the  far  field  backscattering  pressure  as  a 


function  of  ka.  These  results  in  turn  are  used  to  better  understand  and 
interpret  the  recently  developed  Resonance  Scattering  Theory  (RST)  as  applied 
to  acoustical  scattering.  It  is  shown  that  surface  wave  contributions  to 
backscattering  nay  be  summed  in  a  way  analogous  to  the  frequency  response  of  a 
Fabry-Perot  resonator. 
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1 .1  Overview 

Many  class  la  physics,  results  derived  from  investigating  "canonical 
problems"  are  important  in  understanding  more  difficult  situations.  A 
particular  example  of  this  is  the  scattering  of  waves  incident  on  a  sphere. 
This  problem  is  of  interest  in  such  areas  as  electromagnetics,  acoustics, 
geophysics,  and  particle  physics.  Because  of  this,  it  has  been  Che  subject  of 
much  analysis.^  As  early  as  1863  Clebsh  addressed  the  problem  of  acoustical 
wave  scattering  from  a  rigid  sphere.  In  the  ensuing  years  from  1863  to  World 
War  II  many  people  examined  various  versions  of  the  problem  in  many  different 
fields.  Reference  1  1s  a  survey  of  many  of  these  analyses  up  to  1941. 

Formal  solutions  to  the  various  problems,  in  terms  of  an  infinite  Partial  Wave- 
Series  (PWS)  of  spherical  functions,  are  obtainable  by  using  the  separation  of 
variables  technique  in  solving  the  differential  equations.  In  many  cases  the 
PWS  results  give  little  physical  insight  into  Che  problem  at  hand.  This  is 
particularly  true  if  the  sphere  is  much  larger  than  the  wavelength  of  the 
plane  wave.  This  is  equivalent  to  saying  the  dlmenslonaless  parameter  lea  of 
the  sphere  is  much  larger  chan  one  where  k  is  the  wavenumber  of  the  plane 
wave  and  a  is  Che  sphere's  radius.  Much  of  Che  work  even  to  this  day 
revolves  around  obtaining  a  better  understanding  of  observable  physical 
phenomena  by  applying  various  mathematical  techniques  to  Che  PWS  results  when 
ka  »  1.  From  this  type  of  analysis  one  can  obtain  results  equivalent  to 
Chose  of  less  rigorous  methods  such  as  geometrical  optics .  For  the  problem  we 
shall  study  (see  below)  we  will  first  use  the  less  rigorous  methods  directly 
since  they  have  less  of  a  tendency  to  mask  the  physics  of  the  situation.  We 
will  then  attack  the  problem  by  beginning  with  the  PWS  results  which  are  also 
given  below.  Finally  we  will  examine  the  interpretation  of  a  recently 
developed  theory  which  uses  the  PWS  directly.  In  the  paragraphs  below  we 


shall  specify  the  problem  we  will  examine,  give  a  short  summary  of  its  PWS 
solution  and  discuss  work  which  had  been  carried  out  up  till  1983  which  was 
important  to  analysis  in  this  dissertation. 


The  problem  we  wish  to  study  is  the  scattering  of  an  acoustical  plane 
wave  incident  from  an  nonvlscous  fluid  onto  an  Isotropic  elastic  sphere.  We 
are  most  interested  in  the  total  pressure  external  to  Che  sphere  in  the  range 
10  <  ka  <  100.  The  PWS  solution  for  this  pressure  was  found  first  by  Faran  in 
1951.“'  Related  work  up  to  chat  point  is  summarized  in  Reference  2. 

To  solve  che  problem,  Faran  wrote  the  vector  displacement  u  in  Che 
elastic  solid  in  terms  of  a  scalar  potential  p  and  a  vector  potential  A. 

In  terms  of  these  potentials  Che  vector  displacement  u  is  given  by 

u  ■  -Vp  +  7  x  A  (1) 


where  7  is  the  dal  operator.  The  potentials  in  turn  must  satisfy  the 
equations 


V2ip 


_1_  ii 

2  3.2 
cT  9t 


7  x  7  x  A 


(3) 


where  and  c,,  are  the  longitudinal  and  shear  elastic  wave  speed 

respectively.  Assuming  a  plane  wave  incident  as  shown  in  Fig.  1  symmetry 


considerations  allow  one  to  show  that  there  will  only  be  an  A^  component  to 


A  and  that  p  and  A^  will  have  no  $  dependence.  Therefore  solving  Eq. 


(2)  and  (3)  by  the  method  of  separation  of  variables  one  finds  that  'p  and 
A,  may  both  be  written  as  an  infinite  PWS  with 


Th«  p,  9,  «  coordiaaea  syseaa  show  abov*  wma  us*d  ia  vrlciag  cfa* 
parti*!  wav*  sarlaa  (?WS)  solution  co  plana  wav*  scattaring  from  an 

alaaclc  sphar*.  Tha  plan*  w*w*  ta  aaauaad  co  ba  eravaling  ia  cha  +z 
diraction. 


i|Kr,9,c) 


l  a  j  (k.a)  P  (cose)#"3"' 
n-0  a  n  T  Q 


<*: 


A^(r,9,c) 


l  h  j  <k  a)  ?hcosB)e~imt 
n-0  n  a  S  n 


(5] 


where  w  is  Che  angular  frequency  of  che  incident  plane  wave,  a  and  b 

n  a 

are  unknown  constants,  jQ  is  Che  spherical  Bessel  function,  k^  and  kg  are 
che  longitudinal  and  shear  wavenumbers,  a  is  the  radius  of  the  sphere,  and 
Pr  and  P  are  the  Legendre  and  associated  Legendre  functions  respectively. 
The  pressure  in  the  liquid  satisfies  the  equation 


V2p 


:2  3t2 


(6) 


where  c  is  che  acoustic  wave  speed  for  Che  liquid.  The  vector  displacement 
in  the  liquid  is  given  by 


^  -  -i_  7p  (7) 

puT 

where  P  is  the  density  of  che  liquid  and  uj  is  che  angular  frequency  of  the 
acoustic  wave  in  che  liquid.  The  total  pressure  p  is  given  by  the  incident 
pressure  p^  of  the  plane  wave  plus  the  scattered  pressure  pg;  i.e., 

P  ■  P.  +  P  •  Assuming  an  unit  amplitude  Incident  plane  wave  p  can  be 

IS  1 

written  as 

00 

p  (r,9,t)  -  l  (i)a(2n+l) j  (kr)P  (cos8)e-iu,C  (8) 

n-0  n  n 

whereas  the  scattered  pressure  is  a  combination  of  outgoing  spherical  waves 
and  is  given  by 

-  I  c  h(1)(kr)P  (cos9)e~iu,t 


P„(r,9,c) 


m 
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In  Eq.  (8)  and  (9)  k  is  cha  wavenumber  for  Che  pressure  waves  In  Che  liquid. 
In  Eq.  (9)  h^1^  is  cha  spherical  Hankel  of  che  first  kind  and  che  are 

unknown  constants. 

Now  to  find  che  a  ,  b  ,  and  c  we  must  apply  Che  appropriate 
na  ti 

boundary  conditions.  The  boundary  conditions  are;  che  displacement  uf  (the 
radial  component  of  u)  is  continuous  across  che  boundary,  che  pressure  in  che 
liquid  must  equal  che  normal  component  of  Che  stress  in  the  solid  sphere  at 
the  boundary,  and  the  tangential  (shear)  stress  components  for  Che  solid  muse 
vanish  at  che  boundary.  The  mathematical  form  of  these  boundary  conditions 
may  be  found  in  Ref.  2.  We  shall  not  carry  out  Faran' s  analysis  to  find  che 
an*  <5a  ’  an<*  cn*  w®  simply  quote  che  result  for  che  total  pressure  external 
co  che  sphere  in  a  form  appropriate  for  our  purposes.  The  total  pressure  is 


p(r,0) 


go 


l 

n-0 


itl(2n+l)[jQ(kr) 


B  (ka)  "  n  s 

+  r7TThu;(kr)]P  (ccs0) 
0  (ka)  n  n' 


(10) 


In  the  equation  above,  B  and  D  are  3x3  determinant  whose  elements  are 

Q  a 

given  in  Appendix  A  of  Chapter  3.  In  Eq.  (10)  the  time  dependence  e  i(iJC  is 
not  shown  explicitly. 

In  1962  Robert  Hickling,  using  Faran' s  results,  published  an  article 

in  which  he  carried  out  a  computer  analysis  of  che  far  field 

backscactered  pressure  from  elasdc  spheres  made  of  several  different 
4 

materials.  The  spheres  had  ka  in  che  range  0  £  ka  £  60.  The  case  of 
backscacterlng  (i.e.,  3  "it)  is  particularly  important  in  such  areas  as  sonar 
and  non-destructive  testing.  The  analysis  of  this  dissertation  emphasizes  che 
development  of  an  understanding  of  the  scattering  pressure  near  0  *  it. 

Hickling  obtained  the  steady  state  backscacterlng  pressure  as  a 
function  of  the  dimensionaless  parameter  ka  of  the  spheres.  A  diagram 


■  »  '  >  '  i''  ^ 
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equivalent  Co  hie  steady  state  results  for  Armico  Iron  is  given  in  Fig.  2a. 

Re  used  these  steady  state  results  along  with  Fourier  transform  techniques  to 

determine  che  backscattsring  echo  from  an  Armico  iron  sphere  which  was 

subjected  to  a  tone  burst.  Redrawn  in  Fig.  2b  is  Sickling's  result  for  che 

% 

case  of  a  5  cycle  pulse  incident  on  an  Armico  iron  sphere  with  ka  *  24.5.  The 
figure  reveals  a  decaying  pulse  train.  About  a  year  before  Hickling's 
analysis  Hampton  and  McKinney  had  published  an  experimental  paper  in  which, 
among  other  things,  they  had  subjected  an  aluminum  sphere  to  a  short  acouaclc 
pulse  of  the  form  used  by  Hickling  in  his  computer  analysis.  A  pulse  train 
type  echo  was  clearly  evident.^  Also  evident  after  reading  both  articles  is 
the  fact  that  much  of  the  phenomena  seen  were  incompletely  understood .  In 
fact  the  PWS  solution  is  not  a  particularly  well  equipped  vehicle  through 
which  to  understand  the  physical  origins  of  che  echoes.  We  turn  now  to  work 
that  allows  an  alternate  avenue  through  which  to  understand  the  backscattered 
pressure  from  che  elastic  sphere.  It  should  be  reaffirmed  chat  we  are 
examining  work  which  was  relevant  to  the  methods  used  in  the  current 
dissertation  and  make  no  claim  that  the  present  summary  is  exhaustive. 

In  the  middle  1950's  Joseph  Keller  and  associates  developed  a  new 
6  7 

theory  of  scattering  *'  which  they  called  the  Geometrical  Theory  of 
Diffraction  (GTD).  The  theory  is  based  on  the  postulate  that  fields  propagate 
along  rays  as  in  geometrical  acoustics,  but  introduces  diffracted  waves.  Like 
geometrical  acoustics  it  Is  applicable  to  scattering  from  a  sphere  with 
ka  »  1.  An  article  by  Levy  and  Keller  in  1959  is  of  particular  interest 
since  it  applied  GTD  to  acoustical  scattering  from  impenetrable  spheres  and 
cylinders.  For  the  problems  discussed  in  the  1959  paper  the  diffracted  rays 
propagated  in  part  along  the  surface  of  the  scattering  object  and  were 
therefore  called  surface  waves.  The  GTD  surface  wave  picture  is  developed  in 


Calculations  of  eha  backscattarad  pressure  from  an  Armico  Iron 
sphere:  a)  The  backseat taring  form  function  as  a  function  of  ka. 
This  figura  was  calculated  vising  cha  program  CHIVES,  given  in  eha 
appendix,  b)  Using  cha  rasulcs  vised  co  generate  figura  Fig.  2a  and 
Fourier  eransf oraatlon  techniques ,  Sickling  obtained  cha 
backseattaring  response  to  a  3  cycle  eone  burst  shorn  hare.  The 
horizontal  claa  scale  is  rsfarencsd  eo  eha  canter  of  a  pulse 
traveling  in  uatar  froa  (r-»,  9*O"*(r«0)-»(r«»,  9  ■  r) 
the  aondimanslonal  cima  unit  t  »  ct/a. 
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Levy's  paper  by  first  using  ideas  similar  to  those  of  geometrical  acoustics. 

A  topical  trajectory  for  a  diffracted  ray  in  going  from  a  source  point  Q  to 
a  receiver  P  is  shown  in  Fig.  3*  The  surface  ray  is  launched  at  the  point 
whose  location  can  be  determined  by  phase  matching  conditions  (cf.  Chapter 
2) .  The  surface  wave  reradiates  energy  into  the  fluid  continuously  as  it 
travels  around  the  sphere  (indicated  by  the  dotted  rays  in  Fig.  3).  At  the 
point  Q2  energy  is  reradiated  into  the  liquid  in  the  direction  of  the 
receiver  P.  The  surface  wave  picture  developed  by  use  of  geometrical 
arguments  was  recalculated  more  rigorously  later  in  Levy's  paper  by  using  the 
Sommerfeld-tfatson  transformation  (SWT). 

The  SWT  consists  of  rewriting  the  PWS  solution  of  the  problem  under 
study  in  terms  of  a  contour  integral.  The  contour  is  then  deformed  to 
surround  poles  in  the  complex  plane  which  eventually  lead  to  surface  wave 
contributions.  The  contour  can  be  also  be  separated  and  further  deformed  to 
pass  through  saddle  points  of  the  complex  integral.  A  saddle  point  analysis 
then  gives  contributions  to  scattering  which  are  familiar  from  geometrical 
acoustics;  for  example,  specular  reflection.  The  SWT  was  first  used  by  Watson 
and  then  further  developed  by  many  others.^ 

In  Levy' s  paper  the  SWT  analysis  occupies  only  a  few  pages  and 
contains  little  detail  since  it  was  not  Che  primary  goal  of  Che  paper  Co 
understand  all  of  the  many  effects  which  can  be  seen.  In  contrast  Co  chis 

g 

Nussenzveig,  in  a  1965  paper,  embarked  on  a  very  derailed  analysis  of  scalar 
plane  wave  scaccarlng  from  a  sphere  using  boundary  conditions  appropriate  for 
an  acoustically  soft  sphere.  He  used  a  modified  SWT  analysis  to  investigate 
Che  PWS  solution  for  a  wavefunction  which  could  be  interpreted  as  either  the 
velocity  potential  for  an  acoustically  soft  sphere  or  as  the  Schrodinger  wave 
function  for  a  hard-core  potential.  At  the  time  of  his  paper  the  SWT  was 


A  typical  trajectory  for  a  diffracted  ray  of  Che  GTS  in  going  from 
cha  sender  0  to  receiver  ?.  A  surface  wav®  is  launched  at  point 
.  As  indicaead  by  cha  dottad  Linaa  cha  surface  wave  continuously 
rsradiacaa  anargy  into  cha  fluid.  At  cha  point  0^  anargy  is 
radiated  in  cha  direction  of  ?.  As  drawn  hara  cha  ray  is  similar  co 
ones  used  later  in  this  dissertation.  In  Levy's  paper^  rays  ware 
assumed  to  be  launched  onto  cha  sphere  tangentially. 
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b«ing  used  la  particle  physics  In  connection  with  Reggi's  work  on  complex 
angular  momentum.® 

Later,  In  1969,  Nussenzveig  authored  two  papers  which  treated  a  sphere 
with  boundary  conditions  appropriate  for  Che  acoustical  problem  of  a  liquid 
sphere  embedded  In  another  liquid.  The  second  of  chese  papers  was  concerned 
partly  with  understanding  a  scattering  effect  known  as  glory.  The  term  glory 
refers  to  a  weak  axial  focusing  of  scattered  waves  along  either  the  forward  or 
backward  scattering  directions. More  will  be  said  about  glory  scattering 
lacer  in  this  Introduction. 

The  results  of  Nussenzveig1 s  papers  included  surface  wave 

contributions  in  agreement  with  the  ideas  of  the  GTD.  One  obvious  next  step 

in  studying  scattering  from  a  sphere  would  seem  to  be  a  generalization  of 

Nussenzvelg' s  papers  to  Che  case  studied  by  Faran;  i.e.,  an  Isotropic  elastic 

sphere  embedded  in  a  nonvlscous  liquid.  One  could  foresee  a  possible 

profitable  Interaction  between  experimental  results  such  as  Hampton's’’  and 

theoretical  results  of  tha  SWT.  This  dissertation.  In  fact,  represents  a 

partial  fulfillment  of  this  next  step. 

It  may  at  first  seem  curious  chat  the  generalization  and  interaction 

between  experiment  and  theory  described  above  did  not  happen  before  now. 

However,  part  of  tha  reason  for  the  delay  can  be  found  by  realizing  that 

there  is  another  similar  problem  which  is  in  some  ways  easier  to  study.  That 

is  che  problem  of  plane  wave  scattering  from  an  isotropic  elastic  cylinder  in 

a  liquid.  In  Che  same  time  frame  as  the  studies  discussed  above,  other 

12  13 

investigators  were  developing  equivalent  results  for  cylinders.  ’  Indeed, 

in  1968  (a  year  before  Nussenzveig' s  study  of  the  glory  in  Ref.  10)  Doolittle, 

Uberall,  and  Uginclus,  published  a  SWT  analysis  of  scattering  from  an  elastic 
13 

cylinder.  Some  of  the  work  (by  them  as  well  as  ocher  researchers)  leading 
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up  co  this  analysis  is  referenced  in  Che  incroduccion  of  Ref.  12.  Elastic 

cylinders  wich  ka  <  100  (where  a  is  now  che  radius  of  che  cylinder)  became 

Che  objecCs  of  sCudy  chrough  which  people  Increased  chair  underscanding  of  che 

surface  waves  of  che  SWT  (ac  lease  in  acouscics  research).  Thus  che  delay  in 

obtaining  SWT  resulcs  for  an  elastic  sphere.  Doolittle  ec  al.  separaced  che 

surface  waves  they  found  on  cylinders  inco  two  groups;  "Franz-Ilka"  waves 

whose  energy  resided  mostly  in  che  liquid  around  die  cylinder  and  which  were 

present  in  che  rigid  and  sofc  cylinder  cases,  "Rayleigh- like"  waves  whose 

energy  resided  mainly  in  che  elastic  cylinder.  Investigations  carried  ouc 

over  che  next  five  years  showed  chat  che  "Rayleigh- like"  surface  waves  could 

be  further  spile  by  Identifying  one  Rayleigh  wave  whose  phase  velocity  wane 

asymptotically  as  ka  30  to  che  phase  velocity  of  che  Rayleigh  surface  wave 

known  to  exist  on  a  flat  interface  between  a  liquid  and  a  solid  (the  possible 

existence  of  a  unique  type  of  surface  wave  along  che  flat  interface  between  a 

14 

vacuum  and  a  solid  was  first  demonstrated  100  years  ago  by  Rayleigh  in  1885 

hence  Che  name)  and  a  number  of  "whispering  gallery”  surface  waves  whose  phase 

velocities  went  asymptotically  co  the  phase  velocities  of  che  3hear  or 

12 

longitudinal  waves  of  che  solid.  These  whispering  gallery  waves  correspond 
in  the  limit  that  ka  —  *>  co  shear  or  longitudinal  surface  waves  which  are 
found  to  exist  when  a  plane  wave  is  incident  from  a  liquid  onto  a  flat 
interface  between  che  liquid  and  a  solid  ac  che  shear  critical  angle  6  or 
longitudinal  critical  angle  9  .  These  critical  angles  are  related  co  che 

la 

relative  acoustical  velocities  of  waves  in  the  liquid  and  the  solid. 

Reference  12  discusses  che  above  identifications  in  some  detail. 

Before  proceeding  it  is  appropriate  co  further  characterize  che 
Rayleigh  surface  wave  since  it  will  prove  co  be  a  central  character  in  che 
analysis  which  follows.  The  particle  displacement  for  a  Rayleigh  wave  is  a 
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unique  combination  of  shear  and  longitudinal  displacement.  The  amplitude  of 
the  shear  and  the  longitudinal  particle  displacements  decay  exponentially  as 
one  proceeds  away  from  the  liquid  solid  interface.  The  decay  constants, 
however,  are  different.  The  phase  velocity  of  the  Rayleigh  wave  on  a  flat 
interface  is  slightly  less  than  the  shear  wave  phase  velocity  of  the  solid.  A 
discussion  of  the  flat  interface  Rayleigh  wave  can  be  found  in  Ref.  15. 

Figure  3.8  of  chat  reference  shows  Che  combined  particle  displacements  of  a 
Rayleigh  wave  as  a  function  of  the  observer's  distance  from  the  interface. 

That  figure  has  been  redrawn  here  as  Fig.  4.  In  1969  Rulf  published  an 
article  examining  Rayleigh  waves  on  cylinders. ^  His  work  indicated  chat  the 
Rayleigh  wave  velocity  on  the  surface  of  a  cylinder  is  a  function  of  the  ka  of 
the  cylinder.  This  velocity  went  asymptotically  to  the  flat  surface  Rayleigh 
wave  phase  velocity  as  ka  •*  ®.  Similar  behavior  will  be  found  in  Chapter  2 
for  the  Rayleigh  wave  when  seen  on  spheres. 

Experimentally,  backscattered  echoes  resulting  when  a  short  tone  burst 
was  incident  on  an  aluminum  cylinder  showed  a  decaying  pulse  train  similar  to 
that  shown  before  for  the  Armico  iron  sphere.17  The  SWT  results  allowed  one  to 
identify  some  of  these  echoes  as  due  to  specular  reflection  and  bulk  waves 
transmitted  within  the  cylinder  and  others  as  due  to  repeated 
circumnavigations  of  the  cylinder  by  surfaces  waves  which  periodically 
radiated  energy  in  the  backscattering  direction.17  In  short,  the  SUT  analysis 
for  elastic  cylinders  allowed  an  interpretation  of  the  echo  structure  not 
obvious  by  use  of  the  PWS  alone.  The  same  sort  of  ideas  should  be  applicable 
to  elastic  spheras.  However,  one  might  expect  new  phenomena  due  to  the 
increased  symmetry  of  the  sphere. 

The  SWT  analysis  was  also  found  to  be  useful  in  understanding  steady 
state  results.  In  1977  Uberall,  Dragonette,  and  Flax  published  a  paper  on  the 
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relation  between  surface  waves  (also  known  as  creepiag  waves)  and  normal  modes 
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of  vibration  of  a  curved  body.  Each  term  of  the  PWS  solution  to  a  problem 

can  be  identified  as  a  multipole.  Each  term  may  resonate  (i.e.,  dlvergs)  at 

several  ka  values.  Uberall  et  al.  argued  that  for  a  cylinder  the  resonance  of 

the  nC^  term  of  the  PWS  occurred  at  ka  values  where  n  cycles  of  a  particular 

surface  wave  fit  on  the  body.  These  ideas  had  been  useful  in  particle  physics 

where  it  was  termed  Raggi  pole  analysis.  These  ideas  eventually  led  Flax  to 

propose  in  1978  an  acoustical  and  electromagnetic  scattering  theory  which 

examined  PWS  results  directly  and  used  nuclear  scattering  theory  ideas.  This 

19 

theory  was  termed  Resonance  Scattering  Theory  (RST). 

RST  achieves  a  separation  between  rapidly  varying  " resonances" 

portions  of  the  PWS  for  the  scattered  signal  and  a  slowly  varying  background. 

In  particular,  for  an  acoustical  plane  wave  scattered  from  fluid- loaded 
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elastic  bodies  such  as  cylinders  or  spheres  the  RST  was  used  to  interpret 
the  rapid  fluctuations  in  the  backseat tered  pressure  as  the  ka  of  the  body  is 

varied  (cf.  Fig.  2a).  One  of  the  Important  numerical  tasks  in  RST  is  to  find 

the  complex  "resonance"  ka  values  which  locate  the  poles  of  the  scattering 
amplitude.  Two  basic  ideas  are  then  used  in  interpreting  scattered  pressures. 
For  a  sphere  these  ideas  translate  into  the  following  statements"  :  (a)  each 

of  the  many  resonance  ka's  can  be  labeled  with  the  integers  n  and  l  and  at 
the  (n,£)th  resonance  (n  +  1/2)  wavelengths  of  the  2.th  surface  wave  fit 
onto  the  circumference  of  the  sphere,  (b)  that  any  rapid  fluctuation  in  the 

backacattered  pressure  is  due  to  in-phase  adding  of  a  specific  surface  wave 

[(n,l)  resonance]  and  these  fluctuations  are  thus  associated  in  a  simple 
fashion  with  the  set  of  multipole  resonances  for  the  sphere.  Two  of  the  major 
goals  of  the  last  chapter  of  this  dissertation  are  to  examine  the  limitations 


of  these  ideas. 


An  example  of  the  type  of  results  which  can  be  developed  through  the 
use  of  the  RST  Ideas  is  shown  In  Fig.  5.  In  that  figure  (which  is  equivalent 
to  Fig.  8  of  Ref.  21)  the  absolute  value  of  backscattering  form  function 
(related  to  the  backacattering  pressure,  cf .  Chapter  3)  for  a  tungsten  carbide 
sphere  is  shown.  Also  shown  are  some  of  the  resonance  ka  values.  Their 
positions  are  indicated  by  arrows  and  they  are  labeled  by  the  (n,£)  pairs 
whose  meaning  was  discussed  in  the  last  paragraph.  Rapid  fluctuations  in  the 
form  function  are  found  to  occur  at  ka  values  close  to  specific  (n,£)  pairs. 
In  Chapter  4  we  Investigate,  through  the  use  of  the  SWT,  the  specific  case 
depicted  in  Fig.  5  in  order  to  better  understand  the  caution  which  should  be 
exercised  in  using  the  RST  labeling  of  the  backscattering  form  function. 

It  is  Interesting  to  note  in  passing  the  similarity  between  the 
structure  found  in  Figs.  2a  and  Fig.  5.  It  would  seem  that  appropriate 
rescaling  of  the  ka  axis  of  Fig.  2a  would  yield  results  which  could  be 
overlaid  on  the  tungsten  carbide  result  with  good  agreement.  This  might 
suggest  similar  physical  reasons  for  the  structure. 

It  should  be  emphasized  again  that  in  the  case  of  a  sphere  that  the 
RST  attempts  to  Interpret  the  PWS  results  such  as  shown  in  Fig.  2a  in  terms  of 
the  surface  wave  picture  developed  through  the  SWT  without  ever  doing  the  SWT. 
Once  the  reader  has  examined  the  SWT  analysis  for  the  fluid- loaded  elastic 
sphere  given  in  Chapter  3  of  this  dissertation  he  will  probably  appreciate 
this  goal  of  RST.  However,  once  he  has  read  the  last  chapter  it  is  hoped  that 
he  appreciates  the  need  for  caution  when  trying  to  reach  this  goal. 

The  reader  should  not  get  the  idea  that  the  era  between  Nussenzveig' s 
o  10  21 

papers  '  of  1969  and  Gaunaurd's  paper  in  1983  was  devoid  of  fruitful 
research  on  acoustical  scattering  from  elastic  spheres.  For  example,  in  the 
seventies  and  early  eighties  other  researchers  used  Fourier  transform 
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The  baekscattaring  fora  function  of  a  tungsten  carbide  spb  vm  as  a 
function  of  ka.  Also  indicated  b j  the  arrows  ars  chs  real  part  of 
the  RST  rasonancas  for  a  -  0-7  and  l  •  1-8.  the  RST  raaonanca 
locations  wars  taken  fro*  Table  IV  of  Ref.  21  and  the  figure  is 
siallar  to  Pig.  4  of  the  saae  reference. 


Fig.  5 . 
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techniques  to  obtain  steady  state  results  analogous  to  Fig >5  from  transient 
22-24  . 

experimental  data.  However,  we  neglect  discussing  these  and  other  works 

in  any  detail  since  they  are  not  of  direct  relevance  to  the  development  of 
this  dissertation. 

Further  discussion  of  RET  as  applied  to  spheres  can  he  found  in  Refs. 

25  and  26.  These  references  along  with  Ref.  21  were  the  latest  research 

results  available  on  acoustical  scattering  from  elastic  spheres  with 

0  <  ka  <  100  when  the  research  in  this  dissertation  was  initiated  in  1983. 

At  this  point  it  is  appropriate  to  discuss  the  avenue  through  which 

this  dissertation  topic  waa  arrived  at  since  it  influenced  greatly  the 

direction  of  the  research.  The  phenomena  of  acoustical  glory  backscatterlng 

(which  was  mentioned  earlier  in  connection  with  Nussenzveig' s  papers)  had  been 
27-29 

under  study  JLn  the  yearn  preceding  1983.  The  spheres  being  studied 

experimentally  and  theoretically  had  ka  values  typically  much  larger  than  100. 
At  these  ka  values  the  radiation  damping  of  surface  waves  has  for  the  most 
part  caused  them  to  be  of  little  importance  to  backscatterlng.  The  most 
Important  waves  are  those  of  geometrical  acoustics;  i.e.,  reflected  and 
transmitted  waves  whose  prograss  can  be  ascertained  by  ray  tracing  techniques. 
In  the  research  of  Ref.  27-29  Che  term  glory  rays  was  used  to  describe  rays 
whose  corresponding  waves  were  weakly  focused  along  the  backscatterlng  axis; 
i.e.,  glory  scattering.  A  few  glory  rays  are  shown  in  Fig.  6  along  with  a 
couple  of  rays  whoee  waves  are  not  focused.  The  rays  labeled  by  n  ■  3,4,5 
are  glory  rays.  The  value  of  n  gives  the  number  of  ray  chords  within  the 
sphere.  For  instance  the  ray  with  n  -  3  is  transmitted  within  the  sphere 
then  internally  reflected  twice  before  being  tranomitted  along  the 
backacattering  axis  which  gives  three  chords  within  the  sphere.  The  rays 
labeled  0  and  2  are  not  glory  rays.  Ray  0  is  the  specular  reflection  from  the 


Savaral  rays  tracing  cranaaittad  wmva  propagation  within  cha  sphara 
ara  shown.  Tha  paranatar  a  indicacts  tha  aumbar  of  ray  chorda 
within  tha  sphara.  Tha  rays  labalad  a  *  3, 4, .5  ara  assodatad  with 
wavaa  waakly  focuaad  along  tha  backacattaring  axis;  l.a.,  chay  ara 


glory  rays.  Tha  wavaa  assoclatad  with  rays  labalad  n  •  0,  2  ara  not 


focuaad 


front  of  eh«  sphere  vhila  Ray  2  is  Cha  ray  which  Is  transmitted  within  tha 

sphara,  reflected  from  tha  back  of  cha  sphere,  and  than  ratranaaltted  along 

tha  backscattaring  axis.  To  proparly  modal  cha  waak  axial  focusing  of  rays 
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such  as  3,  4,  and  5  ona  must  usa  diffraction  tachnlques.  ’  Ona  finds  chat 
glory  scattaring  is  characterized  by  a  JQ(kby)  angular  dapandance  in 
prassure  amplitude  for  small  angles  y  away  from  cha  backscattaring  axis 
where  JQ  is  cha  zeroth  ordar  Bessel  function.  Tha  parameter  b  is  shown  in 
Fig.  6  for  tha  n  ■  3  glory  ray  and  its  significance  is  discussed  in  detail 
in  Refs.  27-29  as  wall  as  Chapter  2.  Tha  JQ  angular  dapandance  was 
confirmed  experimentally  by  usa  of  backscattaring  echoes  resulting  from  Cone 
bursts  incident  on  a  fused  silica  sphere  with  ka  =>  457.  Because  of  cha 
different  propagation  times  Che  cone  burst  experiments  allowed  isolation  in 
time  of  various  glory  ray  contributions.  A  comparison  between  the  theoretical 
and  experimental  angular  dependence  for  a  glory  ray  with  n  ■  4  is  shown  in 
Fig.  7.  The  model  also  allowed  calculation  of  Che  absolute  pressure  amplitude 
along  Che  backscattaring  axis  (when  y  -  0)  for  any  glory  ray.  These  model 
predictions  were  tested  experimentally  for  several  glory  rays  and  were  found 
to  give  answers  within  5  percent  of  experimental  results. 

It  is  interesting  to  ascertain  whether  surface  wave  contributions 
which  are  important  at  lower  ka's  can  be  shown  experimentally  to  display  weak 
axial  focusing.  Next,  if  surface  waves  do  display  glory  scattering  can  one 
develop  a  model  of  surface  wave  glory  which  allows  one  to  predict  Che  absolute 
pressure  amplitude  along  che  backscattering  axis? 

AC  about  Che  same  time  chat  the  last  reference  on  transmitted  wave 
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glory  was  submitted  tha  RST  article  by  Gaunaurd  was  published.  It  was 
felt  that  research  which  would  be  necessary  to  answer  Che  questions  above  on 
surface  wave  glory  might  also  be  useful  in  developing  a  better  understanding 


■I  0  I 

maaavess*  oa*-ACtM«Nr  •  or  htowothok  «»i 


Thn  rasule  of  a  comparison  of  cha  sxparimancal  angular  dapandanca 
found  aa  a  racaivar  hydrophona  was  diaplacad  away  from  cha 
backacaetaring  axis  vs.  cha  chaoraelcal  JQC’*bY)  dapandanca 
pradicead.  Tha  docs  ara  oxparlmaucal  rasulcs  and  cha  curva 
Chaoraelcal  rasulcs.  tha  glory  ray  balng  eaatad  in  chls  casa  had 
a  »  4.  Along  cha  horizontal  axis  x  -  0  corraaponds  Co  Y  «  0  and 
(x|  -  2  corraaponds  eo  Y  3  >055  radians. 


Fig.  7 
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of  RST.  In  cha  course  of  answering  cha  quaaciona  on  surface  wave  glory 
scattering  and  crying  co  undaratand  RST  cha  hiscorlcal  perspective  which  haa 
formed  a  larga  pare  of  Chis  lncroduccion  was  davalopad. 

Tha  next  thraa  chapters  addraas  cha  lasuas  puc  forch  above.  Soma 
commanes  on  cha  ganaral  scruceura  of  aach  chapear  ara  useful.  Each  chapear  la 
in  raallcy  a  paper  which  haa  baan  aubmiccod  co  cha  Journal  of  cha  Acouscical 
Society  of  America.  As  such,  aach  has  iCs  own  incroductlon  section, 
popendices,  rafarancaa,  figura  numbers,  and  equation  numbers.  The  papers 
which  form  Chapters  2  and  3  have  baen  accepted  for  publication  at  cha  present 
time  while  cha  paper  which  forces  Chapter  4  has  only  recently  been  submitted. 
Also  included  as  a  separate  appendix  ara  Che  computer  programs  which  ware 
developed  for  the  calculations  carried  out  in  Chapters  2  through  4.  This 
appendix  contains  text  which  Introduces  each  program  and  briefly  discusses  any 
cautions  deemed  appropriate. 

In  Chapter  2  a  theoretical  model  for  surface  wave  glory  is  developed 
via  cwo  different  routes.  Tha  first  uses  GTD,  RST,  and  some  of  cha  ideas  and 
results  used  in  che  transmitted  wave  glory  model.  This  method  allows  a  simple 
intuitive  picture  co  be  formed  of  Che  production  of  surface  wave  glory.  The 
second  rneehod  uses  a  slightly  different  form  of  Che  (PWS)  in  Eq.  (9).  The  use 
of  the  SWT  is  intentionally  avoided  in  Chapter  2.  However,  it  is  through  the 
use  of  Che  SWT  (which  is  carried  out  in  Chapter  3)  that  one  can  obtain  a 
theoretical  prediction  of  the  on- axis  pressure  amplitude  of  the  surface  wave 
glory.  In  Chapter  2  Che  model  developed  for  angular  dependence  of  the  surface 
wave  glory  is  experimentally  tested  using  tungsten  carbide  spheres  with 
30  <  tea  <  100.  The  agreement  between  experiment  and  theory  is  good. 

In  Chapter  3  the  SWT  is  performed  and  the  relation  between  SWT  results 
and  those  of  Chapter  2  are  discussed.  The  SWT  of  Chapter  3  concentrates 


asialy  on  contributions  from  the  specular  reflection  end  surface  waves  at 
small  backscattering  angles.  The  SWT  Is  a  more  rigorous  method  through  which 
one  can  predict  the  possible  experimental  observation  of  glory.  It  confirms 
the  physical  model  used  in  Chapter  2  and  for  the  first  time  predicts  the 
absolute  backscattering  amplitude  associated  with  one  or  more 
circumnavigations  of  the  surface  wavs  around  the  sphere;  i.e.,  it  predicts  the 
amplitude  of  various  surface  wave  contributions  to  the  decaying  pulse  train 
seen  in  tone  burst  experiments.  This  prediction  is  tested  experimentally  in  ' 
Chapter  3  using  tungsten  carbide  spheres.  Again,  the  agreement  is  good. 

Finally  in  Chapter  4  the  RST  Ideas  discussed  previously  In  this 
introduction  are  reexamined  using  the  SWT  results  of  Chapter  3.  The  chapter 
uses  the  specific  case  of  the  backscattering  form  function  of  a  tungsten 
carbide  sphere  (cf.  Fig.  5).  The  main  effect  of  this  analysis  is  on  the 
interpretation  of  RST.  The  results  show  that  care  is  needed  when  using  the 
labeling  of  Fig.  5  to  understand  fluctuations  in  the  backscattering  form 
function.  It  is  also  3hown  that  RST  ka  resonances  only  approximate  the  ka 
values  at  which  surface  waves  add-in  phase.  However,  the  error  in  this 
approximation  i3  usually  negligible. 
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2.1  Introduction 

An  understanding  of  the  echo  structure  and  angular  dependences  for 

scattering  from  spheres  (or  other  elastic  objects  ensonlfled  along  an  axis  of 

symmetry)  may  be  facilitated  by  modeling  the  focusing  of  sound  along  the 

backward  axis.  We  have  previously  studied  this  axially  focused  scattering 
12  3 

from  large  elastic  *  or  fluid  spheres  where  the  focused  waves  were 
associated  with  rays  transmitted  through  (and  Internally  reflected  from)  the 
sphere's  surface.  The  cases  studied  had  ka  >_  100.  In  the  present  paper  we 
demonstrate  axially- focused  scattering  from  tungsten  carbide  spheres  in  water 
having  ka  as  low  as  30.  For  this  choice  of  material  and  low  range  of  ka  (30 
<  ka  <  100),  it  is  shown  that  surface  waves  on  the  sphere  which  radiate  sound 
back  into  the  water  produce  axially  focused  backscattarlng.  The  specific  type 
of  surface  wave  associated  with  ehe  echoes  studied  Is  the  one  usually  referred 

4 

to  as  the  Rayleigh  wave;  however,  our  method  of  analysis  should  be  useful  for 
modeling  the  angular  dependence  of  echoes  associated  with  other  waves  such  as 

4 

“whispering  gallery"  waves.  Indeed,  the  form  of  the  angular  dependence 
demonstrated  for  the  surface  wave  contribution  to  the  pressure,  as  a  function 
of  backseat taring  angle,  should  be  applicable  to  spherical  elastic  shells  and 
to  certain  other  shapes  ( such  as  spheroids  or  cylinders  having  hemispherical 
caps)  when  the  object  is  ensonlfled  along  the  axis  of  symmetry.  The  angular 
iependence  found  differs  from  that  expected  for  surface  wave  contributions  to 
the  backscattarlng  from  right  circular  cylinders (for  which  the  direction 
of  incident  sound  lies  perpendicular  to  the  axis  of  symmetry) . 

The  backscattarlng  echo  in  Fig.  1  was  produced  by  a  four  cycle 
sinusoidal  tone  burst  incident  on  a  tungsten  carbide  sphere  with  ka  =  49.1. 

The  dominant  contributions  to  the  signal,  after  the  specular  reflection,  are 
due  to  Rayleigh  surface  waves.  We  experimentally  and  theoretically 


An  oacilloscopa  erne*  of  chm  backscactaring  acho  from  a-  cungscan 
carbida  sphara  ansooifiad  by  a  cona  burse.  Tha  ka  of  cha  iphara  Is 
approxiaaealy  49. L.  Tha  Individual  achoas  ara  Labalad  as  follows:  A 
-  spaeular  caflaccion,  B  -  Raylaigh  surfaca  wavs  achoas,  C  -  achoaa 
who a a  spacific  origin  waa  one  identifiad  buc  which  ara  conjacturad  eo 
ba  ochar  surfaca  wavaa  such  aa  “Whispering  Gallery"  waves. 
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investigated  cha  angular  dependence  of  these  Rayleigh  achoas  as  ona  novas  away 

from  the  backscattering  direction.  Our  chaoratical  analysis  uses  the  high 

frequency  mat hods  of  tha  Geometrical  theory  of  Diffraction  (GTD)^  and 

4 

Resonance  Scattering  Theory  (RST).  The  angular  dapandanca  found  can  be 

approximately  modeled  by  Jg(gy)  where  8  depends  on  tha  frequency  of  tha 

tone  burst  (or  equivalently  tha  ka  of  tha  sphere)  and  y  is  Che  angle 

relative  to  backscattering .  This  angular  dependence  Is  due  to  axial  focusing 

1-3 

and  is  characteristic  of  what  has  been  called  acoustical  "glory.”  The  term 

"glory”  comes  from  the  enhancement  of  backscattered  light  from  cloud  droplets 

3  8  9—12 

due  to  axial  focusing.  *  Previous  research  on  tungsten  carbide  spheres 

does  not  directly  address  the  affects  of  axial  focusing. 

In  Section  I  we  present  a  theoretical  model  of  the  Jq(8Y)  angular 

dependence  of  Rayleigh  echoes  via  two  different  routes.  The  first  uses  GTD 

and  RST  together  with  results  of  our  previous  analysis1  ^  of  virtual  ring-like 

sources.  This  method  allows  a  simple  intuitive  picture  to  be  formed  of  the 

production  of  Rayleigh  echoes.  The  second  method  uses  the  Partial-Wave  Series 

(PWS )  resulting  from  solving  the  steady-state  scattering  from  a  sphere.  The 

approximate  dependence,  Jq(Sy),  also  follows  from  the  form  of  the  PWS  for 

spheres  by  way  of  a  Sommerf eld-Watson  Transformation  (SWT)  or  a  “modified"  SWT 

(see  e.g.  Eq.  (5.17)  of  Ref.  8).  Our  analysis  given  in  Sec.  I. A  has  the 

advantage  that  it  relates  8  to  the  phase-velocity  of  the  Rayleigh  wave 

without  further  use  of  sophisticated  analytical  methods.  In  addition,  our 

analysis  may  be  extended  to  nonepherical  objects  (subject  to  restrictions 

noted  in  Sec.  Ill)  for  which  the  PWS  may  not  be  known.  In  Sec.  IC,  the  model 

is  specialized  to  the  particular  case  of  tungsten  carbide  spheres. 

We  intentionally  avoid  the  uae  of  the  SWT  in  Sec.  I.  However,  it  is 

through  the  use  of  the  SWT  that  the  ka  dependence  of  the  Rayleigh  echo  at 
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fixed  y  may  b«  determined.  W«  return  co  cha  SWT  analysis  In  a  subsequanc 

papar  slnca  cha  analys.-*  rcr  elastic  spharaa  diffars  in  significant  datails 

8  13 

from  cha  analysis  for  othar  spharaa  or  alasclc  cylindars. 

In  Saction  II  we  east  Cha  angular  dapandanca  modal  on  tungstan  carbida 

spharaa  with  30  <  ka  <  100.  Tha  axpariaancal  setup  is  similar  Co  thac  used 

1  2 

previously  in  transmitted  wave  studies  on  alasclc  spheres,  *  and  so  will  be 

discussed  only  briefly.  It  is  appropriate  to  comment  on  our  use  of  tone 

bursts- for  comparison  with  tha  steady-state  analysis  given  in  See.  I.  This 

1  14 

issue  was  addrassed  in  earlier  work  *  on  axially- focused  transmitted  waves. 
It  was  found  chat  chough  cha  scattered  tone  bursts  should  be  distorted,  the 
peak-co-peak  amplitudes  of  the  central  cycle  of  a  four-cycle  burst  should  be 
given  by  a  steady  state  analysis.  It  is  anticipated  chat  a  detailed  analysis 
of  the  shapes  of  transient  Rayleigh  echoes  would  lead  to  similar  conclusions. 

Tungsten  carbide  was  particularly  well  suited  as  a  target  material  for 
these  experiments  because  (i)  the  Raylelgh-echo-amplltudes  were  significant 
for  the  range  of  ka  examined  and  (11)  transmitted  bulk-wave  echoes  were  either 
insignificant  in  amplitude  or  well  spaced  in  time  from  the  Rayleigh  echoes  of 
interest.  These  properties  are  in  part  a  consequence  of  tungsten  carbide's 

3 

large  density  -  14  g/cm  . 
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2.2  Modal  of  angular  dapendence 


using  GTD ,  RST .  and  Virtual  Ring-Lika  Sources 


Yha  foundation  of  GTD  was  laid  by  Keller. Tha  GTD  is  usaful  in 
casas  vhara  tha  scattaring  body  fulfills  tha  condition  kd  ?~>  1  whara  d  is 
a  charactarlstic  dimansion  of  tha  scattersr.  It  is  baaed  on  the  postulate 
that  fialds  propagate  along  rays,  but  Introduces  diffracted  rays.  In  the  case 
of  Rayleigh  surface  waves  on  a  sphere  it  allows  us  to  form  the  picture  In  Fig. 
2.  Tha  "ray"  AA'  gats  onto  tha  sphere  at  angle  8  .  It  proceeds  around  the 

A 

sphere,  along  its  surface,  shadding  energy  continuously  back  into  the  water  at 
angle  8^  as  indicated  by  the  dotted  lines.  Of  particular  interest  is  the 
fact  that  at  ana  point  <B')  anargy  is  shad  along  B'A'.  This  process  continues 
and  after  circumnavigation  of  tha  sphere  again  more  energy  is  shed  in  tha 
backward  direction.  This  process  can  ba  pictured  for  any  number  of  trips 
around  tha  sphere.  A  smaller  amplitude  echo  will  be  seen  after  each  trip 
because  of  tha  anargy  shadding,  i.e.  radiation  damping. 

The  angle  9  in  Fig.  2  can  ba  determined  by  the  following  phase 

A  ,  * 

matching  criteria.  For  continual  reinforcement  of  the  Rayleigh  wave  the  pi  ise 

velocity  of  the  wave  c^  must  match  the  component  of  tha  phase  velocity  in 

water  c  along  tha  sphere's  surface.  This  component  is  c/sind  .  Thus  we 

can  determine  9  through  the  equation 
ft 


sin  9t 


which  is  actually  Just  Snell's  law. 

Tha  Above  ideas  explain  qualitatively  Fig.  1  though,  as  yet,  they  give 
no  clue  as  to  the  angular  dependence  that  will  be  seen  as  one  moves  away  from 


Surface  wave  ’ray"  diagram  ualng  Che  Geometric  theory  of  Diffraction 
(GTD) .  An  incoming  plana  wave  represented  by  AB  allows  the  launching 
of  a  surface  wave  at  3  which  circumnavigates  the  sphere  while 
reradlacing  back  into  the  surroundings*  At  point  B '  energy  is 
radiated  in  Che  backward  direction.  This  accounts  for  the  echoes 
labeled  3  in  Fig.  1.  5^  is  cha  virtual  point  source  from  which  che 

ray  A’B'  and  cha  dashed  rays  co  either  side  appear  co  originate. 

When  cha  diagram  is  rotated  around  che  C'C  axis,  cha  point  F  cracas 

K 

out  a  virtual  ring- like  source. 


ch«  backscattarlng  axis.  To  obtain  the  angular  dapandanca  of  Cha  Rayleigh 
achoas  ona  must  axamina  tha  daahad  linas  to  aithar  side  of  B  *  A ' .  Thasa 
rapreaant  radiation  damping  slightly  away  from  tha  baclcacattering  direction. 
By  tracing  thasa  raya  as  wall  as  B'A'  backward  wa  arrive  at  7  from  which 

n 

they  appear  to  originate.  (In  tha  Appendix  tha  horizontal  position  of  F_ 

R 

behind  the  vertical  Una  through  C'  in  Fig.  2  is  found  to  ba  equal  to  a, 
the  sphere's  radius.)  Next,  because  of  tha  spherical  symmetry  of  tha 
situation,  ona  can  rotate  tha  picture  around  tha  line  CC'  in  which  casa  tha 
point  F^  traces  out  a  circle.  Thus  one  seas  that  the  Rayleigh  waves  appear 
to  emanate  from  a  ring- like  source  with  radius  b  -  a  sin8^.  When  extended, 
the  F  D  ray  intersects  tha  backscattarlng  axis.  Rays  such  as  F  D  are 

iv  R 

responsible  for  the  axial  focusing  whan  the  CC*  rotation  is  performed.  (In 

the  limit  that  the  observation  point  goes  to  infinity  F  D  goes  to  FA' .) 

R  R 

Ring-like  sources  have  been  discussed  in  tha  previous  work  on  acoustical  glory 
due  to  transmitted  wavas  in  spheres.1  ^  These  sources  lead  to  an  angular 
dependence  for  pressure  amplitude  which  can  be  approximated  as 

P(Y)  -  p(Y-O) | JQ(kbY) |  ,  y  «  1  radian  (2) 

where  Y  is  the  angular  position  of  the  observation  point  in  Che  Fraunhofer 

or  Fresnel  region,1  relative  to  the  backscattarlng  axis  and  Che  sphere's 

center.  (See  Eq.  (16)  of  Ref.  1  where  in  the  present  case  z  *  z  +  a, 

n 

Y  3  tanY  •  h/zn,  an<*  h  corresponds  to  x  in  Fig.  4  of  the  present  paper.) 

Tha  results  to  this  point  lead  one  to  expect  Che  following  form,  in 
the  backscattering  region,  for  the  pressure  amplitude  of  individual  Rayleigh 
echoes  in  Fig .  2 

Pm  "  Am(^)  V^Y)  (3) 


1 "* ,  ^ ■  *  '***’<.'•*<.  ^ "I*.  -"o".  ’*■  ■' *  "*«  •"‘a,  '*•  "* 
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where  the  argxiaanc  of  the  Bessel  function  combines  Eqs.  (1)  and  (2),  and 
8  -  kac/cR,  and  n  ia  a  natural  number  (0,1,2  .  .  .)  .  In  thia  axpraaaion  m 
indexes  which  Raylaigh  echo  we  are  examining,  l.e.  a  »  0  impliea  the  first 
Rayleigh  echo  after  the  specular  reflection  etc.  A^  is  the  on-axis  amplitude 
of  a  particular  Rayleigh  echo  and  aay  be  a  function  of  ka.  The  phase 
velocity  ratio  may  also  be  a  function  of  ka. 

Calculation  of  the  on-axis  amplitudes  A  (ka)  and  the  phases  of 
Rayleigh  echoes  requires  use  of  the  SWT.  However,  knowledge  of  this  parameter 
is  not  required  in  order  to  teat  the  JQ  dependence  of  the  Rayleigh  echo.  We 
therefore  defer  to  a  subsequent  paper  calculations  of  the  A^ ' s  and  of  the 
phases  of  the  Raylaigh  echoes. 

To  obtain  angular  dependence  results  to  compare  with  experiment  one 

must  know  the  phase  velocity  ratio  c/cR  as  a  function  of  ka.  This 

dependence  can  be  found  through  the  results  of  RST.  The  form  function  of  an 

9 

elastic  sphere  is  given  by 

00  3  (ka) 

f  -  (2/ lka)  l  (2n+l)  r^-r—  P  (cose)  (4) 

a-0  Da(ka)  n 

where  9  denotes  the  scattering  angle  relative  to  forward  scattering  and  B_ 

and  D  are  determinants  of  3  x  3  matrices  whose  elements  are  combinations  of 
n 

9 

spherical  Bessel  and  Hankel  functions.  In  RST  one  examines  the  PWS  of  Eq. 

(4)  term  by  term.  Resonances  in  each  term  are  determined  numerically.  One 
finds  the  resonant  ka  values  for  the  nth  term  by  solving  the  equation 

W  '  0  (5> 

where  is  the  1th  resonance  of  the  nth  term.  Knowing  the  value  of  a 

particular  resonance  in  the  nth  term  allows  one  to  approximate  the  phase 


v.v.y 


1 *i  v 


* 


velocity  of  the  surface  wavs  corresponding  to  chat  resonance.  The  ralaclon 

A  17 

uaad  co  obtain  Cha  phase  velocity  at  the  resonance  values  X  ,  is  * 


cp/c  a  Re(Xaa)/(n  +  1/2) 


The  Rayleigh  phase  velocity  cR  is  given  by  c?  when  l  •  1.  For  the 
specific  case  of  a  tungsten  carbide  sphere  the  resonances  corresponding  co 

Q 

Rayleigh  waves  have  been  found  for  n  *  0  to  7;  c  /c  as  a  continuous  function 

ll 

of  ka  can  be  found  by  curve  fitting  Che  discrara  values  found  via  Eq.  (6). 
Equation  (3)  with  cR(ka)/c  found  using  Eq.  (6)  represents  Che  final  results 
for  Che  angular  dependence  predicted  for  Rayleigh  wave  echoes.  It  is  evident 
that  Che  reasoning  above  can  also  be  applied  to  whispering  gallery  waves. 

B.  Angular  Dependence  using  Legendre  polynomial  Expansions 
It  is  appropriate  to  discuss  briefly  an  alternate  route  to  the  angular 
dependence  found.  If  the  forrr  function  of  Eq.  (4)  is  rewritten  using  y  ■  tt  - 
9,  where  Y  is  Che  angle  relative  to  backscattering  and  here  9  is  the 
scattering  angle  relative  to  the  forward  direction,  it  becomes 


(2/ika)  l  (-l)a(2a+l) 
n*0 


3  (ka) 

—-7.-7  P  (cosy) 
D  (ka)  n 
a 


8  18 

For  small  angles  Y  one  can  use  the  following  expansion  derived  by  Szego  ’ 

Pn(cosY)  *  (Y/sinY)1/2  JQ((n  +  l/2)y)  +  0(n"  3/2)  (8) 

4 

Close  to  a  Rayleigh  resonance  ka  s  Re  X  ,  and  the  Y  dependence  of  the 

n,  l 

nth  partial  wave  amplitude  is  similar  to  Chat  of  Eq.  (3).  This  follows  by 
noting  Eq.  (6)  gives  c_/c  3  ka/(afl/2)  and  hence  chat  S  -  kac/c_  3  n+1/2. 

K  R 

Taking  y/slny  =  1  reduces  the  y  dependence  of  Eq.  (8)  to  the  form  of  Eq.  (3), 
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Pq(cO«Y)  »  Jq(8Y) 


(9) 


C.  Specializing  Angular  Dependence  Model  eo  a  Tungsten  Carbide  Sphere 


Earlier  it  was  stated  that  the  Rayleigh  resonance  tea' s  were  known  for 


a  tungsten  carbide  sphere  for  the  specific  cases  of  n  «  0  through  7. 


Unfortunately,  for  use  in  comparison  with  our  experiment  we  need  the  Rayleigh 
resonant  tea's  for  a  ■  10  through  30.  Since  these  were  not  available,  we 
developed  a  computer  program  which  would  find  them.  (The  program  was  verified 
with  the  results  given  in  Ref.  9  for  n  ■  4  through  7.)  The  Rayleigh 
resonances  for  a  ■  0  through  30  are  given  in  Table  I.  The  elastic 
parameters  used  in  the  program  are  given  in  the  caption  to  the  table.  The 
results  for  n  *  0  through  3  are  from  Ref.  9,  the  rest  of  the  resonances  are 
from  our  program.  Using  these  resonant  ka's  and  Zq.  (6)  we  obtained  the  phase 
velocity  ratio  as  a  continuous  function  of  ka  for  ka  from  0  to  85.  We  also 
calculated  the  group  velocity  ratio  for  the  Rayleigh  wave  using4'17 


cg/c  -  d(Re(Xnl))/dn  (iq) 

Curves  for  both  the  phase  and  group  velocities  are  given  in  Fig.  3.  Note  that 
the  group  velocity  ratio  is  essentially  constant  for  ka  >_  30.  We  can  use 
this  result  and  integrate  Zq.  (10)  obtaining 

(cg/c)n+G’  -  Re(Xal) 

where  G'  is  an  integration  constant.  Then  writing  G'  as  G'  ■  G  +  cg/(2c) 

we  have 

(cg/c)(n  +  1/2)  +  G  - 


Re(Xnl)  =  ka 


Table  I.  Roots  of  -  0  corresponding  to  the  Rayleigh  wave  resonances  of  a  tungsten  carbide  sphere. 
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for  n  =  A  to  30  have  uncertainties  of  10.0005 


Rayleigh  velocity  ratios  for  a  tungsten  carbide  sphere  as  a  function 
of  lea .  The  individual  curves  are:  A  -  group  velocity  ratio,  B  - 
phase  velocity  ratio  using  Zq.  (11),  C  -  phase  velocity  racio  using 
Eq.  (6),  0  -  Che  phase  velocity  ratio  for  a  Rayleigh  wave  propagating 
along  a  tungsten  carbide-water  plane  interface. 
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Finally,  using  Eq.  (6)  and  soma  algebraic  manipulation  we  have  Che  relation 

c^/c  -  (c  /c)(l  -  G/ka)"1  (11) 

between  the  phase  velocity  ratio  and  Che  group  velocity  ratio.  G  in  thia 

expression  is  a  dimensionless  constant.  We  emphasize  that  this  is  valid  only 

if  c  /c  is  independent  of  lea. 
g 

A  comparison  between  the  exact  curve  for  c  /c  and  the  approximate 

£ v 

formula  of  Eq.  (11)  is  shown  in  Fig.  3.  For  Che  case  represented  in  Fig.  3, 

c  /c  3  2.60  for  ka  >  30  and  G  was  chosen  as  3.1.  For  ka  >  20  agreement 
S  ~  *“ 

between  Eq.  (11)  and  Che  exact  curve  is  good. 

Also  shown  in  Fig.  3  is  the  phase  velocity  ratio  for  a  Rayleigh  wave 
propagating  along  a  flat  tungsten  carbide-water  Interface.  This  ratio  was 
obtained  by  numerically  solving  a  secular  equation  (below)  for  c^,  the 

? 

complex  phase  velocity  along  a  flat  solid- liquid  boundary.  Let  u  -  (c^/c^)” 

2 

and  v  m  (c^/c^)  where  c^  and  c^  are  the  transverse  and  longitudinal 

1 9  20 

phase  velocities  respectively.  The  secular  equation''  ’  is 

(l-2u)2  -  4u2(u-l)I/2(u-v)L/2  -  -i(Pz/Ps)(u-v)1/2[(cT/c)2  -  u]“  1/2 

where  0a  and  P  are  the  densities  of  the  liquid  and  solid  respectively. 

*•  s 

For  the  material  parameters  given  in  Table  1,  this  gives  c^  ■  (3.8154  - 
10.0192)  x  10^  cm/sec  for  the  Rayleigh  root.  The  flat  surface  ratio  shown  in 
Fig.  3  is  Re(c^/c)  3  2.60  and  both  c^/c  and  c^/c  were  found  to  approach 
this  value  as  ka  ■*  «. 

Using  Eq.  (11)  and  0  -  kac/c^  one  has 


S 


(c/c_) (ka  -  G) 


(12) 
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which  can  b«  substituted  into  Eq .  (3)  to  give  the  normalized  pressure 
amplitude  of  the  mch  Rayleigh  echo  when  c^/c  1*  constant. 

In  the  experiment  to  be  described  p(y)  is  determined  as  the  receiver 
Is  moved  away  from  the  backscattaring  axis.  These  data  were  fitted  to  give 
the  location  of  the  first  pressure  null.  This  procedure  is  repeated  for 
several  ka  values.  Examining  Eq.  (3)  one  sees  that  the  pressure  null  location 
corresponds  to  the  point  where  the  Jq  angular  dependence  first  equals  zero, 
l.e.  whan 

(c/c  )<ka  -  G)Ym  -  2.4048  (13) 

Y  is  the  angle  in  radians  corresponding  to  the  pressure  null.  Equation  (13) 
can  be  rewritten  as 

ka  -  (c  /c)  U^OAajY”1  +•  G  (14) 

Since  Cg/c  and  G  are  constants  this  relation  implies  that  a  plot  of  ka 
vs  Y  should  be  a  straight  line.  Furthermore  the  slope  of  that  line 
depends  on  the  value  of  c^/ c.  The  theoretical  prediction  above  will  be 
tested  experimentally  in  the  next  section. 

2.3  Experimental  test  of  angular  dependence 

In  this  section  the  angular  dependence  model,  as  specialized  in 

Sec.  I.C  to  tungsten  carbide,  will  be  tested.  The  experimental  apparatus  has 

1  2 

been  discussed  previously.  ’  Here  we  note  only  the  significant  alterations 
to  the  apparatus  which  facilitate  implementation  of  the  current  experimental 
test,  and  then  proceed  to  Che  experimental  procedure  and  results. 

•  e  • 

Figure  1  of  Reference  1  la  shown  again  in  Fig.  4.  This  tigura  gives  a 
conceptual  description  of  the  experiment.  A  piston-like  transducer  produces  a 
short  sinusoidal  tone  burst  4  cycles  in  duracion.  This  cone  burst  is 
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Simplified  diagram  of  ch«  scat caring  experiment.  tha  hydrophone  may 

be  3 canned  along  a  llna  transverse  to  tha  symmetry  axis  defined  by 

the  source  and  tha  sphere.  y  is  tha  angle  relative  to 

baclcacattering .  Tha  diagram  is  not  drawn  to  scale.  In  all  tha 

experiments  described  in  Section  II,  x  3  160  cm,  while  z  3  33  can 

P 

for  experiments  with  the  2.34  cm  diameter  tungsten  carbide  sphere  and 
z  3  19  cm  for  experiments  witH  the  1.27  cm  sphere. 


?i«.  4 
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scattered  from  a  tungsten  carbida  sphara.  The  resulting  echo  Is  picked  up  by 
a  needle  hydrophone  which  nay  be  positioned  at  small  angles  Y  relative  to 
the  backacattering  axis  (tany  ■  x/(z+a)).  In  tha  present  experiment  the 
Panamatrlcs  model  7309  sending  transducer  of  Ref.  1  was  replaced  by  either  a 
Panamatrlcs  modal  73261  or  a  Panamatrlcs  model  73260.  These  broadband 
transducers  have  lower  resonance  frequencies  than  the  V309  (the  73261  has  a  1 
MHz  resonance  and  tha  73260  a  2.25  MHz  resonance  compared  with  5MHz  for  the 
7309).  In  this  experiment,  which  uses  a  1.27  cm  and  2.54  cm  diameter  tungsten 

carbida  sphere,  these  transducers  allow  investigation  in  the  ka  range  from  30 

to  ICO.  The  only  other  significant  changes  in  apparatus  from  Ref.  1  are  the 
use  of  a  digital  signal  averager  (described  in  Ref.  2)  and  a  new  target 
mounting  system.  The  main  concern  in  designing  a  mounting  system  is  to 
eliminate  spurious  echoes  in  the  time  window  of  the  experiment  ( cf .  Fig .  1 ) . 

The  procedure  used  to  test  the  model  was  to  average  the  echo  received 
at  any  angle  y  over  256  bursts  using  the  digital  scope.  Then  tha  amplitude 
in  volts  of  the  central  cycle  of  the  first  Rayleigh  wave  was  determined.  This 

voltage  and  the  angle  y  were  recorded  in  a  data  file.  The  amplitude  of  the 

first  Rayleigh  wave  was  found  via  the  process  for  at  least  25  values  of  y  on 
either  side  of  Y  ■  0.  The  computer  normalized  all  data  points  by  dividing  by 
the  amplitude  at  Y  ■  0  and  plotted  these  normalized  points  on  an  amplitude 
vs  Y  graph.  The  computer  next  generated  plots  of  the  function  J^CSy)  for 
several  values  of  tha  constant  3  to  find  the  value  yielding  the  closest  fit 
to  the  experimental  results.  An  example  of  the  results  of  this  analysis  is 
shown  in  Fig.  5.  One  immediately  sees  at  least  qualitative  agreement  between 
experiment  and  the  model  of  Eq.  (3). 

The  procedure  above  was  repeated  for  29  frequencies  corresponding  to 
ka  values  from  30  to  100.  The  29  values  of  ka  gave  29  different  8  values. 


45 


These  S  values  ware  usad  co  find  the  smallest  angle  Y_  for  which  Jn(0y  ) 

q  u  m 

equaled  zero,  i.s.  By  •  2*4048.  This  gave  29  values  of  y  which  were 

n  n 

plotted  on  a  graph  of  ka  vs  y 

To  compare  the  experimental  results  for  tea  vs  Y with  the 

theoretical  prediction  of  Eq.  (14)  one  needa  the  group  velocity  ratio  for  the 

sphere  being  uaed.  The  group  velocity  ratio  can  be  obtained  from  the  time 

spacing  between  Rayleigh  waves.  This  spacing  remained  essentially  constant  at 

ail  lea  values  examined  as  it  should  if  the  group  velocity  remains  constant  as 

predicted  in  Fig.  3  for  the  tea  range  of  interest.  The  experlrnenmX  gr..-up 

velocity  ratios  found  were  2.46  for  thii  2.54  cm  sphere  and  2.60  for  the  1.27 

cm  sphere.  One  caveat  should  noted.  The  constant  0  was  found  via  an 

analysis  using  the  elastic  parameters  in  Ref.  9  as  in  Table  I.  The  elastic 

parameters  of  our  sphere  should  be  close  but  not  necessarily  the  same.  This 

difference  will  probably  have  a  small  affect  on  G. 

In  Fig.  6  the  experimental  results  for  ka  vs  Y  ^  are  compared  with 

m 

the  model  prediction  of  Eq.  (14),  The  theoretical  lines  used  5  ■  3.1  and  the 

approximate  experimental  values  of  c  /c.  Figure  6a  shows  the  comparison  for 

8 

the  2.54  ci  sphere  and  Fig.  6b  shows  the  comparison  for  the  1.27  cm  sphere. 

agreement  between  model  aud  experiment  is  good.  This  agreement  confirms 
the  theoretical  model  developed  in  Section  I. 

Figure  7  shows  Che  data  obtained  at  the  lowest  ka  values  examined; 
i.a.,  ka  «  32.3.  The  data  appear  to  have  a  small- amplitude  structure 
superimposed  on  the  Jq(8y)  dependence.  The  cause  of  the  structure  was  not 


determined 


7  ( radians ) 


Experiment  vs  ch«or7  resales  analogous  co  Fig.  5  buc  at  chs  lowest  lea 
value  esscsd;  l.a.  lea  »  32.3.  Tbs  does  art  experimental  rasules. 

Tha  curve  Is  given  by  |Jg(Sy)|  with.  3  adjusted  eo  fit  cha  data. 
Froa  Fig.  6(b)  it  Is  evident  that  this  3  gives  a  value  for  v! 

agrees  with  theory. 


Fig.  7 


2.4  Discussion 


We  have  prassnccd  a  modal  of  cha  axially-focused  scattering  due  to 
surfaca  waves  ganaratad  on  spheres.  Wa  have  varifiad  that  cha  angular 
dapandanca  of  cha  prassura  amplituda  is  3  Jq(By),  which  may  have  been 

Q 

antlcipatad  from  the  SWT  appliad  to  other  surface  wave  scattering  problems. 

Wa  have  related  3  to  the  surfaca  wave's  phase  and  group  velocities  and  have 
verified  tha  modeled  3  for  the  range  of  ka  examined .  The  experiment 
demonstrates  chat  axial  focusing  can  be  Important  down  to  ka's  of  at  least  30. 
These  results  complement  our  earlier  work  on  acoustical  glory  for  large 
spharas  (ka  >_  100)  where  the  axial  focusing  was  due  to  transmitted  waves  not 
surface  waves. ^  ^ 

The  model  of  axial  focusing  was  derived  using  two  methods.  Both 

approaches  led  to  the  same  angular  dependence  prediction;  however,  neither 

method  gave  a  quantitative  prediction  of  the  on-axis  amplitude  of  the  surface 

wave  echoes.  This  shortcoming  can  be  eliminated  via  the  use  of  the  SWT  as 

will  be  shown  in  a  subsequent  paper. 

As  an  indication  of  the  difference  between  the  angular  dependence  of  a 

focused  and  unfocused  echo,  data  on  the  y  dependence  of  the  unfocused 
3 

specular  reflection  from  the  2.54  cm  diameter  tungsten  carbide  sphere  were 
obtained.  These  data  are  shown  as  squares  in  Fig.  5.  Only  small  fluctuations 
In  the  amplitude  of  the  specular  reflection  are  seen.  This  can  be  contrasted 
with  the  amplitude  variation  of  the  focused  Rayleigh  wave.  (The  origin 

of  the  small  fluctuation  in  the  specular  reflection  has  not  been  determined.) 

The  arguments  used  here  in  our  method  of  predicting  the  high-frequency 
angular  dependence  of  surface  wave  echoes  can  be  extended  to  non— spherical 
elastic  objects  possessing  a  symmetry  axis.  When  the  object  is  ensonified 
along  its  axis  of  rotational  symmetry,  axially- focused  surface  wave  echoes 


should  b«  present  in  Cha  naar  backward  scattering  dlraction.  It  is  necessary 
that  tha  phaaa  matching  crir.arla  ba  mat  for  tha  excitation  of  a  surfaca  wava 
and  that  tha  attanuation  asaociatad  with  circumnavigation  ba  sufficiently 
small.  Tha  consequences  of  focusing  ara  avidant  in  pravioua  maasuramants  of 
tha  backscattarlng  of  tona  bursts  from  a  spharoid  in  air  (Fig.  6  of  Raf.  10). 
Tha  cramping  wava  (i.a.,  Franz)  acho  was  significantly  largar  whan  cha 
spharoid  was  ansonifiad  along  its  symmetry  axis  than  whan  ansonifiad 
parpandlcular  to  tha  symnatry  axis. 
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Appendix;  Location  of  the  Focal  Circle 


In  this  appendix  we  calculata  eh*  horizontal  diatanc*  a  of  eh* 
point  in  Fig.  2  behind  eh*  vertical  line  through  C'.  In  Fig.  A1  a 

portion  of  Fig.  2  la  redrawn.  The  aurface  wave  la  launched  at  point  B. 
Radiation  in  the  backacattarlng  direction  la  repreaented  by  ray  A'B'  while 
radiation  at  an  angle  Y  relative  to  backacattarlng  la  repreaented  by  A"B” . 
The  intersection  of  the  backward  extrapolated  A"B”  and  A'B*  in  the  limit 
Y  -*■  0  defines  the  location  of  ehe  virtual  focua  F  . 

Inspection  of  Fig.  A1  shows  chat  5  ■  y.  In  the  limit  y  0,  it  is 
evident  that  the  distance  between  B'  and  8",  denoted  as  (B'B"),  becomes  a6 
where  a  is  the  sphere's  radius.  Inspection  of  Fig.  A1  also  shows  that 
(B"C")  *  (B'B")coa6R  as  Y  *  0  and  that  (B"C">a'Y*  Combining  these  results 
gives  a'  -  a  cos  0R.  By  definition,  a  -  a'  +  a(l-cos8R)  which  gives  a  -  a. 

The  result  a  -  a  could  have  been  anticipated  from  the  calculation  of 
a  in  Ref.  2  of  for  the  case  of  transmitted  waves  within  an  elastic  sphere. 

If  w*  examine  Eq.  (10)  of  Ref.  2  in  the  limit  that  n  (the  total  number  of 
chords  of  Che  transmitted  wave  within  the  sphere)  goes  to  °°  we  find  a  -*■  a. 


i." 


C'-JU'w  \."*>  ■*?' 
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3.1  Introduction 


Understanding  the  physical  nature  of  the  echo  structure  of  elastic 

objects  having  simple  shapes  is  an  important  step  toward  understanding  echoes 

from  complex  elastic  bodies.  Studies  on  acoustical  scattering  from  spheres 

and  cylinders1”11  data  back  to  Rayleigh.1  The  plane  wave  scattering  of  an 

elastic  sphere  or  cylinder  may  be  written  as  an  infinite  partial-wave  series 
2  3 

(PWS).  *  At  high  frequencies  the  PWS  converges  slowly  and  the  Sommerfeid- 

Watson  transform  (SWT)  may  be  used  to  convert  the  PWS  to  a  more  rapidly 

convergent  form.^  Furthermore ,  the  results  of  the  SWT  may  be  interpreted  in 

terms  of  reflected  and  transmitted  bulk  waves  and  the  scattering  contributions 

of  Franz,  Rayleigh,  and  whispering  gallery  waves.  In  this  way,  the  physical 

origins  of  tha  echo  structure,  are  elucidated.  Detailed  analysis  using  the  SWT 

has  bean  carried  out  on  fluid- loaded  elastic  cylinders.  *  The  "modified"  SWT 

has  been  used  to  study  the  case  of  scalar  plane  wave  scattering  from  a 
12  13 

transparent  sphere.  ’  "he  transparent  sphere  analyses  may  be  used  as  an 

aid  in  understanding  acoustical  scattering  from  fluid-loaded  fluid  spheres. 

In  the  present  paper  we  carry  out  the  SWT  on  a  fluid- loaded  elastic  sphere. 

The  analysis  for  elastic  spheres  differs  in  significant  details  from  the 

analysis  for  ocher  spheres 1^’1^  or  elastic  cylinders 4 

Our  SWT  analysis  concentrates  on  the  specular  reflection  and  Rayleigh 

cave  contributions  to  scattering  at  small  backscattering  angles.  Previously 

we  have  measured  and  modeled  the  angular  dependence  of  the  Rayleigh 

14 

contributions  to  near  backward  scattering.  The  SWT  confirms  the  physical 
picture  used  and,  for  the  first  time,  predicts  the  absolute  backscattering 
amplitude  associated  with  one  or  more  circumnavigations  of  the  Rayleigh  wave 
around  the  sphere.  Explicit  expressions  for  the  damping  of  Rayleigh  waves  on 
an  elastic  sphere  due  to  radiation  into  the  surrounding  fluid  are  obtained. 
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the  analysis  may  b«  extended  to  Include  whispering  gallery  waves.  Though  the 
emphasis  of  this  paper  is  on  the  backward  axially- focused  scattering  from 
spheres,  the  resulting  physical  picture  should  be  applicable  to  to  the  objects 
of  revolution  ensonified  along  the  symmetry  axis. 

To  test  the  SWT,  tungsten  carbide  spheres  In  water  were  ensonified  by 
tone  bursts  having  central  frequencies  such  that  24  <  ka  <  80  where  a  Is  Che 
sphere's  radius  and  k  is  the  wavenumber.  Measurements  were  made  of  the 
first  and  second  Rayleigh  contributions  to  Che  backscattered  pulse  train. 
Plots  of  the  measured  distinct  Rayleigh  amplitudes  as  a  function  of  ka  confirm 
Che  results  of  Che  SWT  and  illustrate  Che  significance  of  radiation  damping 
and  axial  focusing. 

In  Section  I  che  SWT  is  performed.  The  criteria  ka  >>  1  is  assumed 
throughout  Che  section.  The  physical  picture  which  ensues  from  the  analysis 
is  discussed.  The  analytical  prediction  of  Che  Rayleigh  contributions  to 
backscactering  follows  from  che  discussion.  In  Section  II  tungsten  carbide 
spheres  are  used  to  experimentally  confirm  this  prediction.  Throughout  the 
paper  we  relate  results  of  the  present  analysis  to  the  work  of  Ref.  14. 
Section  III  discusses  che  results. 

3.2  The  Sommerf eld-Watson  Transformation 

In  this  section  the  SWT  is  carried  out  on  the  total  pressure  field  of 
a  fluid- loaded  elastic  sphere  ensonified  by  a  plane  wave.  An  expression  for 
the  form  function  f  ensues  from  the  Initial  analysis.  Having  obtained  f, 
further  analysis  will  concentrate  on  its  interpretation.  The  main  emphasis  of 
the  analysis  will  be  on  the  specular  reflection  and  Rayleigh  contribution  to 
f  at  small  backscattering  angles.  However,  extensions  necessary  to  include 
other  waves  are  delineated.  Portions  of  the  present  analysis  parallel  the  SWT 
analysis  of  a  fluid-loaded  elastic  cylinder. 


Using  the  coordinacs  system  shown  in  Fig.  1  and  assuming  a  unit 


amplitude  plana  wavs  traveling  in  the  +z  direction  the  total  pressure  in  the 

11 


fluid  surrounding  the  sphere  is  the  real  part  of 


B 


p(r,9,t) 


ri£l)C  l  ifl(2n+l)(J  (kr)  +  ^  h^1}(kr))Pn(p)  (l) 

a-0  a 


In  this  expression  U  -  cos0,  k  -  w/c,  c  is  the  sound  speed  in  the  liquid,  r 


and  the  scattering  angle  8  are  defined  in  Fig  >1.1  and  h  are  spherical 

a  n 


Bessel  and  Hankel  functions,  and  Pq  is  the  Legendre  polynomial.  and 


are  3x3  determinants  whose  elements  are  given  in  Appendix  A.  They  are 
functions  of  x,  x_,  x_  ,  p  ,  p_  where  x  *  ka;  x  -  xc/c„;  c.  -  shear  sound 
speed  in  the  elastic  sphere;  x^  ■  xc/c^;  c^  *  longitudinal  sound  speed  in  the 
elastic  sphere;  pQ  -  density  of  the  liquid;  pg  -  density  of  the  solid. 

The  SWT  consists  of  rewriting  the  PWS  of  Eq.  (1)  in  terms  of  a  contour 


integral  using  the  relation 


15 


l  «(»  +  1/2) 

n“Q 


-IttA 


g(A) 


costtA 


dA 


(2) 


The  contour  T  is  shown  in  Fig.  2.  The  contour  surrounds  the  positive  poles 

of  (cosrrA)  .  Substituting  Eq.  (1)  into  Eq.  (2),  using 

C  ■  D  1  (kr)  +  B  h  (kr) ,  and  the  substitution  \>  •  A-  1/2  for  compactness 
n  q  ti  n  n 

of  notation,  we  have 


p(r,8) 


*<>) 


Vu) 


a-iTrA/2 

cosirA 


dA 


(3) 


We  proceed  by  deforming  the  contour  so  p  may  be  written  as  on  integral  over 
the  poles  of  1/DV  in  the  first  quadrant  (Section  II  discusses  a  pole  finding 
procedure).  The  general  location  of  the  poles  in  the  first  quadrant  is  shown 


Th«  r,  9,  $  coordinate  system  shown  above  wee  used  In  writing  the 
partial  wave  series  (PWS)  solution  to  plane  wave  scattering  from  am 

elastic  sphere.  The  plane  wave  is  assumed  to  be  tr  ■  iliag  In  Che  +2 
direction. 


in  7ig<  3.  The  location  and  identification  of  tha  poles  aa  Rayleigh, 
vhiaparing  gallary,  or  Franz  was  checked  using  tha  specific  case  of  an 
aluainua  sphere  with  the  same  physical  parameters  as  tha  aluminum  cylinder 


which  was  addressed  in  Raf  3.  Figure  3  also  shows  a  contour 

i  a  h 

7-  (  P  ,  Tq,  r*,  -  T,  r^)«  Thera  are  no  poles  within  the  contour. 

Therefore  tha  integral  of  Eq.  (3)  over  this  contour  would  vanish.  In  Appendix 
B  we  show  that  the  contributions  to  tha  integral  over  7  due  to  the  paths  at 
<*  are  negligible.  With  this  result  we  can  write  p  as  p  »  p^  +  where 


PI,II 


(4) 


and  the  respective  contours  are  ■  rQ  and  «  F*  .  We  limit  the 

analysis  which  follows  to  p^.  since  it  contains,  among  ocher  things,  the 

contributions  of  main  concern  to  us;  l.e.,  Che  specular  reflection  and 

Rayleigh  surface  waves,  p^  is  the  "background  integral”  which  has  been 

3  4 

argued  to  be  small  in  tha  cylindrical  case.  ’  We  will  return  briefly  to  p 
in  Section  III. 

To  evaluate  p^  we  first  break  F^  into  contours  F^  and  shown 

in  Fig.  4.  This  must  be  done  as  a  preliminary  to  separating  out  the 

3 

contributions  to  p^.  from  the  Incident,  and  bulk  waves.  Then,  In  the 

9  15 

integral,  we  use  Che  relations  ’ 


a11™  +  2i  ei7nJ  costtA  Q^(y)  ■  Py(|i) 


(5) 


Q^Cu)  -  >s[Pv(u)  +  (“)  Q^(p)  3  (6) 

where  and  Qv  are  Legendre  functions  of  the  first  and  second  kind.  The 

auxiliary  function  is  introduced  because  of  useful  properties^  when 


*  Franz  Potoa 


•  Rayleigh  Pole 
a  Whispering  Galery  Poles 


The  contour  T-  (T*  ,  To,  rf,  -r,  is  used  co  rawrica  cha 
concour  intagral  for  cha  praasura.  Tha  prassura  can  ba  written  as  a 
contour  intagral  ovar  TT  and  r  . 


*  Franz  Polaa 
■  RayMgh  Pol* 
o  Whtapartng  Gallary  Polaa 


Tha  coacour  To  of  Figura  3  can  ba  saparacad  inco  contours  f  and 
*r,  shown.  Tha  contour  r#  is  as  ad  to  find  eha  contribution  to  cha 
prasaura  from  cha  lncidanc  plana  wav a  and  cha  spacularly  raflactad 


wava. 


I X |  >>  1  (see  Appendix  C).  Substituting  Eq.  (5)  Into  Eq.  (4)  we  have 


PI  "  P1  +  p2  +  p3  wh*r* 

-13ir/4 


Px  "  21 


Fl(X)  Q^fti)  aiirX/2  dA 


(7) 


J1W/  iffA/2 

•  fi<»  *— r  « 


(8) 


-lirA/2 


-iir/4  . 

a  Fjtt)  Pv(y)  dA 


(9) 


and  F^A)  -  A  /Dy  .  Note  chac  the  integrand  of  no  longer  has  the  term 

-I 


(costtA)  and  therefore  it  has  no  poles  along  the  real  axis. 

2 


We  now  examine  the  expression  for  r  »  ka  (the  far  zone)  and  show 
chat  it  leads  to 

ikr 


lkz  ,  aa 
a  +  — - -  f 


2r 


(10) 


where  f  is  a  form  function  descriptive  of  the  scattering  for  ka  »  1. 


To  aid  in  manipulating  the  expressions  for  p^,  p, ,  and  we  rewrite 


F^(A)  using  the  relation 


Jv(kr)  -  H  (h<1}(kr)  +  hj2)(kr)) 


(11) 


and  the  result  from  Appendix  A  (Eq.  (A4))  that 


,(2) 


B  /  D- 

V  V 


ar/(x) 

-  1,(1  +  -^ -  >T) 


,d) 


(12) 


where  £)  ”,  and  .£)  *  are  given  by  Eq.  (A5)  and  H,  is  a  cylindrical  Hankel 


function.  Substituting  Eq.  (11)  and  (12)  into  *  D^J^(kr)  +  hj  (kr) 
and  using  this  rssult  in  ths  ralation  for  F^(X)  gives 


FX(X) 


|  <hf  «r> 


<2)<»)  $>:  ,  m 
'h<»<„  a; .» 


(kr)) 


Whan  this  expression  is  substituted  into  Eqs.  (7),  (8),  and  (9)  we  find  there 

(2) 

is  no  contribution  to  these  integrals  due  to  the  h^  (kr)  in  Eq.  (13)  since 

(2 )  (2) 

rt  and  r2  contain  no  poles  of  hr  (kr)  and  since  along  r_,  hj  (kr) 
vanishes  exponentially  as  a  function  of  X  over  the  portions  of  which  go 

to  ®.  (The  latter  assertion  may  be  established  using  the  methods  of  Appendix 
B.  See  also  Ref.  4  and  references  therein  for  the  cylindrical  case.) 

The  above  results  give  for  p^,  and  p^ 


-  -i 


e**i3ir/4  (kr)qU>  (y)  ^±irX/2  dA 


-  H 


-i3ff/4  r  i  >  iirX/2 

*  F2(X,x)h^'(kr)Pv(-u)5_dX 


P3  “  -  4 


■i  it /a  /i  \  *  iTrA/2 

*  F2(X,x)h^(kr)Pv(u)  ^-dX 


where  F2(X,x)  •  XH  X(2)(x)^  /hx(1,<x)£>;  .  Since  the  integrand  of  p^  contains 
no  poles  along  the  real  axis  we  may  deform  T  into  the  path  r  shown  in 
Fig.  4.  r  goes  through  the  two  saddle  points  found  to  exist  for  the 

integrand  of  p^.  The  right-hand  saddle  point  is  found  in  Appendix  C  to  lead 
to  the  incident  wave  exp(ikz).  This  and  the  approximation 
h,,  y(kr)ce  (kr)  exp[i(kr  -  vrr/2  -  ir/2)],  which  is  valid  for  kr  -1*  30 ,  allows  f  in 
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In  f,  the  contour  T  is  only  a  portion  of  the  contour  r  _  namely 
1  L  S 

that  portion  which  goas  ovar  tha  left-hand  saddle  point  since  the  right-hand 

saddle  point  has  already  yielded  the  incident  wave. 

Tha  problem  now  is  to  evaluate  f  and  interpret  the  results.  In 

examining  f  we  will  be  particularly  concerned  about  its  form  at  small  angles 

Y  relative  to  backscattering  where  it  -  9  ■  Y»  I Y I  <<  1  rad.  It  can  be 

verified  via  the  analysis  which  follows,  that  as  9  Increases  from  0  to  tt, 

the  left-hand  saddle  point  approaches  the  origin  in  the  complex  \  plane. 

This  implies  that  for  small  angles  relative  to  bsckscatterlng  the  contour 

encloses  no  poles  and  f^  ■  0.  Therefore  near  backscattering  we  need  only 

examine  f.  and  f_. 

1  2 

In  the  cylindrical  case  an  analysis  of  an  integral  similar  to  f^  has 

been  carried  out.*  In  Appendix  A  we  use  the  methods  of  Ref.  4  to  rewrite  ~ 

v 

a*  *  Vy  This  allows  one  to  write  F2(Atx)  -  Fg  +  F£w  with 

Fs  -  UsAhJ2)(x)/R1(1)(x)  and  -  -AUX  B^2  )  (  x) /Hx(1 )  (  x)  .  With  this 

separation  we  have  f^  -  fg  +  f£w  where  fg  and  f  are  given  by  Eq.  (17) 
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with  ?„(X,x)  replaced  by  F„  and  F  respectively.  The  motivation  for 
Z  s  tv 

this  separation  will  be  evident  from  the  results  which  follow. 

In  Appendix  D,  f  is  evaluated  via  the  saddle  point  method  giving 

fS  ■•*SA  ."12ka",('f/2)  (20) 

3 

where  R  ,  is  given  by  evaluating  Eq.  (A8)  at  the  saddle  point 
3 

X  ■  ka  siny/2.  In  the  limit  ka  -*■  »  one  can  show  that  (-R„,  )  becomes  the 
s  SAs 

reflection  coefficient  of  a  plane  wave  incident  onto  a  plane  solid** liquid 
boundary  at  angle  y/2  with  raapect  to  the  normal,  la  particular  for  y  ®  0 
one  finds 


<'Rs<V°>)  *  l0^ 


poe>/(Vi  *  °oc) 


R 


P 


(21) 


For  finite  ka  values  (-R^^)  i*  therefore  the  coefficient  for  specular 
reflection  from  the  sphere.  Equation  (20)  also  gives  the  propagation  phase  of 


the  specular  reflection  measured  relative  to  a  wave  traveling  in  liquid  along, 
the  path  (r  »  *,  0  ■  i)  •*  (r  ■  0)  -*•  (r  ■  «,  9  •ir-y).  We  have  previously 
obtained  results  analogous  to  Eq.  (20)  by  different  methods. *  (In  these 


previous  results  the  phase  of  the  specular  reflection  was  given  relative  to  a 
wave  traveling  in  the  liquid  along  the  path  (r»<»,  9»it)  -*•  (r»a,  0-tt)  -*■  ( 
9«rr-y)). 


Through  further  analysis  which  would  involve  a  combination  of  the 
methods  of  Ref.  4,  16-18,  one  should  be  able  to  recapture  from  f  the 

tW 

contribution  to  scattering  from  bulk  waves  transmitted  within  the  sphere. 
Since  f  is  not  of  direct  relevance  to  the  goals  of  the  present  paper  we 
will  not  examine  it  in  detail;  analysis  using  the  techniques  of  Refs.  16 


and  18  suggests  |ftw|  is  small  for  Cha  experiments  presented  in  Sec.  II  (a 

short  discussion  of  this  analysis  is  given  in  Sec.  III). 

A  caveat  concerning  the  evaluation  of  f^  should  be  noted.  The 

asymptotic  expressions  used  in  the  saddle  point  analysis  of  Appendix  D 

breakdown  near  6*0  and  9  ■  ir.  For  9  it  the  problem  lies  in  the 

asymptotic  expression  used  for  Q^^(u)*  Nuasensveig  has  shown  in  Refs.  12 

and  15  that,  for  scattering  from  a  fluid  or  soft  sphere,  the  saddle  point 

result  for  specular  reflection  remains  uniformly  valid  up  to  9  •  ir  even 

though  the  saddle  point  method  is  not  valid  near  9  ■  ir.  Using  Nussensvaig'  s 

methods  one  may  show  the  same  to  be  true  for  evaluation  of  f^.  For  9  -*■  9 

one  cannot  use  either  the  asymptotic  expressions  for  Q^(u)  or  the  Debye 

expressions  for  the  Hankel  functions.  In  the  neighborhood  of  9  "  0  one  must 

use  Airy  type  expansions  for  Che  Hankel  functions  and  new  effects  arise.15 

The  above  evaluation  of  f  Is  therefore  Incomplete  near  forward  scattering. 

3 

3  19  20 

Using  cha  mechod  of  residues  ’  ’  cha  integral  expression  for 

can  be  written  as  a  sum  over  Che  Rayleigh,  Franz,  and  whispering  gallery  poles 

»  • 

which  we  designate  as  A^,  A^,,  A^  respectively.  Hare  A  «  A^  is  a  root  of 
Che  equation  ,  *0.  From  the  residue  analysis  one  obtains 

f2  •  E  £i.  W) 


,  ,  2.i 

'  4  ta  00*- 


where  ,  is  the  derivative  of  §)  *  ,  with  respect  to  A  evaluated  at 

A  £-4  A 

A  ■  A^  and  the  relation  it  -  9  ■  y  v»s  used.  The  procedure  used  to  obtain 
Eq.  (22)  and  (23)  is  analogous  to  that  of  Ref.  7.  Further  manipulation  of 
f,  depends  on  the  type  of  pole  being  examined.  The  primary  difference  in 
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further  analysis  is  in  th*  asymptotic  expansions  used  for  the  Hsmal  functions 
in  Eq.  (23).  For  cha  Rayleigh  and  whispering  gallery  poles  the  Debye 
expansions 


,(1.2) 

'X 


(x) 


,1.1/2 


axp[-  i(x2-X2)1/>2  +  IX  arccos(X/x) 


lir/4] 


(* 


X2)i/4 


(24) 


may  be  used  since  the  poles  are  between  x  and  -x  in  the  X  plane  and 
|X„-x|  »  (X.)  (cf.  Fig.  3,  also  Appendix  A  of  Ref.  15  for  a  relevant 

jL  X* 

summary  of  asympeotic  expansions).  The  Franz  poles,  however,  are  in  a  region 

of  Che  X  plane  where  combination  of  asymptotic  expansions  must  be  used  (Ref. 

7,  Appendix  A).  Since  our  main  interest  is  in  the  Rayleigh-wave  form  function 

f  =  f  we  particularize  further  analyses  to  chat  case.  The  whispering 
XR  R 

gallery  contributions  to  have  the  same  form  and  interpretation  as  we  will 

find  for  fR. 

Q 

In  rewriting  fR  the  relation 

OD 

(coairXj)  1  -  2i  £  axp (i7r(2m+l)  (X^  -  1/2)]  (25) 

m**0 

is  useful.  Also,  for  the  near  backscactering  region  (|y|<<l  rad)  of  most 
interest  to  us,  we  can  use  the  following  approximation  for  the  Legendre 
function 12,13,21  valid  when  |  y| <<1  rad,  |  X|»if  |ReX|»|lmX| 


pxr*i  (coaY)  3  VV} 


(26) 


where  is  the  zeroth  order  Bessel  function  and  the  requirement 
| ReX I >> | Im X I  has  been  verified  numerically  for  the  Rayleigh  zeros  of 
aluminum,  tungsten  carbide,  and  fused  silica.  We  assume  |  Re  A  |  >>  1 1  m  X.)  for 
the  remainder  of  the  analysis.  Using  Eqs.  (24-26)  gives  Che  following 


approximation  for  the  Rayleigh  wave  form  function 


This  figure  show  eha  physical  picture  which  ensues  f tom  'ha  SWT 
analysis  of  eha  Raylaigh  contributions  to  scattaring.  Tha  incoming 
plana  wave  raprasancad  by  A B  allowo  eha  launching  of  a  Rayleigh 
surfaca  wava  at  3  which  clrcumnavlgatas  eha  sphars  whlla  raradlaclng 
back  Into  eha  surroundings.  At  point  3'  anargy  is  radiatad  in  eha 
backward  direction.?^  is  eha  virtu#!  point  sourca  from  which  ray  A'B 
and  eha  dashad  rays  eo  aithar  sida  appaar  eo  originate,  Whan  eha 

diagram  is  rotated  around  eUa  C'C  axis,  .ha  point  F  cracaa  out  a 

& 

virtual  ring- like  sourca. 
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fR  *  -<3r  «in  J0(Ary)  l  exp[iir(2rrH)(AR-i1)] 


where 


a  .  n  M  * 

R  R  ka 


n  -  2[-<(ka)2  -  X2)1/2  +  XR  cos-1  ^1  (29) 

The  Rayleigh  pole  satisfies  Che  condidon  0  <  Re  A  <  ka  and  ic  trill 

R 

be  convenient  co  define  Che  parameters:  SR  ■  Im(XR),  k^a  *  Re(XR)  and 
0R  ■  arcsin(k^k)  .  We  may  then  write 


and  our  assumption  jReAR|»|lmAR|  becomes  |kRa|»j6R|.  Using  Eq.  (30)  in 
the  exponentials  and  y  dependence  of  Sq.  (27)  and  ignoring  t»ms  of  0(B  R/  ) 
and  O(0R/ka)  we  find  that 


fR  =  "S  VkbR^  8 


-8r(2tT“28r)  inR  ®  io(2irl^a-iT)  ~27rm6R 

«  l  e  ‘  «  (31) 

m*0 


nR  a  -2  ka  cos0R  +  k^a  (2tt  -  20R)  -  tt/2  (32) 

where  we  define  b  ■  a  sin8  ,  This  expression  for  f  may  be  interpreted  in 
R  R  R 

terms  of  a  Rayleigh  weve  coupled  onto  the  sphere  at  a  local  angle  of  in>  Idence 

0 .  This  wave  repeatedly  circumnavigates  the  sphere  while  shedding  energy 

R 

bavk  into  the  liquid.  Figure  3  llluc  'races  the  interpretation  of  f  ,  The 

R 

ray  AA'  in  Fig.  3  may  be  attributed  to  the  m  •  0  term  of  f  .  AA'  includes  a 

R 

ray  traversing  Che  surface  of  the  sphere  from  B  to  B '  .  Evidently  k6  is  the 

R 
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propagation  constant  of  BB'  and  la  tha  phase  delay  of  AA'  relative  to 

a  hypothetical  ray  following  the  path  in  water  C"  -►  C '  -*■  C  -*•  C '  -*■  C  '  '  ; 

includes  a  ir/2  phase  advance  due  to  the  caustic  at  K.  The  exponential 

decay  exp{-8 _(2'n,-20_)]  is  due  to  continual  reradiation  of  energy  back  into 

the  liquid;  l.a.,  radiation  damping.  Reradiation  at  any  point  along  BB' 

occurs  at  angle  0R  with  respect  to  the  normal  to  the  sphere's  surface  at 

that  point.  To  understand  the  y  dependence  of  the  m”>0  term,  one  must 

examine  the  dashed  rays  to  either  side  of  A'B'.  These  dashed  lines  represent 

radiation  slightly  asav  from  backscaftaring.  Tracing  the  dashed  lines  and 

A'B*  backward  locates  tha  point.  7  from  which  they  appear  to  originate.  The 

spherical  symmetry  of  the  situation  th<  n  allows  the  rotation  of  the  figure 

around  the  C'C  axis  in  which  case  F_  traces  out  a.  ring  of  Mdlus  b  .  In 

R  F- 

Eq.  (31)  the  JQ(kbRY)  angular  dependance  followed  f-  on>  <■  ^  SWT;  however,  we 

previously  derived  this  radiafon  pattern  of  \l*tual  ring-like  sources  by 

other  methods  of  approximation. ^ ^ 

The  m  >  0  contributions  tc  f  art  due  to  the  continuation  of  the 

R 

ray  BB'  around  the  sphere  with  r. diction  in  the  backscattering  direction 
each  time  it  arrives  at  3'.  The  exp.-Z^mS  ]  term  in  Eq.  (31)  accounts  for 
radiation  damping  associated  with  m  circumnavigations  of  the  sphere.  T~\:. 
factor  exp(-imiT)  accounts  for  the  phase  advances  due  to  caustics  at  C  .a 


K.  The  on-axis  magnitude  of  the  contribution  to  f  by  the  afh  term  in  Eq. 

R 

(31)  is 


This  interpretation  of  the  SWT  is  consistent  with  the  analysis  and 
mcasuraments  of  the  axially- focused  Rayleigh  wave  contributions  to  the 
scattering  presented  in  Ref.  14.  The  absence  of  a  Y-dependent  phase  factor  in 
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Eq.  (31)  indicates  that  ehe  virtual  focal  clrcla  lias  in  a  plana  which 

contains  che  center  of  the  sphere  (see  discussion  below  Eq.  (37)  of  Ref.  17 

for  the  case  in  which  a  -a).  This  agrees  with  the  location  given  by 

n 

14 

geometric  methods.  The  on- axis  amplitudes  were  left  unspecified  in  Ref.  14. 

The  SWT  yields  Che  following  high-frequency  approximation  to  the  total 
form  function  away  from  the  forward  direction 


f  =  f. 


+  4  f 

A-WG.F  X, 


where  Che  summation  is  over  all  whispering  gallery  and  Franz  poles, 
will  be  of  che  same  form  as  Eq.  (31)  except  chat  parameters 


The 


are  to  be 


calculated  from  che  X^,  . 


3,3  Experimental  confirmation  of  Rayleigh  backacattering  amplitudes 

In  this  section  tungsten  carbide  spheres  are  used  to  test  che  Rayleigh 

backacattering  amplitudes  of  Eq.  (33)  fo c  the  cases  m  ■  0,1.  We  first 

outline  che  methods  used  to  obtain  numerical  results  from  Eq.  (33).  We  then 

briefly  summarize  the  experimental  procedure.  Last,  we  compare  analytical  and 

experimental  results.  The  reader  Is  referred  to  Refs.  14,  16,  and  18  for  a 

description  of  che  experimental  apparatus. 

The  backacattering  amplitudes  An  of  Eq.  (33)  are  functions  of  ka. 

As  a  first  step  in  calculating  Ag(ka)  and  A^(ka)  for  tungsten  carbide  we 

determined  the  Rayleigh  pole's  ka  dependence.  A  computer  program  wa? 

3  20 

developed  based  on  the  "winding  number  formula"  '  to  find  the  complex 
zeros  of  in  Eq.  (12).  The  program  was  used  to  fino  XR  for  a  tungsten 

carbide  sphere  for  23  values  of  ka  between  19  <  ka  <  83.  The  ka  values  usui 
and  values  obtained  are  given  in  Table  I.  Also  given  In  Table  I  are  the 

material  parameters  used  for  che  tungsten  carbide  sphere.  The  motivation  for 

14 

Che  particular  ka  values  chosen  in  Table  I  comes  from  our  previous  work 


have  uncertainties  of  1-  (0.0005  -1  10.0005) 
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using  Resonance  Scat  caring  Theory  (RST).  In  RST  one  deals  with  Che  partial 

£ 

wavs  series  of  Eq.  (1)  directly.  In  Ref.  14  ve  found  the  complex  ka  zeros  of 

D  in  Eq.  (1)  for  integer  vslues  of  a  from  n  ■  0  to  n  *  30.  The  ka 
n 

values  used  in  Table  1  are  the  real  parts  of  Che  complex  ka  zeros  found  in 

Ref.  14.  It  is  not  necessary  chat  the  ka  values  be  chosen  In  this  manner; 

however,  we  have  done  so  to  facilitate  discussion  in  a  subsequent  paper.  The 

quantities  Re  X^  ■  k^a  and  1^^  X^  •  0^  were  found  as  continuous  functions 

of  ka  by  curve  fitting  the  discrete  values  of  Table  1. 

The  fitted  curves  for  k_a  and  g  facilitated  the  numerical 

K  R 

evaluation  of  Aq  and  in  Che  range  20  <  ka  <  80.  The  evaluation  of  GR 

in  Eq.  (28)  requires  the  evaluation  of  Bessel  functions  of  complex  order  (see 
Eq.  (A5)),  The  procedure  used  in  this  evaluation  was  based  in  part  on  the 
analysis  given  in  Ref.  22. 

Figure  6  gives  a  conceptual  diagram  of  the  apparatus  used  to  test  the 
SWT  predictions.  A  piston- like  transducer  produces  a  sinusoidal  tone  burst  4 
cycles  in  duration.  This  tone  burst  is  scattered  from  either  a  1.27  cm  or 
2.54  cm  diameter  tungsten  carbide  sphere.  In  this  experiment  Che  needle 
hydrophone  which  picks  up  the  resulting  echo  was  sac  at  y  ■  0;  i.e.,  along 
Che  backscattering  axis.  Figure  7  shows  the  structure  of  the  backscattared 
echo  for  the  case  ka  ■»  49.1.  The  figure  displays  the  amplified  voltage  from 
the  hydrophone  averaged  over  repeated  echoes.  After  the  specular  echo,  the 
principal  contributions  are  the  Rayleigh  echoes.14  The  frequency  of  the 
incident  burst  could  be  varied  from  8CC  kHz  up  to  1.5  mHz.  The  corresponding 
range  of  accessible  ka  is  from  25  to  80. 

The  experimental  procedure  was  to  measure  the  peak-eo-poak  voltage  of 
Che  central  cycle  for  each  echo  structure  of  interest.  Those  voltages  will  be 
designated  as  V  for  the  specular  echo  and  as  V  for  the  first  (m  •  0)  and 
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Simplified  diagram  of  cha  scattering  experiment.  The  backacattaring 

angla  Y  la  zaro  ia  cha  preaeac  experiment.  Tha  diagram  Is  aoc 

drawn  co  acala .  In  all  cha  experiments  daaerlbad  in  Sacclon  II , 

z  3  160  at,  while  z  3  33  cm  for  axparimanca  with  cha  2.54  as 
P 

dlametar  cungacan  carbida  aphara  and  z  a  19  cm  for  axparimanca  with 
1.27  cm  aphara. 


Fig.  6 


in  oscilloscope  erne*  of  :ha  backacaetariag  echo  from  a  cungacan 
carbide  sphara  ansoolflad  by  a  cone  burse.  The  Sea  of  cha  spbara  la 
approximately  49.1.  The  individual  achoaa  ara  labeled  as  follows: 

-  specular  reflection,  3  -  Rayleigh  surfaca  wave  achoaa,  C  -  achoaa 
whose  specific  origin  was  hoc  identified  but  which  ara  conjectured 
ba  other  surfaca  waves  such  aa  "Whispering  Gallery"  waves. 
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second  (m  -  1)  Rayleigh  echoes.  The  voltage  measurements  were  made  on  the 
averaged  echo  from  256  bursts.  An  absolute  calibration  of  the  apparatus  would 
require  that  signal  attenuation  associated  with  propagation  through  the  water 
be  accounted  for.  Instead  of  absolute  scattering  amplitudes,  we  will  concern 
ourselves  with  the  Rayleigh  echo  amplitudes  relative  to  that  of  Che  specular 
reflection  for  the  particular  value  of  ha  under  consideration;  these  will  be 
designated  as  RQE  -  VQ/VS  and  R^  -  V1/Vg. 

In  order  to  compare  RqE  and  R^  with  the  SWT  prediction  for  the 
and  A^ ,  it  was  necessary  to  also  evaluate  the  SVT  predictions  for  Che 
specular  amplitude  at  Y  *  0  as  a  function  of  ha.  The  required  contribution  to 
the  form  function  is  I  with  -  0;  see  Eqs.  (20)  and  (A3).  We 

designate  this  quantity  as  .  Evaluation  of  Eq.  (A8)  shows  chat  as  ha 

increases  from  18  to  83,  R^  increases  monoconically  and  approaches  the 
specular  reflection  amplitude  for  a  plana  tungsten-carbide  water  interface, 

|R  |  *  0.9693  from  Eq.  (21).  At  ha  -  20,  R  was  only  0.1X  below  |R  | . 

P  5T  p 

Since  RS>J,  varies  by  only  about  0.1Z  in  the  ha  region  of  interest,  the 

theoretical  echo  amplitude  ratios  R__  -  A„/R„_  and  R,_  »  A, /R are 

0T  0  ST  IT  1  ST 

dominated  by  the  ka  dependences  of  AQ  and  .  These  ratios  ara  plotted  as 
the  solid  line  in  Fig.  3  together  with  the  experimental  amplitude  ratios  R 
and  R^  .  The  general  agreement  appears  to  confirm  the  predictions  of 
Eq.  (33). 

A  couple  of  comments  on  the  relation  between  the  experimental  and 
theoretical  results  are  in  order.  First,  the  theoretical  prediction  uses 
material  parameters  from  Ref.  11  since  the  elascic  parameters  of  the  actual 
spheres  used  could  not  be  determined.  The  material  parameters,  of  the  spheres 
used  in  the  experiment  may  be  slightly  different.  For  example,  the  densities 
of  the  spheres  used  were  measured  to  be  14.75  gm/cm  and  14.95  gm/cm  for  the 


Experimental  and  theoretical  results  for  the  backscactering 
amplitudes  of  cha  first  and  saeond  Rayleigh  contributions  as  a 
function  of  lea.  The  amplitudes  ara  normalized  Co  cha  specular  echo 
amplitude  at  cha  indicated  value  of  lea.  Tha  does  ara  experimental 
results  for  eha  2.34  cm  diameter  sphere,  Che  triangles  ara  results 
for  cha  1.27  cm  diameter  sphere,  and  the  solid  lines  theoretical 
resales  using  Eq.  (33). 


I 
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large  and  small  spheres  respectively.  This  may  account  for  some  of  Che 

difference  becveen  experiment  and  theory.  Second,  the  SWT  analysis  Is  a 

steady  state  method.  In  the  experiment,  however,  pulses  must  be  used  so  chat 

the  Rayleigh  contributions  can  be  isolated  in  time.  The  validity  of  applying 

a  steady-state  analysis  to  experiments  of  the  type  described  above  has  been 
16  23 

examined  *  in  connection  with  waves  transmitted  within  the  sphere.  It  was 
found  that  though  focused  Cone  bursts  could  be  distorted,  the  peak-co-peak 
amplitude  of  the  central  cycle  of  a  four-cycle  burst  should  be  well 
approximated  by  that  of  a  steady  state  analysis.  It  Is  anticipated  that  a 
detailed  analysis  of  the  shapes  of  transient  Rayleigh  echoes  would  lead  to 
similar  conclusions. 

It  should  also  be  noted  that  the  receiver  hydrophones  were  not  always 

in  the  far-zone  of  Che  scattering.  For  example,  in  the  far-zone  condition 
2  2 

r  »  ka  ,  a  «  1.27  cm  and  ka  ■  50  gives  ka  *63.5  cm.  Fortunately,  however, 

this  is  not  a  serious  drawback.  This  is  because  the  ratios  of  Rayleigh  and 

specular  echo  amplitudes  should  be  well  approximated  by  |a  |/R._  for 

m  ST 

2 

distances  much  closer  than  the  condition  r  >>  ka  .  The  radiation  from  a 

virtual  ring- like  source  at  radius  b_  is  modeled  in  Section  V  of  Ref.  16  for 

R 

1/3 

Che  region  r  >>  bR  (kbR/2  )  =  r^,  the  Fresnel  region.  This  analysis  may 

be  applied  to  Rayleigh  echoes  in  the  present  problem  since  these  appear  to 
originate  from  virtual  ring- like  sources;  ka  ■  50  and  a  -  1.27  cm  correspond 
to  kbR  3  18  and  r^  3 0.7  cm.  The  measurements  were  taken  with 
r  -  a  +  z  3  33  cm.  The  near- field  corrections  required  to  predict  r.h--t  f air¬ 
field  amplitude  ratios  from  Che  conditions  of  this  measurement  wiuru  eat  La-actd 


to  be  <  IX. 
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3.4  Difcusaion 

As  was  noted  in  Saetion  II,  tha  theoretical  amplitude  ratios  RQT  and 

R^  plotted  in  Fig*  8  have  dependences  on  ka  which  are  dominated  by  the 

Rayleigh  echo  amplitudes  AQ  and  A^ ;  the  R^  may  be  approximated  to  within 

0.12  by  A  /0.9693.  The  qualitative  features  of  these  plots  may  be  understood 
in 

14 

as  follows.  We  have  previously  demonstrated  that  tha  Rayleigh  echoes  are  an 

example  of  axially  focused  scattering.  Consequently,  it  is  to  be  anticipated 

14 

that  over  a  range  of  ka  (consistent  with  the  phase  matching  condition  ),  the 
amplitude  coupling  factor,  GR  in  Eq.  (38),  increases  with  ka.  Thus  AQ  and 
A^  are  expected  to  increase  with  ka  until  the  exponential  factors  in  Eq.  (33) 
are  dominant.  For  tungsten  carbide  the  damping  parameter  also  increases 

with  ka.  Tha  dominance  of  radiation  damping  is  evident  in  the  plot  of  R^ 
for  ka  40.  Tha  effects  of  damping  are  more  significant  in  A^  then  in  A^ 
since  A1  -  AQ  exp(-2Tr8R). 

The  qualitative  features  mentioned  above  are  clearly  evident  in  the 
data.  Some  discussion  of  tha  discrepancies  between  observation  and  theory  are 
in  order.  As  noted  in  Sec.  II,  one  plausible  cause  of  error  was  the  use  of 
tabulated  elastic  parameters  for  tungsten  carbide  instead  of  the  unknown 
parameters  for  the  actual  sphere  studied.  (There  is  other  evidence  that  the 
parameters  may  depend  on  the  manufacturing  process  of  the  tungsten  carbide.) 
The  differences  between  RQR  for  the  large  and  small  spheres  in  the  region  of 
overlapping  data  (ka  close  to  43)  is  suggestive  of  tha  magnitude  of  systematic 
experimental  uncertainties. 

A  further  source  of  discrepancy  may  be  inferred  from  the  structure  of 
the  data.  Tha  deviation  from  the  thsoretlcal  curve  is  noc  random  but  seems  to 
have  a  structure  (especially  RQ2  for  50  £  ka  £  80).  This  is  suggestive  of 
possible  interference  of  Che  Rayleigh  wave  with  ocher  smell  amplitude  waves. 
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We  have  used  the  methods  of  Ref.  16-18  to  find  the  amplitudes  and  times  of  the 
transmitted  wave  contributions.  The  most  significant  transmitted  wave 
contributions  have  amplitudes  which  are  between  0.1  and  Q.15  times  that  of  the 
specular  reflection  and  are  delayed  by  about  17.5  and  32.7  usee  relative  to 
the  specular  reflection.  The  timing  suggests  these  waves  may  be  in  part 
responsible  for  the  two  left  most  pulses  labeled  C  in  Fig.  7.  The  largest 
amplitude  transmitted  wave  in  the  time  windows  of  the  first  two  Rayleigh 
contributions  had  an  amplitude  of  about  .07  times  the  specular  reflection  and 
could  be  a  source  of  the  structure  of  Fig.  8. 

A  related  source  of  error  may  be  the  omission  of  whispering  gallery 

waves.  As  noted  earlier,  the  analysis  of  Sec.  II  for  Rayleigh  waves  is  also 

applicable  to  whispering  gallery  waves.  We  have  not  carried  out  the  somewhat 

tedious  computation  of  the  whispering  gallery  amplitudes  for  tungsten  carbide 

since  an  alternate  analysis  suggests  their  amplitudes  are  small  in  the  ka 

region  of  the  experiment.  This  alternate  analysis  involves  approximating  the 

exact  total  form  function  f  for  a  tungstan  carbide  sphere  by  f  •  f  +  f  • 

™ ™ ™ l>^™,  £  ASK 

see  Eq.  (34).  We  have  compared  jf  |  with  |f  |  computed  using  the  PWS^ 

A  1 

and  find  that  If^l  correctly  reproduces  | f ^ J  except  for  a  fine  structure 
(superimposed  on  If^l)  due  co  whispering  gallery  waves  or  other  waves. 
Besides  suggesting  the  whispering  gallery  contributions  are  small,  this 
comparison  is  relevant  to  the  interpretation  of  resonance  scattering  theory 
and  a  detailed  discussion  is  beyond  the  scope  of  the  present  paper. 

Another  poasible  source  of  the  experimental  and  theoretical 
differences  as  well  as  diffsrencss  between  the  overlapping  data  for  the  two 
.spheres  is  the  sphere  mounting  system.  Considerable  time  and  effort  was  spent 
finding  a  mounting  system  which  had  a  minimum  effect  on  the  scattering 
experiment.  Even  with  this  one  must  assume  that  the  effect  of  the  mounting 
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varied  between  the  two  spheres  and  had  an  impact  on  the  attenuation  if  the 
Rayleigh  waves  at  least  locally  in  the  region  of  contact  between  sphere  and 
mount . 

Finally,  it  is  important  to  remember  that  below  Eq.  (4)  we  separated 

the  pressure  into  contributions  P  and  P  .  The  analysis  after  Eq.  (4) 

3  4 

concentrated  on  P^.  is  known  as  the  background  pressure  ’  and  has  been 

4 

shown  to  be  negligible  in  the  fluid- loaded  cylinder  case.  The  comparison  of 
jf^j  and  |f^.|  discussed  above  implies  P  is  small  for  tungsten  carbide 
spheres.  However,  one  cannot  discount  a  priori  the  possibility  of  P 
contributions  being  responsible  for  part  of  the  difference  between  experiment 
and  theory.  P^  is  also  important  when  interpreting  f  of  Eq.  (10)!  f 
must  be  interpreted  as  a  high  frequency  approximation  to  fT  which  excludes 
the  contribution  from  P^. 

We  have,  within  the  limitations  discussed  in  Sac.  I,  performed  a  SOT 

on  a  fluid-loaded  elastic  sphere.  We  paid  particular  attention  to  the 

specular  reflection  and  Rayleigh  contributions.  Further,  subject  co  the 

discussions  in  Sec.  II  and  this  section,  we  have  experimentally  tasted  the  SOT 

prediction,  Eq.  (33),  for  the  amplitude  of  the  first  two  Rayleigh 

contributions  to  backscattering .  The  general  agreement  between  experiment  and 

theory  appears  to  confirm  this  particular  SOT  result.  The  Bessel  function 

angular  dependence  given  by  the  SOT  (cf.  Eq.  (31))  was  experimentally  confirmed 
14 

in  previous  work.  It  may  be  shown  that  kb  y  is  equivalent  to  the  argument 

tv 

of  JQ  which  was  predicted  and  confirmed  in  Ref.  14. 
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Appendix  A 

In  this  appendix  wa  give  Che  elements  of  the  3x3  determinants  Bq 
and  Dq  of  Eq.  (1).  We  also  manipulate  the  ratio  (where  v  ■  X  -  1/2) 

to  obtain  the  results  used  in  Eqs.  (12)  and  (20) - 
From  Raf.  11  we  have  chat 


12 


‘22 


-  5  \  Vx) 

d12 

^a 

d13 

*  j' (!) 

n 

d22 

d23 

1 

(Al) 

0 

d32 

d33 

p  X2  h^ (x) 

4 

4 

-  X  h(1)*  (X) 

n 

d22 

4 

> 

(A 2) 

0 

d32 

d33 

[ 2n (n+1 )  -  x2]j  ( 
s  n 

V  “  4  *L 

w  * 

2n(n+l)[xs  j^(xs) 

-  w] 

> 

XL  :!n('xI.)  ’  d23 

-  n(n+l)j 

n(xS}  * 

(A3) 

d32  "  2Cja(xL}  ~  \  ^n^L^ 


‘32 

in 

33 


2  xs  j;(xs)  +  [x|  -  2n(n+l)  +  2]  j^)  . 


(1) 


In  these  expressions  J  and  h  are  spherical  Bessel  and  Hankel  functions 

n  a 


respectively.  The  primes  denote  differentiation  with  respect  to  the  argument. 
The  definitions  of  all  other  variables  are  as  listed  below  Eq.  (1)  with  the 
additional  definition  that  p  m  P  /p_. 


The  expression  By/Dy  is  given  by  the  ratio  of  Eq.  (Al)  and  (A 2)  with 
n  replaced  by  v.  By  using  the  relations^4  j^(x)  *  l/2(h^(x)  +  h^  (x)  ] 
and  h^’^(x)  ■  (tt/2x)^^  H^'^(x)  one  can  rewrite  the  ratio  B  /DV1  as 


B  /D 

v'  v 


hx2)(x)S)^ 

^  (1  +  H,(1)(x)^) 


"  Zl,2  +  0<1 


where 

q  *  xs(d22d33  “  d32d23)/(d12d33  "  d32d13) 

(A6) 

zt  -  x  h^i)'(x)/h^i)(x) 

Equation  (A3)  Is  a  convenient  fora  for  residue  analysis  leading  to  Eq.  (22)  and 

(23).  However,  for  the  saddle  point  analysis  of  Section  must  be 

rewritten.  The  object  of  developing  an  alternate  expression  for  £)v/£*  is 

to  separate  out  a  taro  which  reduces  (in  the  limit  ka  -  ®)  to  the  reflection 

coefficient  of  a  plane  solid,  liquid  interface.  A  related  analysis  was  carried 
4  . 

out  by  Brill  tor  the  cylindrical  case.  Following  Section  II  of  Ref.  4  one 
can  rewrite  (after  considerable  algebra)  as 


-  RsX  -  UX 


RsA  "  (z2  B22  +  ^A22^/(-zl  B22  +  ^A22^ 


g(z2  ~  zl} 


A22Y2  •  B22Y1 


(*tB„  +  PA.,)  1  +  (z,I.  +  pr,)/(*,B„  +  pA*_) 


Aij  ■  ziL[(zjs +  l^2x?* 1/43  * 

Bij  -  [v(vH.)/x*  -  m  xj)-  1/4]  +  2v(v+l)/x4  -  2j^(2x^)  + 

21L[v(xh-1)/x4  -  2/x4  -  1/xj]  +  ZjgZ^IMv+D/Xg  -  2/x^ ] 

h  “  *12*1  +  A217S  +  *UVS 

Y2  “  B127L  +  B21yS  +  B117L7S 

7A  *  hi1)(xA)/hv2)(xA)’  ZU  "  XA  }  <*A)/hi4>(xA) 


where  A  ■  L  or  S.  The  significance  of  Eq.  (AS)  is  discussed  after  Eq.  (20) 
and  in  Appendix  D. 
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In  chia  appendix  we  wane  to  ahow.  chat  the  integrals 


p(l,2,3) 


JrU,2,3) 


...  C,  -IitX/2 

•  X(^)  »„(**>  *5^*-  dX 


(B1 


are 


negligible.  are  the  three  contours  at  infinity  shown  in  Fig.  3 


.(1) 


and  Cy  *  Dv  jy(kr)  +  By  h^  ' (kr) .  We  first  examine  B^/D^  in  the  limit 


|  X  |  ->  «.  To  examine  By/Q^  we  use  Che  relations 


>1/2 


jv(x)  -  (tt/2x)  Jx(x)  ,  (v  -  X  -  1/2) 

(i) 


(B2; 


sr (x) 


(7T/2x)1/2  U^Cx) 


(B3: 


between  spherical  and  cylindrical  functions  and  the  result  from  Eq.  (A6)  of 
Ref.  15  that 

-I/O  l 

(B4] 


Jx(x)  -  (2ttX)*  1/2(ex/2X)A  as  |X|  -« 


where  e  ■  (2.718 


j^(x) 


Vx) 


v 

—  as 
x 


).  From  Eq.  (B2)  and  (B4)  it  follows  that 

I X  I  ■*  « 


(B5) 


where  the  prime  denotes  differentiation  with  respect  to  x.  This  result  can 


be  used  in  the 


j 


of  Eq.  (A6)  to  find  their  limiting  values  as 


X 


V  V  V  V  V 

Using  these  limiting  values  we  find  that  ^.jd^  -  <i32ci33'>  0  aa 


jXj  ■*  «.  Finally,  substituting  this  last  result  into  Eqs.  (A4)  and  (A5)  we 
find 

3  . 


X  - 


(B4; 


for  all  regions  in  the  complex  X  plane. 


Using  Eq.  (B2),  (B3) ,  (B4)  and  Fig.  15  of  Ref.  15  one  can  also  show 
thee  j^(x)/h^(x)  ■*  J^(x)/H^  (x)  as  |X|  ■*  ®.  This  and  Eq.  (B2)  and 
(B3)  allow  us  Co  write 


,(1,2,3) 


-iw/4  ,  ir  J,LX<*kr) 


r (1,2,3) 


(21cr}1  (1)  '  MX(w)dA 

Hx  (x) 


where 


Lx(x,kr)  -  H<lV(x)  Jx(kr)  -  J' (x)Hj[1)  (kr) 


MX(U) 


XP  (u)  «xp(rl^/2) 
v '  '  cosirX 


(B5 


Now,  using  Appendix  A  of  Ref.  15  along  with  Eq.  (2.21)  and  (2.24)  of  che  same 
reference  and  che  substitution  X  -  R  exp (10)  one  can  establish  on  che  con¬ 
tours  (where  R  ■*•  *)  that:  and  behave  no  worse  than 

(1) 1 

exp(R),  and  11^  (x)  behaves  as  axp(RinR).  Thus  the  integrand  in  Eq.  (B5) 

vanishes  on  and  therefore  *  0.  In  establishing  this 

behavior  it  is  important  to  note  that  che  curve  h^  in  Fig.  3  is  che  curve 
along  which  the  zeros  of  (x)  <.  re  located.^ 


In  this  appendix  we  examine  the  integral 


inc 


-i 


-i3ir/4  _  . 

e  F2<1,x) 


h^1)(kr)Q^1)(cos0)elirX/2 


dA 


(Cl) 


M 


where  the  contour  Is  the  portion  of  in  Fig.  4  which  goes  through  the 

right  hand  saddle  point.  In  identifying  r_,  the  definition  of  I*.  below 

ri  L 

Eq.  (19)  was  used.  We  show,  using  the  saddle  point  method,  chat  Eq.  (Cl) 
represents  the  contribution  to  p  (cf.  Eq.  14)  from  the  Incident  plane  wave. 

To  evaluate  Eq.  (Cl)  we  must  first  find  the  saddle  point  location.  To 
do  so  we  rewrite  the  integrand  of  Eq.  (Cl)  in  a  form  valid  whan  x  <  ReA  <  kr 
(specifically  part  of  region  A  in  Fig.  15  of  Ref.  15).  In  this  region  we  can 
use  the  relation  (Appendix  A  of  Ref.  15) 

^2)00  si  -»x1)(x)  (C2) 

This  relation,  Eq.  (A5)  and  (A6) ,  and  the  definition  F^d.x)  - 
(x)5^  allows  one  to  show  that  -  1  and 

F2(A,x)  »  -  A  (C3) 

Also,  when  |  A  |  »  1,  e  <_  0  <_  tt  -  e,  )A|e  »  1  we  can  use  the  relation15 

Q«)(„.e,  ,  zJlBXL  .  {ct) 

(2irAsin0) 

We  return  briefly  to  the  requirement  j A | e  »  1  at  the  end  of  this  appendix. 
The  requirement  j  A |  »  1  is  fulfilled  in  the  region  RaA  >  x  since  at  the 
outset  of  the  paper  we  assumed  x  >>  1.  Finally,  since  we  are  examining  the 
region  where  ReA  <  kr,  we  can  rewrite  h^(kr)  using  Eq.  (B3)  and  (24). 


With  chase  results,  and  Che  change  of  variables  \  •  kr  sinw,  p 


becomes 


5(w)  -  e  [kr  sinw  cosw/2tt  sin8]^^ 

where  s(w)  ■  8) sinw  +  cosw  -  (Jjrr  -  w)sinw]  and  is  the  trans¬ 

formed  contour  of  F  which  passes  through  the  saddle  point  w  *  w  .  One 

n  3 

20  25 

can  evaluate  Eq.  (C5)  using  the  saddle  point  method.  ’  We  find  Che  saddle 

point  at  w^  *  8  (X  ■  Icr  sin8)  which  gives  a  value  for  Pinc 


[kr  s(w3)]  , 

* 


-27 r  .1/2 

i"(w  )J 


S(w  ) 

9 


2  2 

where  s''(w  )  is  d  s/dw  evaluated  at  w  ■  w  .  The  right  side  reduces  to 
S  3 

exp(lkz)  which  corresponds  to  the  Incident  plane  wave  as  required. 

The  condition  Chat  Eq.  (C4)  is  applicable  at  the  saddle  point  is 
|X  |e  -  kr  sin@  »  1.  Consequently,  for  points  close  to  the  backscactering 
axis,  9  close  to  it,  this  proof  is  only  strictly  valid  for  large  values  of 


Ia  this  appendix  we  use  Che  saddle  point  method  Co  approximate  Che 


contain  no  exponentials  and  is  a  slowly  varying  function  of  X.  Since 

20  25 

ita  »  1  the  Integral  can  be  evaluated  via  the  saddle  point  method.  ’  In 

terms  of  the  baclcscattering  angle  y  ■  tt  -  9,  the  saddle  point  is  at  '  w  y/2 

s 

foe  which  X  •  lea  sinw^  =  Xg.  The  approximation  of  the  integral  in  Eq.  (D3) 
is  given  by  the  right  side  of  Eq.  (C6)  with  s,  r,  and  S  replaced  by  t,  a, 
and  T,  respectively;  the  resulting  approximation  for  fg  is  given  by  Eq.  (20) 
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Chapter  4 

RESONANCE  SCATTERING  THEORY  REVISITED  VTA  THE  SOMffiRFELD-WATSON 


TRANSFORMATION  FOR  SCATTERING  FROM  SPHERES 


4.1  Introduction 


Resonance  Scattering  Theory  (RST)  originated  from  the  application  of 

1  2 

nuclear  scattering  theory  to  classical  scattering  problems.  ’  It  has  been 

used  in  addressing  a  number  of  different  scattering  situations.  (Reference  2 

discusses  several  specific  problems  which  have  been  examined  via  RST) . 

Because  of  these  applications  of  RST  an  examination  of  the  limitations  of  some 

of  the  basic  concepts  of  RST  is  justified.  We  stress  at  the  outset  that  these 

limitations  do  not  affect  the  validity  of  most  results  obtained  with  RST  but 

they  do  affect  interpretation  of  these  results. 

RST  achieves  a  separation  between  rapidly  varying  "resonance”  portions 

of  scattered  signals  and  a  slowly  varying  background.  In  particular,  for  an 

acoustical  plane  wave  scattered  from  fluid- loaded  elastic  bodies  such  as 
3  4 

spheres  or  cylinders,  RST  has  been  used  to  Interpret  the  rapid  fluctuations 
in  the  backscattered  pressure  as  the  ka  of  the  bodies  are  varied  (where  k  is 
the  wave  number  of  the  acoustical  plane  wave  in  the  fluid  and  a  is  the 
radius  of  the  sphere) .  One  of  the  important  numerical  tasks  in  RST  is  to  find 
the  complex  "resonance"  ka  values  which  locate  poles  of  the  scattering 
amplitude.  Two  basic  RST  ideas  are  then  used  in  interpreting  scattered 

3 

pressures.  For  a  sphere  these  ideas  translate  into  the  following  statements  : 
(a)  each  of  the  many  resonance  ka's  can  be  labeled  with  the  integers  n  and 
l  and  at  the  (n,2)th  resonance  n  +  (1/2)  wavelengths  of  the  £Ch  surface 
wave  fit  onto  the  circumference  of  the  sphere,  (b)  that  any  rapid  fluctuation 
in  the  scattered  pressure  is  due  to  the  in-phase  adding  of  a  specific  surface 
wave  and  these  fluctuations  are  thus  associated  in  a  simple  fashion  with  the 
set  of  modal  resonances  of  the  target.  Two  of  the  major  goals  of  this  paper 
are  to  examine  the  limitations  of  these  ideas. 


Our  examination  uses  the  particular  problem  of  an  acouatical  plane 
wave  scattered  from  a  fluid- loaded  elastic  sphere.  We  use  a  recently 
completed  Sommerfeld-Watson  transform  (SWT)  analysis  as  an  aid  in 
investigating  RST.  The  SWT  allows  separation  of  the  contributions  to 
scattering  from  reflected  waves,  transmitted  bulk  waves,  and  surf Ice  waves. 

Section  I  contains  an  examination  of  the  relation  between  the 
resonance  ka's  and  the  rapid  fluctuations  in  the  backscattering  form  function 
of  an  elastic  sphere.  To  perform  the  examination  we  use  the  case  of 
backscattering  from  a  tungsten  carbide  3phere  with  10  ka  80.  We 
approximate  the  exact  partial  wave  series  (PWS)  solution  for  the 
backscattering  form  function  of  the  tungsten  carbide  sphere  by  using  a  portion 
of  the  SWT  results.  This  process  allows  a  check  of  the  SWT.  Once  checked  the 
SWT  results  are  useful  in  understanding  the  relation  between  resonance  ka's 
and  the  backscattering  form  function  fluctuations. 

Throughout  the  paper  we  define  physical  resonances  as  the  real  ka 
values  for  which  an  Integer  plus  one-half  wavelengths  of  any  surface  wave  fit 
on  the  sphere.  In  Section  I  we  assume  that  the  real  part  of  the  RST 
resonances  give  the  physical  resonance  values  of  the  sphere.  In  Section  II  we 
show  that  these  RST  results  actually  only  approximate  the  physical  resonance 
values.  The  error  inherent  In  this  RST  approximation  is  small  enough  for 
tungsten  carbide  that  It  has  no  effect  on  the  discussion  of  Section  I.  We 
conclude  in  Section  II  that,  when  one  approximates  the  physical  resonance 
values  by  using  RST  results,  the  numerical  and  conceptual  significance  of  the 
approximation  should  be  understood  and  related  to  the  problem  at  hand. 

A  special  note  on  notation  is  appropriate.  In  past  RST  work  the 

parameter  1  was  an  Integer  where  1-1  for  the  Rayleigh  surface  wave  and 

o 

1  -  2,3,  etc.  for  whispering  gallery  surface  waves.”  In  previous  SWT  analyses 
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l  took  on  the  value  l  -  R  for  the  Rayleigh  wave  and  1-  WG1 ,  WG2,  etc.  for 
whispering  gallery  waves. *  Throughout  this  paper  we  shall  use  both  the  RST 
and  SWT  definition  of  l  and  imply  from  the  previous  work  chat  1  •  1  is 
equivalent  to  l  ■  R  and  1-2  is  equivalent  to  l  -  WG1  for  tungsten 
carbide  spheres. 

4.2  The  backscattering  form  function  via  SWT  and  RST 

In  this  section  we  discuss  Che  relation  between  the  resonances  of  RST 

and  fluctuations  in  the  magnitude  of  the  backscattering  form  function  of  a 

fluid  loaded  elastic  body  as  ka  is  varied.  This  relation  is  examined  through 

Che  specific  case  of  backscattering  from  a  tungsten  carbide  sphere  in  water. 

Section  IA  is  devoted  to  comparing  Che  PWS  with  the  SWT  approximations  for  the 

3  5  6 

backscattering  form  function  of  the  tungsten  carbide  sphere.  ’  ’  Besides 
this  comparison  several  intermediate  SWT  results  are  discussed.  Some  of  these 
intermediate  results  introduce  terminology  and  ideas  useful  in  Section  IB. 
Section  IA  serves  the  dual  purpose  of  partially  checking  the  SWT  of  Ref.  5  and 
laying  the  groundwork  for  the  discussion  of  RST  and  form  function  fluctuations 
in  Section  IB.  A  comment  is  appropriate  on  Che  computer  programs  used  in 
generating  many  of  the  figures  of  this  section.  In  evaluating  SWT  results,  we 
use  subroutines  which  calculate  Bessel  functions  of  complex  order  and  argument 
(see  Ref.  5  for  further  discussion  of  subroutines).  These  subroutines  have  an 
error  of  less  than  (0.00001  +  10. 00001)  when  the  modulus  of  the  Bessel 
function  order  |v|  is  greater  chan  3.  With  |v|  <  3  the  error  increases 
rapidly.  SWT  results  for  the  tungsten  carbide  sphere  below  ka  *  20  are 
affected  by  this  inaccuracy.  In  some  figures  SWT  results  below  ka  -  20  are 
plotted  since  the  qualitative  features  agree  with  PWS  results  even  though 
quantitative  accuracy  cannot  be  assured.  The  PWS  calculations  have  no  such 
limitations.  The  maximum  ka  range  plotted  for  SWT  results  is  10  <  ka  <  80  and 


100 


for  PWS  results  is  1  <  ka  ^  80.  The  methods  used  to  obtain  SWT  results  are 
similar  to  those  discussed  In  Ref.  5.  The  material  parameters  used  for 
tungsten  carbide  ware:  density  -  13.80  g/cm  ,  longitudinal  wave 

speed  -  6.860  x  10'  ca/s,  sheer  wave  speed  ~  4. 185  x  10*  cm/s.  Those  for 

3  5 

water  were:  density  -  1  g/cm  ,  longitudinal  wave  speed  •  1.4760  x  10  cm/s. 


A.  PWS  and  SWT  form  functions 

When  a  plane  acoustic  wave  Is  scattered  from  a  fluid- loaded  elastic 

sphere  and  wa  use  the  coordinate  system  shown  In  Fig.  1  the  total  pressure  In 

3  5 

the  far  field  of  a  fluid- loaded  sphere  can  be  written  as  the  real  part  of  ’ 

tk_  a.ikr 

PT(x,0)  -  •  +  e(x’9)  (1) 


where  x  =  lea.  On  the  right-hand- side  the  first  term  represents  a  unit 
amplitude  incident  plane  wave  and  the  second  term  the  scattered  pressure.  The 
f  in  this  expression  is  a  form  function  descriptive  of  scattering  when 
r  »  ka2. 

The  PWS  backscattering  form  function  f(x,iT)  -  f^WS(x,ir)  Is2 


B_(x) 


*•"  ■  £  i  SS 

n-0  n 


where  Bq  and  Da  are  3x3  determinants  given  in  Ref.  3.  The  complex 

resonances  of  RST  fulfill  the  condition  D  (x  „)  ■  0.  It  Is  important  for  the 

a  ni  r 

discussion  in  Section  II  to  note  that  the  RST  resonances  (the  x  „ )  are 

nV 

"located"  in  a  complex  x  plane.  The  value  x^  corresponds  to  a  Rayleigh 
resonance  of  the  nC^  term  of  Eq.  (2)  while  the  values  x  with  l  >_  2 
correspond  to  any  of  the  many  whispering  gallery  resonances  of  the  nth  term. 
In  the  RST  of  a  sphere2,7  it  is  asserted  that  when  the  ka  of  the  sphere  equals 


Thm  r,  3,  p  coordinae*  syscaa  shown  above  vaa  u**d  j,a  arlciag  eh* 
partial  «*v«  s*rl«*  (?WS)  solution  eo  plan*  vav*  scattering  from  an 
•laaelc  sphere .  Th*  plan*  wav*  ij  aaaumad  co  b*  traveling  in  eh*  -z 
dir action. 


Fig.  1. 


,  f.  ■ 
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Ra( xq^)  chan  a  +  (1/2)  wavelengths  of  Cha  l  surfaca  wave  fic  on  cha 
aphara . 

Tha  SWT  nay  ba  used  to  find  an  alcarnaca  repraaancaclon  for  f.  In 
cha  SWT  ona  ravritaa  cha  PWS  for  cha  cocal  pressure  in  terms  of  a  contour 
lncagral  in  a  complex  y  plana. ^  Tha  contour  uaad  can  ba  daformad  Co  surround 
cooplax  polas  whose  raalduaa  given  cha  contribution  from  Franz,  Rayleigh,  and 
whispering  gallery  waves.  Furthermore,  Cha  contour  passes  through  saddle 
points  whose  evaluation  gives  contributions  from  specular  reflection  and 
transmitted  bulk  waves.  After  a  considerable  effort  one  can  wrlta  the  far 
field  prassure  as  shown  in  Eq.  (1)  with  f  a  f3^  where  ka  »  1  and  f3^ 
can  be  separated  into  the  sum  of  several  contributions.^  Hera  we  consider 
only  the  contributions  from  the  specular  reflection,  Rayleigh  wave,  and  cha 
slowst  whispering  gallary  wave  which  we  label  f  ,  f  ,  f  respectively. 

The  analysis  of  Raf.  5  gives  Che  essential  details  of  tha  derivation, 
interpretation ,  and  numerical  implementation  of  tha  results  for  fg,  f^,  f^G^ . 
We  will  only  briefly  summarize  Chair  form  for  tha  case  9  «  t;  i.a., 
backscatterlng . 

When  9  ■  tt  the  specular  reflection  contribution  fg  has  the  form 

'  -  as«  -0)(l!>  «‘123‘  (3> 

S 

where  Rg^  is  the  coefficient  of  reflection  from  the  front  of  the  sphere  and 
s 

(-2x)  is  the  phase  of  the  specular  reflection  relative  to  a  ray  traveling  in 
the  liquid  to  and  from  a  reference  point  corresponding  to  the  spheres  center. 
The  full  expression  for  R^^  ,  which  was  used  in  the  calculations,  is  derived 
in  Appendix  A  of  Ref.  5. 

The  contributions  f^  and  fyG^  *r*  found  from  a  residue  analysis 
using  the  Rayleigh  pole  and  the  appropriate  whispering  gallery  pole 


'-V-V.V- 


*.  %  •_  *•  L'.  *- 
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VWG1*  expressions  for  fR  and  f.^  have  cha  sama  form  and  at  8  ■  tt 

they  baeoma 


f^Cx.TT)  ■  -Gj  • 


inl  “2(7r“9i)0i  r  -imr  i2TO(Gt£44*) 

«  l  a  a  a 

B-0 


(4) 


whara  2,  equals  R  or  WG1  for  the  Rayleigh  or  whispering  gallery 

contribution  and  wa  have  used  the  substitution  +•  iS^.  We  note  chat 

9Z*  ^2,'  ®2.»  a£»  Gl  ar*  a22  Unctions  of  x  and  that  f^  la  actually  cha 

SWT 

appropriate  form  for  any  whispering  gallery  contribution  to  f  ,  Equation 
(4)  can  be  interpreted  In  terms  of  surface  waves  in  which  case  tha  variables 
and  summation  have  tha  following  physical  significance^:  G^  is  a  complex 
amplitude  factor  that  accounts  for  cha  coupling  efficiency  of  Che  2,ch 
surface  wave  onto  Che  sphere  as  well  as  the  effect  of  axial  focusing,  5^  Is 
Che  local  angle  of  incidence  whara  Che  surface  wave  and  incident  acoustic  wave 

a 

are  phase  matched,  the  cerms  exp<irii)  and  axpf-ZCTT-Q^)  S^]  are  the  phase 
delay  (relative  to  the  same  reference  wave  as  used  for  Eq.  (3))  and 
attenuation  due  co  radiation  damping  the  first  time  the  wave  radiates  energy 
in  the  backscatterlng  direction,  the  sum  over  m  accounts  for  the  surface 
wave  circumnavigating  the  sphere  an  infinite  number  of  times,  the 
exp(i2Tm(a^  +  Jj)]  and  axpI-^TimB^]  are  the  propagation  phase  delay  and 
attenuation  of  the  surface  wave  for  m  circumnavigations  the  sphere,  the 
expj-innr]  accounts  for  phase  shifts  due  to  caustics  at  0-0  and  9  ■  ir. 

We  note  for  future  reference  that  the  results  above  show  that  the  propagation 
phase  delay  for  one  complete  crip  around  Che  sphere  is  exp[i2ir(a^  +  4)1* 

This  Implies  chat  is  related  co  Che  number  of  cycles  of  the  surface  wave 

which  fit  on  the  sphere-.  For  instance  if  the  value  of  x  is  such  that 
equals  an  integer  n  than  n  +  (1/2)  cycles  of  the  surface  wave  would  fit  on 
Che  sphere.  Explicit  expressions  for  and  can  be  found  in  Ref.  5. 


[•] 


By  inspecting  Eq.  (4)  one  can  sss  that  the  backseat tering  amplitude  of 

th _ ___  „th 


the  d'“‘  term  of  ths  1"“  surface  wave  is 

-2  Or -2irm84 


-  |fl£|  • 


(5) 


The  ability  to  calculate  the  coupling  efficiencies  G^  and  beckscattering 

amplitudes  A  can  prove  useful  in  transient  scattering  situations  where 

one  may  want  to  determine  the  physical  origin  of  various  echoes  seen  in  the 
3  9 

scattered  return.  *  The  coupling  efficiencies  may  also  be  useful  in 

investigating  more  complex  shapes  such  as  cylinders  with  hemispherical 

endeaps.  With  these  types  of  applications  in  mind  Figure  2  gives  |G_| , 

|G„_, | ,  A_  A_  for  a  tungsten  carbide  sphere.  Throughout  the  ka  range 

WG1  Q  jK  0  y WG1 

of  Fig.  2,  J  f s  J  3  0.9693  (see  Ref.  5)  and  comparison  with  AQ  R  from  Fig.  2 
suggests  that  fR  will  strongly  influence  the  total  form  function.  This 
supposition  proves  to  be  correct  in  what  follows* 

An  alternate  expression  for  f^(x,rr)  can  be  found  by  applying  a 
relation  used  in  the  analysis  of  optical  Fabry-Peroc  resonators^’ ^ 


CD 


l 


m-0 

Using  Eq. 


(6)  we  can  write  Eq.  (4)  as 


(6) 


-Gz  expI-ZOr-Q^Bj  +  in4] 
fi(x,7r)  ”  {1  +  e*p[-2irgA  +  i27r(aJl"+  *s)]> 


(7) 


This  closed-form  result  is  more  convenient  for  steady  stats  computations. 
Examination  of  the  form  of  Eq.  (7)  shows  that  if  the  x  value  of  the  sphere 
is  such  that  Re  equals  n  (where  n  is  an  integer)  the  magnitude  of 

the  denominator  is  close  to  a  minimum  and  we  have  a  resonance  type  behavior. 

We  have  discussed  previously  the  idea  that  when  ■  n  we  have  n  +  (1/2) 


V.T  V 


/  ..••• 
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Th«  Rayleigh  and  slowest  whispering  gallary  coefficient a  CR  and 
GVG1  f0r  couPlla*  oaCo  »  cungatan  carbide  sphere  ara  shown  (of.  zq. 

4).  Alao  shown  ara  cha  backscactering  amplitudes  A  and  A 

0,3  a0,WG1 

tor  cha  first  cime  cha  Rayleigh  and  whispering  gallery  waves  return 
energy  (cf.  2q.  (3)).  The  different  scales  needed  in  chis  ploc  ara 
an  indication  of  cha  ouch  larger  coupling  and  backscattered  amplitude 
of  cha  Rayleigh  wave  in  Che  ease  of  a  Cungstan  carbide  sphere. 


g.  2. 
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cycles  of  eh«  2.'“  surface  wave  on  Che  sphere;  a  physical  resonance  as 

defined  In  che  Introduction.  One  way  to  understand  the  resonance  when 

n  +  (1/2)  cycles  fit  on  the  sphere  is  through  the  use  of  Eq.  (4).  Each  time 

che  surface  wave  completes  a  trip  around  the  sphere  (each  of  the  individual 

terms  in  sum  of  Eq.  (4))  it  is  in  phase  with  contributions  from  previous 

circumnavigations  ( the  other  terms  of  che  sum)  .  In  phase  addition  occurs  when 

n  -f  (1/2)  cycles  fit  on  the  sphere  instead  of  n  as  would  be  the  case  for  a 
2 

cylinder  because  of  che  phase  advance  due  to  caustics  at  9  -  0  and  9  »  tt. 
To  see  che  resonance  behavior  of  Eq.  (7)  we  have  plotted  |f  (x,it)|  for  a 

IV 

tungsten  carbide  sphere  in  Fig .  3 .  The  resonance  behavior  is  clearly  evident . 

We  have  also  plotted  as  vertical  dotted  lines  the  real  part  of  the  RST 

Rayleigh  resonances  Ra(xa^)  for  5  £  n  £  28.  The  values  of  Re(x^)  were 

taken  from  Table  I  of  Ref.  8.  The  RST  concept  of  n  +  (1/2)  cycles  of  the  1th 

surface  wave  fitting  on  che  sphere  when  the  ka  of  the  sphere  equals  Re(xn^) 

is  confirmed  to  within  che  resolution  of  this  figure.  Indeed  we  will  show  in 

Section  II  that  the  approximation  inherent  in  Che  use  of  this  concept  is  not 

detectable  in  any  of  che  figures  of  this  section  for  the  particular  case  at 

hand.  Therefore,  for  the  remainder  of  this  section  we  use  che  Re(x  ,  )  to 

ni 

indicate  when  n  +  (1/2)  cycles  of  the  Rayleigh  wave  fit  on  the  tungsten 
carbide  sphere . 

The  results  to  this  point  can  now  be  used  to  obtain  curves  for  che 

magnitude  of  the  beckscattaring  form  function  of  a  tungsten  carbide  sphere 

PWS 

immersed  in  water  as  a  function  of  ka.  The  PWS  result  for  |f  (x,ir)|  =  f. 

D 

was  calculated  using  computer  programs  referenced  in  previous  work.^  We  first 

compare  f^  with  a  SWT  result  using  only  the  specular  reflection  and  Rayleigh 
SWT 

contributions  to  f  and  given  by  |f  (x,ir)  +  f  (x,ir)|  =  f  .  Figure  4 

5  R  SR 

shows  f^  in  the  range  1  £  ka  £  80  and  fgR  in  che  range  10  <  ka  <  80. 


In  this  figure  f^  Is  Che  absolute  value  of  cha  backacactering  fom 

function  for  a  cungsten  carbida  spbara  aa  calculated  using  cha  PWS 

jnrt  g  la  an  approxiaation  co  f.  obtained  using  cha  SWT  raaulcs 
3R  ® 

for  cha  specular  raflaccion  and  Rayleigh  surface  wave  contributions 
to  backacactering .  The  vertical  scala  for  haa  bean  shifted 

sinca  cha  evo  curves  are  aoc  easily  distinguished  otharvisa. 


We  have  offset  f  vertically  since  when  Che  curves  are  superimposed  they 
cannot  be  easily  distinguished.  In  Fig.  5  the  lea  scale  has  been  expanded  and 
f^  and  fgg  are  shown  without  offset.  This  facilitates  a  better  view  of  the 
agreement  between  the  two  curves  which  Is  typical  throughout  the  range 
10  <  ka  <  80.  These  figures  show  that  f_  correctly  reproduces  f.  except 
for  a  fine  structure  (superimposed  on  f^R)  due  to  other  waves. 

Though  fgg  is  the  final  result  needed  for  the  discussion  of  Section 
IIB  It  is  appropriate  to  further  test  the  SWT  results  of  Ref.  5  by  adding  in 
the  contribution  from  the  slowest  whispering  gallery  wave.  The  addition  of 
this  wave  gives  an  approximation  |fs(x,7r)  +  fR(x,ir)  +  fWG1(*»ff)l  =  fSRM(,. 
Figures  6  and  7  show  f^  and  £SRWG*  From  thasa  figures  one  sees  that  some 
of  the  fine  structure  abeent  from  f  has  been  recaptured.  The  increased 
agreement  Is  especially  obvious  in  Fig.  7.  Presumably  the  rest  of  the  fine 
structure  could  be  obtained  by  adding  In  contributions  from  other  whispering 
gallery  waves.  This  method  of  adding  one  surface  wave  at  a  time  to  the 
approximate  form  function  allows  ons  to  ascartain  tha  effect  of  each  surface 
wave  by  Inspection. 

It  is  appropriate  to  make  some  further  comments  on  these  results.  One 
can  understand  the  why  G_  is  spproxlmetsly  linear  in  ka  over  the  range 

tv 

displayed  in  Fig.  2  by  the  following  argument.  Tha  ka  dependence  results  from 
the  product  of  two  factors.  One  factor  Is  proportional  to  /ka  because  of 
axial  focusing.  The  other  /ka  factor  can  be  attributed  to  the  Increased 
coupling  efficiency  of  surface  waves  with  increase  of  ka .  This  last  /ka 
factor  can  be  obtained  usirg  energy  conservation  arguments  similar  to  those  of 
Ref.  13.  One  can  flitd'  this  /ka  factor  due  to  Increased  coupling  efficiency 
in  cylindrical  SWT  results.  Casting  the  SWT  for  cylinders^-4  In  s  form  similar 
to  that  for  spheres  givjs  a  coupling  factor  for  Rayleigh  waves  with  s  form 


10  20  30  40  30  80  70  80 

ka 


This  Uigur*  la  similar  eo  Fig.  4  axcapc  Chat  f^^,  Includes  eha 
coacrlbucioaa  co  baclcscactaring  from  eh*  slows c  whisparing  gallary 
wav*  as  wall  as  eha  contributions  from  spacular  raflactlon  and  eha 
Raylaigh  wava  lncludad  la  f_.  tha  Inclusion  of  eha  whlsparing 
gallary  wava  in  f  adds  soma  fin*  scruceur*  eo  f_  chua  giving 


baccar  agraamant  with  f.  . 


11 


ka 


Shown  on  an  expanded  seal*  la  part  of  f.  and  ff  from  Fig.  6 

0  5RWG 

but  without  vertical  seal*  off sat.  As  In  Fig.  5  this  allows  a  bactar 
vlaw  of  eha  agraamant  baevaan  f.  and  cha  SWT  result. 
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Ilka  G^  but  divided  by  /lea  dua  co  Cha  face  there  Is  no  axial  focusing. 

This  implies  an  approximate  /lea  dependence  for  the  cylindrical  counterpart 

to  GR.  The  surface  wave's  phase  velocity  dispersion  may  alter  this  ka 

dependency  for  coupling  efficiency.  The  above  arguments  are  therefore  only 

valid  in  ka  regions  where  the  surface  wave  of  interest  is  weakly  dispersive. 

Finally  we  should  discuss  why  our  synthesis  of  f^  in  Figs.  6  and  7 

works  as  well  as  it  does.  It  was  not  necessary  to  include  transmitted  bulk 

waves  such  as  Chose  described  in  Ref.  12  for  a  fused  silica  sphere.  This  is 

because  for  tungsten  carbide,  the  contributions  due  to  transmitted  waves  are 

estimated  to  be  negligible  in  comparison  to  chose  due  to  specular  reflection, 

the  Rayleigh  wave,  and  cha  leading  whispering  gallery  waves  near  resonance  in 

a 

the  ka  region  of  interest.  It  is  significant  in  this  regard  that  the  density 

of  tungsten  carbide  is  much  greater  than  that  of  water. 

The  reformulation  of  Che  m  series  of  Eq.  (4)  into  the  closed  form  of 

Eq.  (7)  is  similar,  in  principle,  to  Felsen's  "Hybrid  Synthesis"  technique.^ 

The  form  of  Eq.  (4)  is  useful  for  transient  analysis^  while  chat  of  Eq.  (7)  is 

•  « 

more  convenient  for  describing  the  steady  state  scattering;. 

B.  Interpreting  RST  using  f„^ 

The  physical  nature  of  is  relatively  simple.  The  discussion  of 

this  section  elucidates  Che  physical  interactions  which  produce  the  structure 

of  f^  and  in  doing  so  allows  an  excellent  opportunity  for  examining  RST. 

The  main  topics  of  discussion  in  this  section  are  the  significance  of  the 

phase  between  the  speculer  reflection  and  first  Rayleigh  wave  contribution  in 

determining  the  structure  of  f  and  the  relation  between  Che  Rayleigh 

resonances  x  ,  of  RST  and  the  maxima  and  minima  of  f  „.. . 

nl  SR 

The  phaae  between  the  specular  reflection  and  the  first  Rayleigh 
contribution  can  be  determined  from  Eq.  (3)  and  the  m  ■  0  term  of  Eq.  (4)  to 


be  n:  =  2x  +  nR  +  arg(GR)  where  arg(GR)  is  defined  as  arccan(InsGp/ReGR) . 

In  Figs.  8  and  9  both  f^R  and  r^®  =  Modulo  2r  are  plotted  as  functions 

of  ka.  Figure  9  is  an  anlargeaent  of  the  region  35  ^  k  <_  63  of  Fig.  8.  In 

addition  we  have  Indicated  values  of  Ra(x  ,)  for  3  <  n  <  28  by  the 

III  —  — 

position  of  vertical  lines.  These  lines  are  labeler  with  the  appropriate 
value  of  n  at  the  cop  of  the  figures. 

The  vertical  lines  allow  one  to  determine  r^®  at  the  point  where  the 
Rayleigh  wave  contributions  all  add  in  phase.  By  examination  of  ,  >.3  and  9 

one  can  make  several  statements  about  Che  and  its  relation  to  the 

structure  of  f^.  We  see  that  the  value  of  n^®  at  the  point  where 
n  +  (1/2)  Rayleigh  wavelengths  fit  on  the  sphere  is  uot  the  3ame  for  all  n. 
Furthermore,  Che  value  of  at  these  points  dictates  what  kind  of 

structure  f„  will  have  in  the  local  region.  For  instance  n,®  •  0  when 

5K.  1 

n  •  13,28  and  f^  has  a  maximum  right  at  these  points.  Likewise  n^®  -  it 
whoa  n  ■  4,21  and  fgR  has  a  minimum  at  these  points.  However,  it  is 
obvious  from  Figs.  3  and  9  that  in  most  cases  n^®  does  not  equal  0  or  tt 
and  the  maximum  and  minlmums  of  f^  do  not  coincide  with  Che  in  phase  adding 
of  the  Rayleigh  contributions. 

The  ideas  and  results  to  this  point  can  be  used  to  reexamine  and 
qualify  the  recent  RST  analysis  of  the  backscatterlng  form  function  of  a 
tungsten  carbide  sphere.  First,  the  Rayleigh  waves  are  responsible  noc  only 
for  the  quasinulls  in  the  (summed)  form  function  but  are  also  responsible  for 
many  of  chs  maximum ■ ,  Indeed,  the  interference  of  ths  specular  reflection  and 
Rayleigh  contributions  are  rasponsibls  for  tha  overall  structure  of  the  form 
function  in  a  much  more  encompassing  way  chan  the  RST  discussion  implies. 
Furthermore,  the  statement  that  the  quasinulls  of  ths  form  function  are  "at 
the  roots  x^,"3  i»  onl?  tru«  ior  n  ■  4,21  which  points  n^®  *  it.  Last, 


The  phase  n1  between  the  specular  reflection  and  first  Rayleigh 

wave  contribution  to  backsca tearing  is  shown  in  cha  upper  part  of  cha 

figure.  In  cha  lower  part  f^  from  Fig.  4  la  given.  As  in  Fig.  3 

cha  vertical  Unas  indicate  the  position  of  the  real  part  of  cha 

Rayleigh  resonances  Re  x  ,  where  a  is  given  at  cha  top  of  cha 

ni 

figure,  to  within  the  resolution  of  chls  graph  Che  vertical  lines 
allow  one  to  find  the  phase  between  cha  first  Rayleigh  waves 
contribution  and  specular  reflection  when  a  +  (1/2)  cycles  of  the 
Rayleigh  wave  fit  on  Che  sphere;  i.e.,  at  resonance.  Note  chat  f 

wJC 

does  aot  always  have  a  maximum  or  minimum  at  Che  Re(x  )  values. 


Figs.  8  and  9  imply  that  tha  Rayleigh  wave  may  be  responsible  for  a  local 
minimum  or  maximum  in  th«  form  function  even  though  the  minimum  or  maximum  may 


not  be  at  one  of  the  xq^  of  the  Rayleigh  waves.  This  implies  that  some  care 
may  be  needed  in  using  the  labeling  in  Fig.  4  of  Ref.  3  to  understand 
fluctuations  in  The  labeling  may  correspond  to  the  numerical  values  of 

x^ ;  however,  to  understand  tha  effect  of  a  resonance  on  the  form  function, 
the  coupling  efficiency  onto  the  sphere  as  well  as  the  interference  between 
the  specular  reflection  and  the  surface  waves  must  be  taken  into  account. 

4.3  Resonance  conditions  from  RST  and  the  SWT 

In  this  section  the  condition  for  resonances  given  by  the  SWT  (the 

"physical  resonances")  are  compared  with  that  of  RST  ("RST  resonances")  .  The 

RST  and  SWT  treatments  of  resonance  differ.  Tha  RST  x  „  are  found  by 

ni  7 

allowing  x  to  cake  on  complex  values  while  n  remains  an  Integer;  i.a.,  RST 
resonances  are  located  in  a  complex  x  plane.  In  the  SWT,  however,  x 
remains  real  and  one  allows  n  to  take  on  continuous  complex  values  (i.e., 

Che  \>  values)  and  the  poles  are  located  in  a  complex  \>  plane.  It  is  the 
poles  in  this  complex  \>  plane  which  lead  to  the  idea  of  damped  surface 
waves.  In  Sec.  IIA  we  use  Hq.  (7)  to  examine  the  SWT  and  RST  resonance 
conditions  and  then  in  Sec.  IIB  we  arrive  at  the  same  results  through  a  more 
formal  channel. 

A.  A  first  look 

Using  again  the  definition  +  iB^ ,  Eq.  (7)  can  be  written  as 

-Gl  exp[in£  -  j 

V*’70  -  - U~"exP(i2Trva)] -  (8) 

This  equation  resulted  from  the  SWT  and  it  is  obvious  that  the  resonance 
behavior  is  related  to  the  complex  value  of  v,  where  we  must  remember  that 


v  la  a  function  of  the  real  ka  5  x.  By  varying  x  we  can  locata  tha 
valuta  vj  in  tha  conplax  v  plana  whara  ■  n  +  ifi^  and  n  +  (1/2) 

wavelengths  of  tha  surface  wave  fit  on  tha  sphara  (cf.  Sac.  I).  Tha  valuo  of 
x  at  which  Ra(v°  )  ■  n  art  tha  rosonanca  conditions  given  by  tha  SWT. 

These  x  are  raal  valued  and  are  designated  as  x*. 

To  examine  tha  RST  rasonanca  behavior  of  f  wa  can  expand 
locally  around  tha  point  x£  giving 

\>l  *  a  +  ifi“  +  v£(x-x°)  +  .  .  .  (9) 

where  -  Reu|  +  i  Imv^  and  the  prime  denotes  differentiation  with  respect 
to  x  evaluated  at  x™.  Substituting  Eq.  (9)  Into  Eq.  (8)  gives 


exp£in^  -  28jl(ff  -  0£)] 

fl(x,7r)  3  {l  -  exp(i2m^(x'  -  "(x“  +  Ax)  +  iV4) TT  (10) 

whare  Ax  -  -8^Q(ImvJ|)/jv^|2  and  »  ZS^Rev^)/!  vjJ  2 .  Note  that  Ax  is 

negative  real.  This  approximation  is  valid  near  resonance  when 

lv£(x  -  x“)2|  «J  y^(x  -  x^) | .  From  Eq.  (10)  wa  can  look  at  the  RST  resonance 

condition  by  allowing  x  to  cake  on  complex  values.  The  pole  in  f^  is 

located  at  x  »  (x°  +  Ax)  -  LH^  ■  xq^. 

Recall  that  the  usual  RST  resonance  condition  is  x  •  Re(x  „).  We  see 

nr 

that  this  differs  from  the  SWT  result  (x  *  x°)  by  an  error  Ax.  The 
magnitude  of  the  error  obviously  depends  on  Che  pole  in  Che  complex  v  plane 
associated  with  Che  resonance.  We  have  calculated  numerical  values  of  Ax 
for  the  Rayleigh  pole  and  one  whispering  gallery  pole  in  both  tungsten  carbide 
and  fused  silica.  The  values  of  ax  at  ka  =  40  ranged  from  Ax  =  -0.0007 


for  che  whispering  gallery  pole  of  tungsten  carbide  to  Ax  a  -0.1  for  the 
Rayleigh  pole  of  fused  silica.  The  Ax  shift  for  the  Rayleigh  pole  of 
tungsten  carbide  was  less  chan  -0.003  throughout  the  ka  range  of  the  figures 
of  Sec.  I  and  is  not  detectable  for  any  of  the  plots  of  Chat  section. 

For  a  given  value  of  x,  the  phase  velocity  c  of  a  surface  wave  may 

be  calculated  using  the  SWT  from  the  location  of  the  v  pole.  The  result, 
which  follows  from  Eq.  (30)  of  Ref.  5  and  the  phase  matching  condition 
sin0£  ■  c^/c,  is  Cg/c  ■  */(<*£  +  (1/2)].  Letting  x  •  x“  so  chat  ■  n 

gives 

cz/c  -  xJ/[n+(l/2)]  (11) 

Since  x™  ■  Re(x^)  -  Ax,  we  also  have 

c4/c  -  [Re(xnjl)  -  Ax]/[n  +  (1/2)]  (12) 

Equation  (12)  allows  one  co  properly  calculate  surface  wave  phase  velocity 

from  che  RST  resonance  condition. 

2  16 

Previous  work  *  in  RST  relating  che  v  and  x  complex  planes  for 
cylinders  and  spheres  started  in  Che  x  plane  and  then  carried  out  an 
expansion  of  x^  equivalent  to  the  expansion  of  In  Eq.  (9).  By 

examination  of  these  expansions  one  finds  an  implicit  assumption  in  Ref.  2 
that  ddmx^J/dv  «  0.  This  is  analogous  to  assuming  Imv^  -  0  in  Eq.  (9)  in 

which  case  Ax  would  equal  zero.  One  should  note  that  Che  relation  between 

17  18 

the  complex  v  and  x  planes  was  previously  examined  ’  for  applications 
of  Che  SWT  to  quantum  mechanical  scattering,  especially  Ref.  18,  pp.  106-113. 
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B.  A  more  formal  approach 

Tha  results  at  S«c.  XIA  can  b«  darlvad  using  cha  ganaral  expressions 

PWS  SWT  PWS 

for  Cha  fora  function  f  and  f  .  Tha  expressions  used  for  f  and 

SWT 

f  in  Chis  section  are  lass  spacific  than  glvan  in  Sec.  I  sinca  detailed 

forms  are  not  necessary  in  cha  derivation.  Tha  approach  of  this  section  is 

similar  to  chat  described  by  Newton  for  quantum  mechanical  scattering.*7 

PWS 

Th«  PWS  rasult  f  can  bm  vrlctcn  as 
00 

fPWS(x,9)  -  l  A  (x)  -P  (y)  (13) 

n-0  n  n 

where  u  -  cos 0,  PQ  is  the  Legendre  polynomial,  and  Aq(x)  is  a  complicated 

3 

function  of  x;  A  diverges  at  the  complex  x  .  of  RST. 

a  nJc 

SVTT 

As  discussed  in  Sec.  1,  f  can  be  written  as  a  sum  of 

contributions  from  specular  reflection,  transmitted  bulk  waves  and  surface 

5 

waves*  (For  completeness  one  must  also  include  a  possible  contribution  from 
a  "background  integral"  and  that  due  to  forward  diffraction  in  the  case  of 
near  forward  scattering.)  One  can  write 

fSWT(x,9)  -  f.  +  l  t  (14) 

Z-R.WG 

SWT 

where  f^  includes  all  contributions  to  f  except  those  due  to  the 
Rayleigh  and  all  whispering  gallery  waves  which  are  accounted  for  by  the  sum. 

For  purposes  of  the  present  discussion  it  is  convenient  to  rewrite  the 
sum  of  Sq.  (14)  in  slightly  more  detail.  Using  Eqs.  (22)  and  (23)  of  Ref.  5 
we  can  write 


fSMT(x.9) 


fl  + 


l 

l- R,WG 


B  (x) 

V 


Pv  (-U) 
i 


sinTrv. 


(15) 


where  P  is  now  a  Legendre  function  of  tha  first  kind.  Whan  x  is  larga 

vt 

_  .  .PWS  .SWT  _ 
we  nay  taka  f  ■  f  .  Thus  wa  hava 


LV*)Fa<U)  '  ‘l*,J 


»v  (*)  Pv  <-P) 

v2 _ vi 


n-0  “  *  2-R.WG  sinirvl 


(16) 


Now  multiply  both  sidas  of  Eq.  (16)  by  Pa,(y)  ware  n'  is  a  nonnagative 

19  20 

intagar  and  integrate  from  y  -  -1  to  y  -  1.  Using  tha  relations  ’ 


'-1 


P  , (y)Pn (y)du 
a  n 


(2n+l)  ^nn' 


(17) 


siniTv„ 


L1  Pn-(w)Pv/*u)du  "  * 


n  —  v„'  it  (vl-a,)(vji+n,+l) 


(18) 


wa  find  that 


A  (x) 

Cl 


l 


l  .  — 

(v.-n) (v.+n+l)  n 


2-R,WG  VV1 


(19) 


,  2n-t-l . 


(^~)  f  p  (u)du 
*•  J  i  in 


'-1 


(20) 


where  C  -  (2afl)B(i  /ir  depends  on  x  and  6  ,  in  Eq.  (17)  is  the 

nn' 

Kroneckar  delta.  It  is  apparent  from  the  form  of  tha  summation  that  A  (x) 

a 

(and  therefore  tha  form  function)  exhibits  resonance  behavior  at  x  near  x”, 
which  is  the  resonance  condition  from  the  SWT.  Recall  that  Re(v(J)  -  n  when 
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x  ■  x°.  (There  may  also  be  some  resonance  behavior  due  co  Eq  but  wc  will 
not  address  this.)  Applying  the  expansion  given  by  Eq.  (9)  to  (v^  *  n)  In 
Eq.  (19)  gives 


VX) 


1-R.MG  Vi(Vl+n+1Hx"(xi+Ax)+il5rJl)  +  E° 


(21) 


where  Ax  and  are  as  defined  below  Eq.  (10).  We  now  see  poles  in  Che 

complex  x  plane  at  the  values  x^  -  (x£  +  Ax)  -l^  as  before  and  the 

discussion  and  cautions  subsequent  co  Eq.  (10)  apply. 

This  paper  has  reexamined  some  of  Che  central  concepts  of  RST.  Though 

the  specific  case  of  a  fluid- loaded  sphere  was  used,  the  results  have  more 

general  ramifications.  The  usefulness  and  power  of  Che  SWT  for  understanding 

acoustical  scattering  problems  has  been  illustrated.  It  is  appropriate  Co 

note  that  Che  present  research  was  completed  on  che  centenary  of  Rayleigh' s 
21 

original  analysis  of  surface  waves  on  Che  flat  surface  of  an  elastic  solid. 
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ZNU  and  ZKA 

The  program  ZNU  calculates  Cha  SWT  poles  of  a  tungsten  carbide  sphere 
In  Che  complex  v  plana  given  Che  real  lea  value  of  Che  sphere.  Parcs  of  Che 
program  were  adapted  from  a  program  wrlccen  by  Brad  Brim  as  pare  of  a  master's 
degree  project  (Ref.  21  of  Chapter  2).  The  program  requires  Chat  you  Input  a 
guess  (GUESS)  of  che  pole  position.  It  also  asks  for  Che  region  in  Che 
complex  plane  where  it  is  supposed  eo  look  (HEIGHT) ,  the  precision  required  in 
pole  location  (ERROR),  and  che  ka  value  of  Che  sphere.  The  program  assumes  a 
complex  value  for  ka  buc  one  need  only  sec  che  Im  ka  ■  0  co  find  SWT  poles. 

The  program  uses  che  winding  number  cheorrn  Co  find  che  poles.  There  Is  a 
commentary  (written  by  Brad  Brim)  within  che  program  on  che  winding  number 
theorem.  Within  this  commentary  Che  other  parameter  you  are  asked  Co  input 
(Npnts)  Is  discussed.  The  program  finds  Cha  SWT  poles  by  seeking  the  zeros  of 
Che  function  Dy  defined  in  Appendix  A  of  Chapter  3 . 

The  subroutine  OLVER  (which  is  called  within  che  program)  is  given 
later  in  this  appendix.  This  subroutine  limits  the  accuracy  of  finding  che 
SWT  poles  Co  no  better  than  0.000001  +  10.00001.  Furthermore,  this  accuracy 
is  only  applicable  under  che  conditions  given  in  Section  4.2. 

This  program  was  used  to  generate  che  numbers  found  in  Table  I  of 
Chapter  3.  That  table  can  be  used  to  check  that  the  program  is  working 
properly.  If  used  for  other  isotropic  spheres  the  material  parameters  must  be 
changed . 

The  program  ZKA  (not  listed  here)  is  like  ZNU  except  that  it  locates 

the  zeros  of  D  (ka)  as  defined  in  Chapter  2.  Here  n  «  Re(NU)  is  a 
n 

specified  integer  and  Im(NU)  ■  0.  The  complex  ka  for  which  Dn(ka)  «  0  are 
obtained.  The  principal  modification  to  ZNU  is  in  the  subprogram  for  the 
function  f.  The  program  was  used  to  generate  Table  I  of  Chapter  2. 
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cccccccccccccccc 
c  c 

C  ZNU  C 

c  c 

cccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C  Tha  calling  program  and  subroutina  ZROOT  waa  originally 

C  writtan  by  Brad  8ria  o f  tha  S2  dapartaant.lt  baa  baan  altarad 

C  to  apply  to  the  acattaring  problaa  at  hand.ZRQOT  locataa  tha 

C  complex  zeros  and  polaa  of  tha  function  T  which  tha  uaar  dafinaa. 

C  At  praaant  F  ia  On  of  tha  fraction  Bn/Dn  diseuaaad  in  tha  text. 

C  This  varaion  of  tha  program  can  ba  uaad  to  find  tha  zaroa  in  tha 

C  coaplax  laada  (or  nu)  plana  givan  a  raal  valua  of  ka.This  program 
c  waa  uaad  to  find  tha  Rayleigh  polaa  givan  in  Tabla  1  of  chaptar  3 

C  .To  do  spharao  of  aatarial  othar  than  tungatan  carbida  you  muat 

C  changa  tha  notarial  paraaatara  in  tha  function  F. 

C  Tha  function  F  calla  tha  subroutina  OLVER  which  ia  givan 

C  in  a  saparata  part  of  this  appandix.I  make  no  claims  that  this 

C  program  aa  altarad  is  "efficient*  in  tha  sansa  of  computar  tins 

C  or  program  langth. 

C  Tha  program  includes  a  discussion  of  how  ZROOT  works. 

C 

cccccccceccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


£******************************* ************************* *******  ******** 

c*  * 

C*  Thia  ia  a  vary  ganaral  calling  program  asaociatad  with  tha  * 

C*  routina  "ZROOT".  It  handlas  all  tha  input  of  paramatars  and  * 

C*  output  of  raaulta.  * 

C*  * 

c********* ******* ******* ************************************************ 

c*  • 

C*  SEE  tha  subroutina  "ZROOT"  for  tha  rsquiramants  on  F(z)  and  tha  * 

C*  input  paramatars.  Containad  tharain  is  a  fairly  thorough  * 

C*  daacription  of  tha  algorithm  and  numarical  considerations.  * 

C*  * 

C*  **  **★★**#**★  +  w  nr  *  *  ******  **  ******  *  ***************  *******************  * 


REAL *8  RADII  ( 10)  ,X,X,  HEIGHT,  ERROR 

INTEGER  ORDER (10) 

COMPLEX* 16  ZEROS ( 10 ), GUESS, KA 
COMMON  KA 


1  WRITE  (1,10) 

10  FORMAT (///// , 1  REAL  PART  OF  GUESS ',T35,'  ') 

READ  (l,*)  X 


WRITE  (1,20) 


oooooonooooooonno 
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20  FORMAT ( '  IMAG  PART  OF  GUESS ',T35,'  ') 

READ  (1,*)  V 

GUISS  -  CMPLX(X,Y) 


WRITS  (1,30) 

30  FORMAT (/, '  DESIRED  HEIGHT' ,T33 , '  ') 

READ  (1,*)  HEIGHT 

WRITE  (1,40) 

40  FORMAT (/,'  DESIRED  ERROR ',T35,'  ') 

READ  (1,*)  ERROR 

ERROR  -  DMAX1 (ERROR, 1.0D-07) 

Npnta  -  100 
WRITE  (1,30) 

30  FORMAT (/,'  POINTS  PER  QUADRANT  SIDE',T3S,'  ') 
READ  (1,*)  Npnta 

WRITE ( 1 , 33 ) 

33  FORMAT (/, 'REAL  PART  OF  XA',T3S,'  ') 

READ(1, *)  X 

WRITE (1, 34) 

54  FORMAT (/, 'IMAG  PART  OF  KA',T33,'  ') 

READ(1, *)  Y 

KA-CMPLX(X,Y) 

WRITE  (1,60) 

60  FORMAT (5X, /////) 


CALL  ZROOT( GUESS, HEIGHT, ERROR, Npnts, 
+  N*«ro, ZEROS, RADII, ORDER) 


IF  (Nzaro  . EQ.  0)  THEN 


WRITE  (1,70) 

70  FORMAT (/,'  **********  THERE  WERE  NO  ZEROES  ********** 

ELSE 

WRITE  (1,80) 

80  FORMAT (/,'  **********  THE  ZEROES  WERE  ********** 


DO  100  1-1, Nzaro 

WRITE  (1,90)  I, ZEROS (I) , RADII (I) , ORDER (I) 
90  FORMAT (/ , / , 14 , 7X, 2E16 . 7 , 7X, E14 . 7 , 7X, 14 ) 

100  CONTINUE 

END  IF 

GO  TO  1 

END 
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SUBROUTINE  2 ROOT ( CENTR , HITS , ERROR , Npnta , 
+  Nxaro, ZEROS, RADII, ORDER) 


C********************************************************************** 

C *****  by  BRADLEY  L.  BRIM  Sinai  raviaiona  7/24/83  **** 

c ********************************************************************** 
C* 

C*  Thia  routina  aaarchaa  tha  aoaplax  plana  for  tha  roota  of  tha 

c*  function  T(z) .  F(z)  auat  ba  an  axtarnally  dafinad  COMPLEX* 16 

C*  function  aubroutina.  Tha  variabla  s  auat  alao  ba  COMPLEX*16. 

C*  Thia  routina  ia  aaaily  aodifiad  for  F(z)  •>  F(z,{Xn))  * 

C*  whara  (Xn)  ia  a  aat  of  n  paraaatara  indapandant  of  z. 

C* 

C******#*****«***********»#******M****#*************M*******»********' 


c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


Thia  routina  ia  baaad  on  tha  "WINDING  NUMBER"  thaoran,  found 
in  any  comp lax  variablaa  book. 

Thia  thaoram  baaically  aaya  ..... 

"Tha  intagral  around  any  aiaply  cloaad  contour  (avaluatad 
in  tha  clockwiaa  diraction)  of  tha  function  [F' (z)/F(z) ] 
will  ba  aqual  to  2*pi*i*(#P-#2) . " 

whara  ...  #z  ■  nuabar  of  zaroaa  within  tha  contour 
#F  m  nuabar  of  polaa  within  tha  contour 
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C*  | 

C*  | 

C*  | 

C*  | 

C*  | 

c*  I 

c*  I 

C*  | 

C*  CORKER  #3 
C* 

C* 

c* 
c* 
c* 
c* 
c* 

C*  This  routina  was  writtan  sspsoially  to  dstaraina  ths  location 

C*  o f  ALL  tha  roots  in  a  spacifiad  squara  ragion,  to  a  spacifiad 

C*  accuracy.  Othar  routinaa  nay  claim  to  ba  mora  officiant,  but 

C*  all  othar  conplax  root  saarch  routinas  this  author  knows  of  will 

C*  tarmlnata  aftar  finding  tha  location  at  only  tha  first  root  that 

C*  tha  routina  happans  upon. 

C* 

C*  Tha  paramatars  of  this  progran  ara  . 

C* 

C*  INPUT  :  CENTS 

C* 

C*  HITE 

C* 

C*  ERROR 

C* 

C*  Npnts 

C* 

c*  _ 

C*  OUTPUT  :  ZEROS 

C* 

C*  RADIX 

C* 

C*  ORDER 

C* 

C* 

C* 

C*  Tha  winding  nunbar  basis  tor  this  routina  allows  ona  to  find 

C*  tha  location  of  both  zaros  and/or  polas.  Tha  problam  with,  this 

C*  is,  if  a  ragion  is  saarchad  which  contains  an  aqual  nunbar  of 

C*  polas  and  zaros,  this  is  intsrpratad  as  a  ragion  with  no  roots. 

C*  Thus,  it  is  raconnandad  that  this  routina  ba  usad  with  functions 

C*  which  has  only  zaros  OR  polos.  Tha  axcaption  is  of  coursa  whan 
C*  ona  has  sons  knowladga  a  priori  of  ths  location  of  both  typas  of 
C*  roots  and  ragions  whara  thay  occur  mutually  axclusivaly. 

C*  This  routina  starts  out  by  calculating  tha  winding  nunbar  of 

C*  tha  givan  "initial  squara  ragion".  If  it's  non-zaro  tha  routina 
C*  naxt  looks  at  tha  quadrants  of  that  ragion  in  ascanding  ordar. 

C*  If  tha  winding  nunbar  is  zaro,  it  goas  on  to  tha  naxt  quadrant. 

C*  If  tha  winding  nunbar  is  non-zaro  and  tha  arror  is  small  anough, 

C*  tha  information  is  racordad.  If  tha  arror  is  too  larga,  this 


(C*18]  Cantor  of  tha  squara  ragion 
in  which  to  saarch  for  roots. 
[R*Q5]  Half  tha  langth  of  a  sida  of 
this  squara.  (Quadrant  haight) 
(R*08 ]  Tha  maximum  quadrant  haight 
of  tha  root  location (s) . 

(1*02 ]  Nunbar  of  points  par  quadrant 
sida  in  numarical  intagration. 

(10)  [C*16]  Cantor (s)  of  squara  ragion(s) 

containing  tha  root ( s ) . 

(10)  (H*08]  Acutal  quadrant  haight  of  tha 
corrasponding  ragion (s) . 

(10)  (1*02]  ordar  of  tha  root(s) 

+N  ->  Nth  ordar  ZERO 
-N  ->  Nth  ordar  POLE 
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C*  quadrant  becomes  tha  nav  "initial  aquara  raqion"  and  the  routina 
C*  will  atart  ovar.  whan  tha  routina  haa  dona  anough  "splits"  to 
C*  looata  a  root  to  tha  daairad  accuracy,  it  bagina  to  work  ita  way 

c*  back  to  largar  ragiona,  by  chocking  quadrants  of  previous  aplits. 

C*  Tha  variabla  Naplit  kaapa  track  of  how  many  times  tha  original 

C*  initial  aquara  ragion  haa  baan  aplit  up.  A  maximum  of  tvanty- 

c*  four  aplita  ara  allowad,  which  ahould  alwaya  ba  aufficiant. 

C*  Sinca  a  aubroutina  in  FORTRAN  aay  not  call  itaalf,  tha  arrays 

C*  CENTER , HEIGHT , QDRNTS , INTGRL, INSIDE  ara  uaad  to  atora  information 

C*  calaculatad  in  pravioua  aplita  ao  it  doaa  not  hava  to  ba  rodona. 
C* 

C*  Tha  function  (*" (z)/r(x) ] *dz  ia  approxiaatad  aa  . 

C* 

C*  { [T(z2) -F(Zl) ]/(z2-Zl] }  /  (0.5*(F(Zl)+F(z2) ] }  *  (Z2-Z1' 

C* 

C*  OR,  by  reducing 

c*  2*(i-a*{r(zi)/tr(zi)+F(z2) ]}) 

c* 

c*  Integrating  this  over  Npnts  point* 

C* 

C*  integral  «•  2  *  (Npnta  -  2*StJM) 

C* 

C*  where  z  «  (zl+z2)/2 

C*  StJM  •  tha  aua  for  each  z  F(zl)/ [F(zl)+F(z2) ] 

C* 

C*  This  author  found  it  both  numerically  accurate  and  officiant 

c*  in  time  to  calculate  tha  quantity  "SUM"  and  than  apply  tha 
C*  above  integration  formula. 

C*  Tha  aubroutina  "EXACT"  ia  called  if  F(.)  happens  to  ba  exactly 

C*  equal  to  zero.  That  fact  ia  printed  and  the  routine  atope. 

C*  The  subroutine  "BETWN"  ia  called  if  Fsua  ■  F(zl)+F(z2)  happens 

C*  to  be  exactly  equal  to  zero.  Thera,  F(z)  [z  as  defined  above] 

C*  ia  calculated.  If  F(z)  is  non-zero,  Fsum  is  replaced  by  twice 
C*  F(z)  and  tha  numerical  integration  continues.  If  F(z)  is  also 
C*  zero,  that  fact  is  printed  and  tha  routina  stopped. 

C*  Tha  reason  that  tha  routina  is  stopped  if  a  zero  is  found  is 

C*  because  numerical  problems  ara  vary  likely.  If  a  root  is 

C*  located  "near"  tha  contour  of  integration  (relative  to  the  step 

C*  and  quadrant  sizes)  than  this  root  will  contribute  approximately 
C*  half  its  normal  value  to  the  vinding  number.  Tha  contribution 
C*  is  almost  alwaya  graatar  than  half  of  normal  if  the  root  is  just 

C*  inside  tha  contour,  and  less  if  juat  outside.  Similarly,  for  a 

C*  root  located  near  a  corner  in  the  contour,  tha  contribution  will 

C*  ba  approximately  one  quarter  normal.  These  tvo  facts  tend  to 

C*  help  the  program  work,  evan  if  tha  root  is  close  to  tha  contour. 

C*  Tha  required  "Npnts"  will  vary  with  the  behavior  of  F(z)  It's 

C*  recommended  that  for  large  "HITE"  and  large  "Npnts"  one  not  uss 

C*  a  small  "ERROR".  Than,  it  savas  time  to  do  tvo  call  saquancas. 

C********************************************************************** 
COMPLEX* 16  ZEROS (10) , CENTR, dSIDE (4) ,C0RNER(4) , 

+  CENTER ( 24 ) , INTGRL ( 24 , 0:3), INSIDE (24 ,5) , 

+  Zcomr ,  Zand , F,  Fzl ,  Fz2  ,  Feum , SOM ,  WIND 

INTEGER  ORDER (10) , QDRNTS (24) , SIDE, HALT, WINDno, CROSS , QUAD, 

+  start, and 

RZAL*8  RADII (10) ,HEISHT(24) , HITE, ERROR, delta 
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EXTERNAL  T 


dSXDZ(l)  - 
dSXDE(2)  - 
dSXDS ( 3 )  - 

d3IDE(4)  - 

CORNER (1)  - 
CORNER ( 2 )  - 
CORNER  ( 3 )  » 
CORNER ( 4 )  - 

Nsplit  -  0 
Nzaro  ■  0 
CZNTER(l)  - 
HEIGHT (1)  - 
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(  0.000,-1.000) 
(-1.000,  0.000) 
(  0.000,  1.000) 
(  1.000,  0.000) 

(  1.000,  1.000) 
(  1.000,-1.000) 
(-1.-000,-1.000) 
(-1.000,  1.000) 


CSMTR 

HXTS 


1  Nsplit  -  Nsplit+1 

dslta  -  HEIGHT (Nsplit) /Npnts 


00  103  SIDE-1,4 

Zcomr  -  CENTER(Nsplit) +HEIGHT(Nsplit)  ‘CORNER (SIDE) 
Fz2  -  F(Zcorar) 

XT  (ABS (Fz2)  .EQ.  0.000) 

+  CALL  EXACT  (Zcomr) 

DO  103  HAL?- 1,2 

SOM  *  (0.000, O.ODO) 
start  -  l+(HALF-l)*Npnts 
and  -  HALF*Npnts 
00  101  LOOPS-atart , and 
Fzi  -  Fz2 

Fz2  -  F (ZCOCT1T+ LOOPS *dalta*dSIDE ( SIDE) ) 

Fsua  -  ?zl+Fz2 


IF  (AES (7z2)  .EQ.  O.ODO) 

+  CALL  EXACT (Zcomr+LOOPS*dalta*dSIDE (SIDE)  ) 

X?  (ABS  (Faun)  .EQ.  0.000) 

+  CALL  aETWN(Zcomr+(LOOPS-0.5)  *dalta*dSIDE (SIDE)  ,Fsum) 

SOM  -  SOM  +  Fzl/Fsum 

101  CONTINOE 

XNTGRL (Nsplit, (SIDE-1) *2+HALF)  -  2* (Npnts-2*SOM) 

102  CONTINOE 

103  CONTINOE 

XNTGRL (Nsplit, 0 )  -  XNTGRL (Nsplit, 8) 


NINO  -  (0.000,0.000) 

DO  200  LOOPS-1,3 

WIND  -  WINO+INTGRL( Nsplit , LOOPS) / (0 . 0D0 , -8 . 28318S308D0) 
200  CONTINOE 

WINDno  -  NINT ( REAL ( WIND ) ) 


I?  (WINDno  .EQ.  0)  GO  TO  6 
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WRITS  (*,10)  Naplit, WIND, CENTER (Naplit)  , HEIGHT (NSPLIT) *. 3 
10  FORMAT(SX,/,X,I4, *  ****** ,2F11. 5, '  ***** ' , 2F11. 3 , 212 . 3) 

XT  (HEIGHT (Naplit)  .LE.  £RS OR)  SO  TO  3 


DO  302  GROSS* 1, 4 

Zand  -  CENTER (Naplit) -HEIGHT (Naplit) *dSIDE (CROSS) 

Fz2  «  F(Zand) 

XT  (ABS (Fz2 )  . £Q.  0.0D0) 

+  CALL  EXACT (Zand) 

SUM  -  (0.000,0.000) 

DO  301  LOOPS-l,Npnta 
Pll  -  Fz2 

T 82  -  F(Zand+ LOOPS *dalta*dSIDE (CROSS ) ) 

Faua  -  Fzl+Fz2 

XF  (ABS (Fz2)  .  EQ.  0.000) 

+  CALL  EXACT (Zand+LOQPS*dalta*dSIDE( CROSS) ) 

XF  (ABS (Faua)  .EQ.  0.0D0) 

+  CALL  BETWN(Zand+(LOQ?S-0.5) *dalta*dSIDE( CROSS) ,Fsua) 

SIM  -  SOM  +  Fzl/Faua 

301  CONTINUE 

INSIDE (Naplit, CROSS)  -  (2 . DO , 0 . DO) * (Npnta-(2 .DO , 0 . DO) *SUM) 

302  CONTINUE 

INSIDE (Naplit, 3)  -  INSIDE (Naplit,!) 


QUAD  -  0 

QUAD  -  QUAD+1 

QDRNTS (Naplit)  -  QUAD 

HEIGHT (Naplit+1)  -  HEIGHT (Naplit) /2 . ODO 

CENTER (Naplit+1)  -  CENTER (Naplit) +HEIGHT (Naplit+1) *CORNER (QUAD) 
WIND  -  (INTGRL(Naplit,2*QUAD-2)+INTGRL(Napllt,2*QUAD-l) 

-INSIDE (Naplit , QUAD) + INSIDE (Naplit , QUAD+1) ) 

/ (0. ODO , -3 . 283 1333 08D0) 

WINDno  -  HINT (REAL (WIND ) ) 


WRITE  (1,20)  QUAD, WIND, WINDNO, CENTER (NSPLIT+1) 
20  FORMAT (6X, 14 , 2X, 2F11. 5 , X, 14 , 2X, 2F11 . 3) 


IF  (WINDno  .EQ.  0)  GO  TO  3 

IF  (HEIGHT (Naplit+1)  .GT.  ERROR)  GO  TO  1 

3  WRITE  (1,30) 

30  FORMAT (3X,///) 

Nzaro  ■  Nzaro+1 
ZEROS (Nzaro)  -  CENTER (Naplit+1) 
RADII (Nzaro)  -  HEIGHT (Naplit+1) 
ORDER (Nzaro)  -  WINDno 

4  QUAD  -  QDRNTS (Naplit) 


5 


IF  ( QUAD . LT . 4 )  GO  TO  2 


1 
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6  Naplit  -  Naplit-1 

17  (Nap lit  .NE.  0)  SO  TO  4 

return 

END 


SUBROUTINE  BETWN ( Z , 7sua) 

COMPLEX* 18  Z , F , Fsua , Fnaw 
WRITE  (1,10)  Z 

10  FORHAT('  BETWEEN  two  points  at ' , 2E21. 14 , 1  tha  function  was  zaro1) 

Fnaw  ■  (2.000,0.000) *7(Z, a, X) 

17  (ABS(Tnav)  .  EQ.  0.0D0)  THEN 
WRITE  (1,20) 

20  FORMAT ( '  and  it  was  also  zaro  AT  THAT  POINT  l!') 

STOP  08 
ELSE 

7sua  -  Tnav 
RETURN 
ENOI? 

END 


SUBROUTINE  EXACT (Z) 

COMPLEX* 16  z 
WRITE  (1,10)  Z 

10  FORMAT ( '  WE  HIT  A  POINT  WHERE  THE  FUNCTION  IS  EXACTLY  ZERO ' , 

+  / , 2E21. 14) 

STOP  07 
END 

COMPLEX* 16  FUNCTION  F(NU) 

COMPLEX* 16  X, XL, XT, EXPO , J , Y , JP, YP, JL, JLP, YL, YLP,  JT, JTP, YT, YTP, 

+  JLP1 , JLPP1 , YLP1 , YLPP1 , JTP1 , JTPP1 , YTP1 , YTPP1 , NU , NUP1 , 

f  Oil , 012 , 013 , D2 1 , 02  2 , 02  3 , 03  2 , D3  3 , JP1 , JPP1 , YP1 , YPP1 

COMPLEX *8  XC , XLC , XTC , JC , JPC , YC , YPC , JLC , JLPC , YLC , YLPC , JTC , JTPC , 

+  YTC , YTPC , JLP1C , JLPP1C , YLP1C , YLPP1C , JTP1C , JTPP1C , YTP1C , 

+  YTPP1C , NUC , NUP1C , JP1C , JPP1C , YP1C , YPP1C 


RSAL*4  SJ , SY , SJL, SYL, SJLP1 , SYLP1 , S JT , SYT , SJTP1 , SYTP1 , SJP1 , SYP1 
REAL* 8  DENO , 0EN1 , SPEEDO , SPEEOL , S PEEDT , PI 
COMMON  X 


13 


Jun  25  08:25  1985  anu.f  Pag*  9 


EXPO*  (0 . 000 , 1.000) 

PI-3.141592554 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  Sat  tha  matarial  paraaatars  a t  tha  sphara  hara.  C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DENO—l . 0 
DEMI-13 . 800 
SPEEDO— 1.4780 
SPEEDL-6 .860 
SPEEDT— 4 . 185 


NU— NU+  ( .  5D0 ,0.000) 
NUP1— NU+ ( 1 . 0D0 , 0 . 0D0 ) 
MUC— (MU) 

MUP1C— (NUP1) 

XC-(X) 

XL— X* (SPEEDO/SPEEDL) 
XLC-(XL) 

XT-X* ( (SPEEDO/SPEEDT) ) 
XTC-(XT) 


CALL  OLVER (NUC,XC,JC,JPC,SJ,YC,YPC,SY) 

CALL  OLVER (NUP1C , XC , JP1C , JPP1C , S JP1 , YP1C , YPP1C , SYP1 ) 

CALL  OLVER (MUC , XLC , JLC , JLPC , S  JL , YLC , YLPC , SYL) 

CALL  OLVER (MUP1C , XLC , JLP1C , JLPP1C , SJLP1 , YLP1C , YLPP1C , SYLP1 ) 
CALL  OLVER (MUC , XTC , JTC , JTPC , SJT , YTC , YTPC , SYT) 

CALL  OLVER (NUP1C , XTC , JTP1C, JTPP1C , SJTP1 , YTP1C , YTPP1C , SYTP1 ) 


J-(JC*EXP(SJ) ) 

Y-(YC*EXP(SY) ) 

J?1-(JP1C*EXP(SJF1) ) 

YP1— ( YP1C*EXP ( SYF1 ) ) 

JL-  (JLC* EXP  ( SJL)  ) 

YL-(YLC*EXP(SYL) ) 

JLP1-(JLP1C*EXP(SJLP1) ) 

YLP1— (YLP1C*EXP (SYLP1) ) 

JT— (JTC*EXP(SJT) ) 

YT— (YTC*EX? (SYT) ) 

JTP1- ( JTP1C*EXP ( S JTP1 ) ) 

YTP1— (YTP1C*EXP(SYTP1) ) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  Tha  alaaanta  below  ara  givan  in  appandix  A  a f  chaptar  3 .  C 

C  Thay  ara  tha  alaaanta  of  On.  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


Dll- ( (DENO*S?EEDO**2 . ) / (DENI* SPEEDT* *2 . ) ) * (X**2 . ) 
D11-011*CSQRT (PI/ ( 2 . *X) ) * ( J+EXPO*Y) 
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NU-NU’ (.500,0.000) 

012- ( 2 . *NU* (NU+1 . ) -XT**2 . ) *JL-4 . * (NU*JL-XL* JLP1 ) 
012—012 *CSQRT ( PI/ ( 2 . *  XL) ) 

013-2 .  *NU*  (NTJ+l. )  *CSQRT  ( PI/  ( 2 .  *XT) ) 

013— 013* ( (NO-1. ) * JT-XT* JTP1 ) 

D21-CSQRT (PI/ (2 . *X)  ) 

021-021* (X* ( JP1+EXP0*VP1) -NU* ( J+EXPO**) ) 

022— CSQRT(PI/ (2. *XL) ) * (NU*JL-XL*JLP1) 

D23— NO* (NU+1. ) *CSQRT(PI/ (2 . *XT) ) *JT 

032—2. *CSQRT ( PI/ ( 2 . *XL) ) *((1.-NU) *JL+XL*JLP1) 

033-2 . *C3QRT ( PI/ ( 2 . *XT) ) * (NU*JT-XT*JTP1) 
D33-D33+CSQRT(FI/(2.*XT)  )  *(XT**2 . -2 . *NU* (NU+1) +2 . )  *JT 

F-D11*D22*D33-011*032*023-021*D12*D33+D21*D32*D13 


RETURN 

END 


SL0PE3 


This  program  generates  a  polynomial  curve  fit.  It  fits  input  data  to 
a  polynomial  of  up  to  and  including  order  5.  The  program  was  obtained  from 
the  Shock  Dynamics  Laboratory  of  Washington  State  University.  It  is  based  in 
part  on  subroutine  LSQF1T  obtained  from  the  Stanford  Accelerator  Center  by 
P.  Mars ton. 

The  program  was  used  to  fit  discrete  values  of  the  phase  velocity, 
group  velocity,  and  attenuation  of  the  Rayleigh  and  slowest  whispering  gallery 
waves  on  tungsten  carbide  spheres.  Polynomial  expressions  were  obtained  for 
these  parameters  as  a  function  of  the  ka  of  the  sphere  (cf.  Chapters  2  and  3). 
It  was  also  used  to  curve  fit  discrete  data  for  the  coupling  parameter  G 

WG1 

discussed  in  Chapcer  3  and  later  in  the  appendix  for  the  program  WS. 

The  program  is  not  listed  hare  since  polynomial'  fit  programs  are 
common.  If  the  reader  wishes  a  copy  of  the  program  ha  should  request  a  copy 
of  Internal  report  SDL  78-02  from  the  Shock  Dynamics  Laboratory  of  Washington 
State  University. 
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WS 

This  program  calculates  the  SWT  approximation  for  the  magnitude  of  the 
backacattering  form  function  of  a  tungsten  carbide  sphere .  It  outputs  the  ka 
value  and  the  form  function  magnitude  for  that  ka  in  the  range  10  £  ka  £  80. 
As  it  is  given  here  the  program  results  can  be  plotted  to  give  f^  of 
Chapter  4.  This  again  offers  a  check  that  the  program  has  been  entered 
properly.  By  changing  the  l  do  loop  limits  to  1-1,1  one  can  obtain  f 

SR 

of  Chapter  4.  By  adding  write  statements  one  can  generate  numbers  for 

plotting  the  coupling  parameters  G^  and  the  backacattering  amplitudes  used 

in  Fig.  8  of  Chapter  3  or  Fig.  2  of  Chapter  4. 

A  comment  is  important  on  the  way  the  program  calculates  the  coupling 

parameters  GR  and  GWG1  of  Chapter  4.  The  GR  is  calculated  using  the 

curve  fit  results  for  the  phase  velocity  and  attenuation  of  the  Rayleigh  wave 

resulting  from  using  the  program  SL0PE3.  Thi3  wa3  adequate  for  the  Rayleigh 

wave  but  when  done  for  the  whispering  gallery  wave  WG1  it  was  inaccurate. 

This  was  because  chose  results  were  much  more  sensitive  to  small  errors  in  the 

values  of  the  phase  velocity.  Therefore,  for  the  whispering  gallery  waves  WS 

was  first  altered  so  that  instead  of  a  do  loop  over  i  and  l  one  Just 

entered  results  for  the  ka  of  the  sphere  and  the  SWT  pole  location  from 

program  ZNU  and  outputed  an  "exact"  value  for  G.,_. .  This  was  done  for  23 

WG1 

values  of  ka  then  a  complex  curve  fit  for  G„  using  the  program  SI  OPE 3  was 
generated.  This  curve-fit  result  is  used  in  WS.  In  retrospect  the  method  for 
obtaining  GyG^  is  more  efficient  from  the  standpoint  of  computer  run  time 
then  the  current  method  for  calculating  G_  within  WS.  From  the  standpoint 
of  accuracy,  however,  GR  was  calculated  both  ways  and  was  found  to  give  the 
same  results  £.0  the  program  was  not  altered. 


The  subroutine  BD  of  this  program  calls  tha  subroutine  OLVER  given 
next  in  this  appendix.  To  do  other  isotropic  spheres  one  must  use  ZNU  and 
SLOPES  to  generate  the  appropriate  curve  fits  for  phase  velocity, 
attenuation,  and  coupling  parameters  and  alter  WS  accordingly. 
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cecccceeccccec 
c  e 
c  ws  c 
c  c 
cceccccccceccc 

ccccccc.cccceccccccccccccccccccccccccacccccccccccccccccccccccccccccccc 


c  c 
c  This  program  calculate*  an  approxiaat*  form  function  for  a  c 
c  tungsten  carbide  sph«ra  using  sommerfeld-watson  thaory.lt  c 
c  appropriataly  includas  the  spacular  reflection, Rayleigh  wave,  c 
c  and  slowest  whispering  gallary  wav*. By  including  further  writ*  c 
c  statements  ona  can  gat  th*  coupling  coefficients  for  the  surface  c 
c  waves  as  wall  as  th*  waves  amplitudes  each  time  they  circum-  c 
c  navigate  th*  sphere. On*  thing  to  note  is  that  th*  coupling  c 
c  coefficient  is  calcultad  differently  for  th*  Rayleigh  and  c 
c  whispering  gallery  wave. The  Rayleigh  calculation  uses  curve  fits  c 
c  of  th*  wave  speed  and  attenuation  to  find  the  coupling  coeffi-  c 
c  cient  .This  procaedura  was  not  accurate  enough  for  the  whisper-  c 
c  ing  gallery  wav*  so  the  exact  result  for  the  coupling  coeffi-  c 
a  cient  was  calculated  at  several  ka  values  and  than  these  results  c 
c  were  curve  fit  to  give  continuous  values.  c 
c  The  subroutine  BO  of  this  program  uses  th*  subroutine  OLVER  c 
a  which  is  given  elsewhere  in  this  appendix.  c 
c  c 


cccccccccsccccccccccccccccccccccccccccccecccccccccccccccccccccccccccc 

dimension  fabc(lOOO) , kan(iooo) ,cupla(1000,2) , speca(iooo) 
dimension  fabst(lOOO) 

reel *4  vel , svel , lvel , theta , beta , theta , thetp , dans , deni , dum , 

4  ray 1 , ray 2 , ray 3 , ray 4 , kan , f aba , ka , rvel , aspsc , betp , betm 

4  , f abst 

complexes  nu , dminua , dplus , phasal , phas*2 , x , xs , xl , num , nup , 

4  ddplus,  dpiusm,  dpiusp,. cn, in, couple, decay l,d*cay2 , 

4  spec , cthatm , cthstp , cthsta , cf abs , cupl , cupla , specs 

pi-3.141592654 

lud-1 

cccccccccccccccccccccocccccccccceccccccccccccccccccccccccccccccccccccccc 
c  Set  material  parameters  for  the  sphere. These  parameters  are  same  as 
a  those  used  by  Saunaurd  in  JASA  Jan.  1983. 
cceccccccccccccceccccccccccccccccccccccccccccccccccccccccccccccccecccccc 

numb- 10 00 
vel-1 . 4760 
svel-4 . 135 
lvel-6. 860 
dens-13 . 3 
deni- 1. 0 
a- 12. 7 

ccccccccccccccccccoccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  Th*  i  do  loop  calculates  th*  required  Sommerfeld-Watson  theory  (SWT)  c 
c  results  for  a  particular  surface  wave  for  a  range  of  ka.Th*  1  do  c 

c  loop  selects  th*  particular  surface  wav*;l-l  for  th*  Rayleigh  wave,  c 


l>  o 
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c  1*2  for  tha  whiaparing  gallary  wava.  c 

cccccccccccccccccccceccccccccccccccceeeeccccccccecccecccccccceccceccccec 

do  301  1-1,2 
do  100  i-l,nuab 

cccccccccecccccecceccceccccceccecccececcccccccececceecccecccccccccccccec 
e  This  part  o f  tha  loop  calculatas  surfaca  wava  spaads  and  attan- 
c  uations.Tha  calculations  uaa  rasulta  of  a  program  which  curva 
c  fits  data  to  a  polynomial  of  up  to  (ka)**5. 
ccccccecceccecccccecccccccccccccccccccccacccccccceccecccccccceeccccccccc 

ka-10.+(  (i-1)  *70 .  )/numh 
kiua(i)  -ka 
x-cmplx(ka,0. 0) 
if  (l.atr.-S)  »<•>  to  211 
bata-.  0;:^2134if  .2X0487a-2*ka 
rval-3  *  4C413-. 24*972a-l*ka+.331109a-3*fc:i<**2 
&  ".157333a-5*}ca**3 

go  to  2.U 

211  bata-.  2<l«;t2e«*  1-.  306817a-2*ka+.  123060a-3*ka**2 

6  •*,  i««263a-5*ka**3+.752335a-8*ka**4 

423  if  {Xm  .It.  45.2)  go  to  417 

rv«l-5. 480848912 *axp (-. 00437092*ka) 
go  to  213 

417  if  (ka  .It.  29.0)  go  to  418 

rval-3. 498 195789*axp(-.00428319*ka) 
go  to  213 

41 A  rval-5.51968062*axp(-.00447460*ka) 

go  to  213 

21-  continua 

ccccccccccecccccccccccccccccccccccccccccccccecccccccccccccccccccccccccc 
c  This  saction  calculatas  tha  coupling  eoaffciants  for  tha  surfaca  c 
c  wava. Kota  again  that  tha  coupling  coafficiant  is  calculatad  diff-  c 
z  arantly  for  tha  two  surfaca  vavaa.In  ratrospact  it  is  probably  c 
:  sora  afficiant  to  calculata  axact  rasults  for  tha  coupling  coaff-  c 
c  iciants  at  savaral  ka  valuas  and  than  curva  fit  thosa  rasults. This  c 
c  aaauras  accuracy  and  dacraasaa  computar  run  tima.Tha  Raylaigh  c 

c  calculation  was  dona  both  ways  and  tha  sama  accuracy  was  found  c 
c  but  this  is  not  ganarally  found  to  ba  tha  casa. Showing  tha  two  c 
c  diffarant  mathods  hara  allows  a  comparison  of  tachniquas . In  any  c 
c  casa  ona  still  naads  tha  curva  fits  for  tha  surfaca  wava  spaads  c 
c  and  attanuations  givan  abova  for  furthar  calculation.  c 

ccccccccccecccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

dum-ka*val/sval 
xs-caplx ( dum ,0.0) 
dum-ka*val/lval 
xl-cmplx (dum , 0 . 0 ) 

div-1000. 

that a-a  s in ( 1 . / rva 1 ) 

thatm-thata-thata/div 

thatp-thata+thata/div 


o  o  o 
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if  (1  .  aq.  3)  go  to  707 
dua~!«a*«  in  (theta) >.5 
nu~cmplx(dua,bata) 
dua-ka*sin  ( theta)  - .  5 
nua-eaplx(dua,bata) 
dua-ka*ain(thatp) -.5 
nup-caplx ( duo , bata ) 

call  bd ( x , xa , xl , num , dana , danl , dainua , dpluam ) 
call  bd ( x , xa , xl , nup , dana , danl , dainua , dpluap ) 
call  bd ( x , xa , xl , nu , dana , danl , dainua , dplua ) 


ddplus-(dplusp-dpluam) / (nup-nua) 

707  continue 

if  (l.aq.2)  go  to  703 

coupla~(4. *pi*caplx(0.0, 1. 0) *(nu+.5) *dainua) / (ddplus*ka) 
go  to  704 

703  dual—  .133352+.  «44851e-2*)«a-.  353597a-4*ka**2 
dua2-.021110-.  84413  «e-3*lca+.S72282a-5*)ca**2 
coupla«caplx ( dual , dua2 ) 

ccccccccccccccccccccceeccececccccccecccececceccccccccccccccccccccccccccc 
c  Tha  axponanial  decay  and  propagation  pbaaa  shifts  ara  calculated  c 
c  hero  using  tha  appropriate  curve  fits.  c 

cccccecccccccccccccccccccccccccccccccccccccccccceccccccccccccccccccccccc 

704  dua-axp(  (2.  *pi-2 .  *thata)  * (-J.. )  *bata) 
decay 1-caplx ( dua ,0.0) 

dum*»exp(  (-1. )  *2.  *pi*beta) 
decay2~caplx ( dun ,0.0) 

al-coe(2  .  *ka*  (1.  -cos (theta)  )  +ka*sin (theta)  * 
t  (2. *pi-2. *thata) -pi/2.0) 

a2*sin(2 .  *ka*(l.  -coa  (theta) )  +ka*s  in  (theta)  * 
i  (2.*pi-2. *theta) -pi/2.0) 

phasel»caplx ( al , a2 ) 

al~coe(2 . *pi*ka*sin (theta) -pi) 
a2«*in(2 .  *pi*lca*sin  (theta)  -pi) 
pfaaaa2-caplx  (al ,  a2 ) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  The  specular  reflection  is  calculated  at  this  point.  c 

ccscccccsccccccceccccccccccccccccccccceceecccccecccccccccccccccccccccc 

mi-caplx (-.5,0.0) 

call  ref 1 ( x , xa , xl , nu , dana , danl , dainua , dplua ) 

apec~dair.ua/dplus 

specs (i) -spec 

asp «c- cabs ( spec ) 

ccccccececcccceccccccececccccccccccccccccecccccccccccccccccccccccceccc 
c  The  SWT  terms  for  tha  surface  wave  contribution*  to  tha  form  c 

c  function  are  calculated  using  the  results  above.  c 


143 


Jun  10  11:35  1985  wsef.f  Pag*  4 


ccecccccccccccccccccccccceccccceccccccccccccccccccccececcccccccecccccc 

cupl-phasel*dacayl*couple/ (l.-decay2*phas*2) 
cupla(i,l)-cupl 


100  continue 

301  continu* 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

a  Th*  approximate  form  function  is  calculated. Sine*  th*  SWT 
c  results  for  th*  whispering  gallary  wav*  ar*  not  accurat*  below 
c  k*«2l  its  contribution  is  not  included  for  10<ka<21. 
ccecccccccccccecceccccccccecccecccceccccccccccccccccccccccccccccccccc 

do  302  i-l,nuab 

if  (kaa(i)  .It.  21.0)  go  to  303 

fabst (i) -cabs ( specs (i) +eupla ( i, l)+cupla(i, 2) ) 

go  to  304 

303  fabst (i) -cabs ( specs ( i ) +cupla( i, 1) } 

304  writa(lud,315)  Jeaa(i) ,  fabst (i) 

315  format(3x,3fl2.5) 

302  continu* 

end 


non 
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ccccceccccccccccccccccccccccccccccceccccccccecccceeccccccccceeceecccccccc 
e  The  subroutine  bd  calculate*  tha  0-  and  D+  of  tha  SWT.lt  usaa  the 
c  subroutine  olvar  which  can  calculate  baaaal  function*  of  complex 
c  order  and  arguaaantt. 
ccccccccGcccoceccccccccccccccccecccccccccccccceeeceeeeeeeececccccccccccce 

subroutine  bd(x,xa,xl,nu,d*n»,denl,da,dp) 

co*plex*8  x ,  x* ,  xl ,  da,  dp ,  rl ,  r2 ,  rs ,  rl ,  nu ,  *2 , aa2 , al2 , z , 

6  ai,aip,cdua, j , jp,y,yp 

real *4  dana,danl,ka*,r,a,kal,dua,sj  ,sy 

e2-(nu+.5) 

call  olvar ( a2 ,  x ,  j ,  j  p ,  a  j ,  y ,  yp ,  *y ) 

jp-jp*axp(aj) 

j-j*exp<sj) 

yp-yp*«xp(sy) 

y-y*axp(sy) 

rl-( jp+cmplx(0.0, 1.0) *yp)/ ( j+cmplx(0. 0, 1. 0) *y)  “(2.  *x) **(-!.) 

r2-(  jp-caplx(0.0,  l.o)  *yp)/ ( j-caplx(0. 0, 1. 0)  *y)  -(2. *x) **(-1. ) 

call  olver(e2,xa, j, jp,sj ,y,yp,sy) 

j-j*axp(sj) 

jp-jp*axp(aj) 

rs-jp/ j-<2.*xa) **(-i) 

call  olvar(*2,xl, j,jp,sj,y,yp,sy) 

j-j*axp(aj) 

jp-jp*axp(aj) 

rl-jp/j-(2.*xl)**(-l) 


da- (danl/dana ) *xs*xs* ( (xl*rl) *  (2. *xs*rs+(xs*xs-2 . ‘nu* 
s  (nu+1. )  +2. )  )-(2.  *(l.-xl*rl)  *nu*  (nu-t-l. )  ) ) 

dn-da+x*r2*(  ( (2 .  *nu* (nu+l. )  -ocs*xs)  -4.  *xl*rl)  *  (2 .  *xs*rs+- 
5  (xs*xa-2. *nu* (nu+1. }«*2. ) ) -(2. * (l.-xl*rl) ) *(2 . *nu* 

i  (nu+l. ) *(xa*ra-l. ) ) ) 

dp-(danl/dana) *xs*xs*( (xl*rl) *(2. *xs*rs+(xs*xs-2 . *nu* 
i  (nu+1. ) +2. ) ) -(2. *(l.-xl*rl) *nu*(nu+l. ) ) ) 

dp-dp+x*rl*( ( (2.  *nu*( nu-t-l. ) -xa*xa) -4. *xl*rl) *(2 . *xa*rs+ 
i  (xa*xa-2 . *nu*(nu+l. )+2 . ) ) -(2 . *(l.-xl*rl) )  *(2. *nu* 

s  (nu+1. ) * (xs*rs-l. ) ) ) 

raturn 

and 

cccccccccccccccccccccceecceeeeecccccccccccccccccceccccccecccccccccccccc 
c  Tha  aubroutina  rafl  calculate*  the  reflection  coefficient  from  tha  c 
c  front  of  tha  sphera.Xt  doa*  not  uaa  olvar  sine*  olvar  is  not  accur-  c 
c  ata  for  order  of  lea*  than  2  and  for  backscattsring  va  need  answers  c 
c  for  order  o.This  aubroutina  uses  large  arguement  expansion*  for  c 

c  tha  besael  and  hankel  function*  and  naglacts  terms  of  0(ka+*>2) .  c 
cecccceccccccccccccecceccccccccceeccccccceccccccccccccccccccccccccccccc 

subroutine  raf 1 ( x , xs , xl , nu , dans , deni , dm , dp ) 


nan 
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coaplax*i  x, xa, xl,nu, da,dp, a2 , zl, z2a, z21,b22,a22, 
*  a2,aa2,al2, j, jp,y,yp 

raal*4  dan* , danl , a j , ay 


a2-(nu+.5)  *(nu-f.5)/(x*x) 
aaa- (nu+ .  3)  *  (nu+ .  9 )  /  ( xs*xa ) 
ala«(nu+.3) *(nu+.3)/(xl*xl) 

a2-x*(caplx(0. , -1. ) *aqrt(l. -a2) -(2 . -«2)/ (2 . *x* 
6  (1. 0-«2) ) ) 

zl»x*<caplx(0. , 1.) *sqr*(l. -«2) -(2. -a2)/ (2. *x* 

4  (1.0-«2))) 

a2a«xa* ( caplx ( 0 . ,-l. ) *aqrt(l. -«a2) -(2.-aa2)/ 

4  (2.*xa*(l.-«a2) ) } 

z21-xl*(caplx(0. ,-l. ) *aqrt ( 1 . -«12 ) - ( 2 . -al2 ) / 

4  (2.*xl*(l.-«12) )  ) 


*22«a21*(l./4.+l./ (2. «xa*xa) +z2s/ (2. *xs*x a) ) 

4  -nu* (nu+1. )/ (2. *xa*xs) 

b22«(-l./4.-l./ (2. *xs*xs)+(nu+l. )  *nu/  (xs*xs) )  + 

4  (2.*nu+(nu+l . )/ (xs*xa*xs*xs) ) - 

4  nu+nu* ( nu+1 . ) * ( nu+1 . ) / (xa*xa*xs*xs ) 

b22-b22- 

4  a2a/(2.*xa*xs)+(nu+(nu+l.  )/(xs+xa*xs*xa) 

4  -1-/ (xa*xa) -2./ (xs»xs*xs*xs) ) *z21 

4  +(nu*(nu+l. )/ (xa*xs*xs*xs) -2 ./ (xs*xa*xs*xa) ) 

4  *z2a*z21 

d»-z2+b22+(danl/dans) *a22 
dp-zl*b22+(danl/dana) *a22 


ratum 

and 


i** 


146 


OLVER 

This  program  calculates  Bassal  functions  of  tha  first  and  second  kind 
J^(x)  and  Ny(x)  as  veil  as  their  derivatives  with  respect  to  x  for 
complex  order  v  and  argument  x.  The  program  was  obtained  from  J.  A. 

Cochran  of  the  Department  of  Mathematics  at  Washington  State  University*  The 
program  Is  asserted  to  have  an  error  of  less  than  (0.00001  +  10.00001) 
wherever  |v|  ^  3  and  x/v  >_  10*"^.  There  are  some  comments  In  Chapters  3  and 
4  on  the  program  as  well  as  a  reference  to  earlier  work  on  the  numerical 
computation  of  Bessel  functions  by  Cochran  at  al.  (Ref.  22  of  Chapter  2). 

Much  more  extensive  comments  can  be  found  In  the  Master  degree  paper  (Ref.  21 
of  Chapter  2)  by  Brad  Brim. 

OLVER  makes  a  decision  of  whether  to  scale  the  values  of  J^(x)  and 
J^(x)  by  exp(jscale)  and  tha  values  of  N^(x)  and  }T(x)  by  exp(yscale) 
for  convenient  output.  In  the  programs  ZNU,  ZKA,  and  WS  tha  scale  factors 
are  included  In  all  calculations.  If  one  desires  to  use  OLVER  for  ocher  work 
the  programs  should  be  checked  to  see  how  to  include  the  scale  factors. 


J'  *• 


vr -  vr< ;  .^.'t 
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cccccccccecccccece 
c  e 
c  o&vra  c 
c  e 
cccccccccecccccccc 


subroutine  olver ( argnu , argw , j , j prim* , j  seal • , y , yprime , yscal • ) 


C* **************************************************** ****************** 

c*  larg*  order  routine  to  generate  bessel  functions  of  first  and  * 

c*  second  kinds  and  their  derivatives  for  complex  values  of  order  * 

c*  argnu  and  argument  argw. written  by  e.j. murphy  1970/71, revised  * 

c*  by  j.a.  cochran  1/72,  revised  by  e.r.  kamgnia  7/1977.  * 

c* ********************* ********************** *************************** 


common  /info/  uc,vc,ac,bc,a,b,c,d,cosl20,cos30,third 

complex* 8  argnu, argw, j , jprime,y, yprime 

complex* 16  w, nu, nuto!3 , nuto23 , nuto43 , nuto2 , nuto4 , nu23xz ; 

>  a(«,2) ,b(<5,2)  ,c(«,2) , d (6 , 2 ) , ar ( 2 ) ,br(2) ,er(2) ,dr(2) , 

>■  u  ( 5 )  ,  v  ( 5 ) ,  a  sum ,  bsum ,  esum ,  dsum , 

i-  z ,  zeta ,  z*tal2 ,  ztam3 2 ,  phiof z ,  tauinv,  tl ,  t2 ,  t3  ,  t4 , 

■  ai , aip , bi , bip , ail , ailp , ai2 , ai2p , 

y  coel20,cos30, 

i-  x,t,terml,term2 

r*al*8  uc(14) ,vc(14) ,ac(5) ,bc(5) , theta, third, si, s2,s3, 
y  realnu , imagnu , magnu 

real*4  j scale, yscal* 


c*  initialization  * 

c  *********************************************************************** 

j  scale  ■  0.0 
yscale  -  0.0 
w  -  argw 

nu  -  argnu 
realnu  ■  real(nu) 
imagnu  -  imag(nu) 

z  -  w/nu 

theta  ■  third*datan2 (imagnu, realnu) 
magnu  “  abs (nu) 

nutol3  ■  (magnu**third) *caplx(dcos (theta) ,dsin (theta) ) 

nuto23  -  nutol3*nutol3 

nuto43  ■  nuto23*nuto23 

nuto2  ■  nu*nu 

nuto4  ■  r.uto2*nuto2 

nterms  -  l 

if  (magnu  .It.  25.0d0)  nterms  ■  2 


v. 


OOQ 
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C** A************************************************************ ******** 

c*  coaputa  z*t*  * 

2*********************************************************************** 

call  azata ( z , z*tal2 , phiof z , tauinv) 

zata  -  zatal2*zatal2 


q  ■*  ★  h  *  *  *  +  ★  **  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  ilr  *  *  *  *  A  *  *  *  ★  *  *  *  ★  *  *  *  *  *  *  *  *  * •*  *  ★  ★  irk  *  *  *  ★  *  ★  ★  ★  ★  *  ★  ★  ★ 

a*  tha- magnituda  of  z  datarainaa  how  v*  coaput a  ar,br,cr,dr  * 

c *********************************************************************** 


if  (abo(z-l.dO)  .gt.  O.ldO)  than 


****************************************************************** 

*  eomputa  ar,br,cr,dr  in  tha  normal  way  ***** 

****************************************************************** 

tl  -  (l.dO, 0.d0)/tauinv 

ta  •  ti*ti 

t  *■  t2 

u(l)  -  (I.d0,0.d0) 

v(l)  -  (I.d0,0.d0) 

u(2)  ■  tl*(uc(l) -uc(2) *t) 

v(2)  «-tl* (vc ( 1) -vc (2 ) *t ) 

u(3)  •  t2*(uc(3)-f*(uc(4)-t*uc(5) ) ) 

v<3)  — t2*(vc(3) -t*(vc(4) -t*vc(5) ) ) 

if  (ntarma  . aq.  2)  than 
t3  -  tl*t2 

ta  -  ta*t2 

U(4)  «  t3* (UC(6) -t* (UC(7) -t* (UC(8) -t*UC (9) ) ) ) 
v(4)  — 13* (vc (6) -t* (vc (7) -t* (vc( 8) -t*vc (9) ) ) ) 
u( 3)  -  t4*(uc(10)-t*(uc(ll)-t*(uc(12)-t*(uc(13) -t*UC(14) ) ) ) ) 
v(5)  —t4* (vc(10) -t*(vc(ll) -t* (vc(12) -t* (vc(13 ) -t*vc ( 14) ) )  ) ) 
and  if 

ztaa32  *  ( 1. OdO , 0 . OdO) / (zata*zatal2) 
do  3  i«l,ntarraa 

ar(i)  •  (0. OdO, o. OdO) 
br(i)  «  (0. OdO, 0, OdO) 
cr(i)  -  (o.odo.o.odo) 
dr(i)  -  (0. OdO, 0. OdO) 
taral  -  (1. OdO, 0. OdO) 
tarma  u  (1. OdO, 0. 0d0)/zatal2 
do  6  k-l,2*i+l 

1  -  2*  (i+1)  -k 

ar(i)  -  ar(i)+bc(k)  *u(l)  *taral 
dr(i)  ■  dr(i)+ac(k) *v(l) *taral 
taral  ■  taral*ztam32 
6  continua 

do  7  k-l,2*i 

1  -  2*i+l-k 

br(i)  ■  br (i) -ac(k) *u(l) *tarn2 

cr(i)  •  cr(i) -bc(k) *v(l) *tara2 
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tarsi  *  tarn2*ztaa32 
7  centlnu* 

cr(i)  -  cr(i)*zata 
3  centlnu* 


•Is* 

^  ***************************** ************************* ************ 

c  *  cosputa  ar,br,cr,4  dr  by  powar  sarlaa  ***** 

c  ****************************************************************** 

do  11  i-l,ntaraa 
ar(i)-a(6,i) 
br(i)-b(«,i) 
cr(i)-e(«,i) 
dr(i)-d(«,i) 
do  10  1-5, 1,-1 

ar(i)-ar(i) *zata+a(l,i) 
br(i)«br(i) *zata+b(i,i) 
er(i)-cr(i) *zata+c(l,i) 
dr(i)«dr(i) *zata+d(l,i) 

10  contlnua 

11  contlnua 

and  If 


c *********************************************************************** 
c*  computa  airy  functions  * 

e*********************************************************************** 

nu23xz  -  nuto23*zata 
x  -  nu23xz 

call  airy(x,ai,aip,s3) 

x  -  coal20*nu23xz 
call  airy (x, ail, ailp,sl) 

x  >  conjg(cosl20) *nu23xz 
call  airy(x,ai2,ai2p,s2) 


c************< 
c*  scaling 
c ************' 


al  -  al-a2 
jacala  -  a3 
yacala  -  *2 

if  (dabs (si)  ,gt.  170. OdO)  than 

if  (si  .gt.  0 . OdO)  than 
/seals  -  yscala+sl 
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ai2  -  (O.OdO, O.OdO) 

al2p  -  (O.OdO, O.OdO) 

else 

all  *  (O.OdO, O.OdO) 

ailp  -  (O.OdO, O.OdO) 

•ndif 


else 


ail  ■  ail*daxp(sl) 
allp  -  ailp*daxp  (*1) 

•ndif 

bi  -  cos30*ail+conjg(cos30) *ai2 
blp  -  -(conjg(coa30) *ailp+coa30*ai2p) 


c*  compute  a,b,c,t  d  sums  * 

e*******************************************************^******wt--‘'-***.'.  » 

aaua.  -  (1. OdO, 0. OdO)  +  ar(l)/nuto2 
baua  -  br(l) 
caua  -  cr(l) 

daua  •  (1. OdO, O.OdO)  +  dr(l)/nuto2 

15  (ntaras  . aq.  2)  than 
aaua  -  aaua+ar ( 2 ) /nuto4 
baua  -  bsua+br(2)/nuto2 
caua  ■  caua+cr ( 1) /nuto2 
daua  -  daua+dr ( 2 ) /nuto4 
•ndif 


C****** ******************************************** ******»•• 
c*  fora  j,y, jpria«,ypria«  t  than  raturn  the  results 
,:*********••*•***•*••*••****•***•*•«:*••••*•**•**•****•***•*• 


tarml  ”  bsua/nuto43 
tara2  -  phiofz/nutol3 
j  ■  tara2*(aip*taral+ai*asua) 

y  «-tara2*(bip*taral+bi*aaua) 

taral  *  (2.d0,0.d0)/nuto23/z/phiofz 
tara2  “  csua/nuto23 
jpriaa  «-taral*(ai*tarm2+aip*daua) 
ypriaa  •  taral* (bi*tara2+bip*daua) 


raturn 


151 
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subroutine  szsta( z , zsta, phiz , r) 

c  ***  double  precision  complex  calculation  of  sqrt(zeta)  ,phi(z*ta)  , 
c  ***  and  sqrt(.l-z*z)  . 

Q  *** 

complex* 16  z , zeta , phiz , r 

double  precision  pi2 ,al3 , zr, zi, zm,az, argz 

data  pi2/6 . 283 183 3 07 17938 6d0/ , al3/ . 333 333 3 33 33 33 3 3 3d0/ 

zr«z 

zi»(0.d0, -l.dO) *z 
if (zi.ne. o.do)  go  to  l 
if (zr.eq. l.dO)  go  to  6 
if (zr.gt.l.ldO)  go  to  3 

c  *** 

c  ***  z  is  not  real  or  is  real  <1.1 

c 

1  r»sqrt(l.dO-z*z) 

if  (r*al(r)  „lt.  O.OdO)  r  -  -r 
argz»datan2 ( zi , zr ) 

Q  *** 

c  ***  compute  zata**1.5 

c 

zeta-1.3d0*(log( (l.d0+r)/z) -r) 
zm-abs ( zeta) 

if (zm.le.l.d-12)  go  to  6 
c  ***  place  zeta**1.3  in  the  proper  quadrant 

Q 

zr-zeta 

zi«(0.d0, -l.dO) *zeta 
az-datan2 (zi, zr) 
if (zr.gt.O.dO)  go  to  4 
if (zi.le.o.do)  go  to  3 
if (argz. It. 0. do )  go  to  4 
az-az-pi2 
go  to  4 


>0  s"  S’ 
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3  if (argz.gt.O.dO)  go  to  4 
az-az+pi2 

Q  *** 

c  ***  computa  aqrt(zata) 

4  az-az*al3 

zata-(zm**al3) *cmplx(dcoa (az) ,dsin(az) ) 
c  ***  computa  phi.  (zata) 

G 

2  pfciz-aqrt(2.d0*zata/r) 

if  (raal(phiz)  .It.  O.OdO)  phiz  <■  -phiz 
return 

G  ha* 

o  ***  z  is  raal  and  '■1.1 

Q  *** 

5  zi-dsqrt(t,r*zr-l.dO) 
r-caplx ( o . do , -zi) 

C  *** 

u  **+  computa  zata+*1.5 

G  *** 

zm-dab»(1.3do*(datan(zi)  -zi) ) 
zata-caplx ( 0 „ do , -zm*+al3 ) 
go  to  2 

o  *•  * 

c  ***  zata-o 
c  *** 

6  phiz— ( 1 . 2599210499dQ , O.dO) 
zata- (0 .do , O.dO) 

r-(0.do,o,  -.0) 

ratum 

and 

C  **  * 

subroutina  airy(::,ai,aip,  faxp) 
dimanaion  iaxp(2) 

daub la  pracisinn  n,o(4, 6) , rapi, tamp, faxp 

compiax+16  x,z,?,g,h,l, tarml, tarm2 , z2 , y (4) , a (4) , f c (2 ) , t , ai, aip 
data  rapi/ .5641895835477563d0/,c/ 

+1513. 439537849732d0, 274 . 0268 600980441d0, -1605. 28 1761546393d0, 
+-293 . 58 9 24 102 2 3 63 OdO , 

+54. 73651753365000d0( 12. 08406751643 030d0, -58. 322383 06494943d0, 
+-13 . 194S2747124167d0 , 

+3 . 01523 7 03 552 022 8d0 , . 3318825830603329d0 , -3 . 30764827S323637d0 , 
+-.»606547161952595d0, 

+. 28 3 034 83 8 7 667 13 ldO , . 1109516996742 113d0 ,-. 3217286160943928d0 , 
+-. 132 69 220464 677 63 d0 , 

+.0350974151234368d0, . 0355259773 662552d0 0670090663580247d0 , 
+-.0470357510238066d0, 

+.0347222222222222d0, . 0694444444444444d0 , -. 0486111111111111d0 , 
+- . 097 22 222 2 22 2 22 2 2d0/ 

c  *+* 

c  *+*  initializa 

C  **  + 

faxp-O.dO 

z-x 

tamp-z+conjg(z) 


153 
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if (taap.gt.3S.dO)  go  to  2 

c  *** 

o  ***  powar  aarias 
c  *•* 

taral»( .3 530280538873 17d0, O.dO) 

tarn2-(.2588194O3792SO7d0,  O.dO) 

f»taral 

l-tarm2 

tara2-z*tarn2 

g-t*ra2 

Z2 -z*z 

h“  ( 0 .  do ,  0 .  dO ) 
n-l.do 

1  n-n+l.dO 
tarml-z2*tarnl/n 
h-h+taral 
n-n+l.do 
tarml-z *t aral/n 
f-f+tarml 
tarm2-z2,»tarm2/n 
l-l+tara2 
n-n+l.do 
tarn2-z+t*rm2/n 
g-g+tarn2 

t*ap«t arml *cen j  g ( tarml ) 

if (tamp.gt.l.d-32)  go  to  1 

ai-f-g 

aip-h-L 

raturn 

c  *** 

c  ***  phaaa  amplituda  axpanaiona 
c  *** 

2  n-l.do 

if (raal(z) .gt.O. )  go  to  3 
n— l.dO 

Zm-Z 

3  h-sqrt ( z ) 

if  (raal(h)  .It.  o.odo)  h  -  -h 
l-sqrt(h) 

if  (raal(l)  .It.  O.OdO)  1  -  -l 

f-z+h 

t-(1.5d0, 0.d0)/f 
g-t+t+n 

3c(l)-caplx(rspi,0.d0)/l 
fc(2)«caplx(rspi,0.d0) *1 

c  *+* 

c  **+  conputa  phasaa  and  aaplitudas 
c  *+* 

do  3  i— 1 , 4 
y(i)-(0.do,o.d0) 

do  4  j-1,6 
y(i)-g*(y(i)+c(i, j)) 

4  continu* 
y(i)-y(i)+l.do 

5  continua 

do  6  i-2,4,2 


Jun  10  11:36  198S  olv*r,f  Pag*  3 


y(i)-y(i)/t 

6  continue 

if (n.gt.O.dO)  go  to  7 

c  *** 

c  ***  sat  values  for  cosine  expansions 
c  *** 

y(2)-y(2)-(0.78S398163397448d0,0.d0) 
y  ( 4 )  **y  ( 4 )  +  ( 3  . 92 69908 16987 2 4 ldO  ,  0 .  dO ) 
n-l.dO 
go  to  9 

c  *** 

c  ***  sat  values  for  exponential  expansions 

G  *** 

7  continue 

do  8  i-2,4,2 
y(i)— y(i) 

3  continue 
n-2 . do 

c  ***  fora  ai  and  aip 
c  *** 

9  continue 

do  10  i-2,4,2 

j-i-1 

k-i/2 

call  func(n,y(i) , f , iexp(k) ) 
a(Jc)  -y  ( j )  *f*fc(k)  /n 
la  continue 

if (iaxp(l) ,eq.0.and.iexp(2) .eq.O)  go  to  11 
f*xp-dble( float (i*xp(l) ) ) 

a(2) -a (2) *dexp (dbl*( float ( iaxp(2) -laxp ( 1) ) ) ) 
11  ai-e(l) 
aip-e ( 2 ) 

if (n.n*.2.d0)  return 

aip— aip 

return 

end 

c  *** 

subroutine  func (n, z , v, i) 
complex* 16  z,w 

double  precision  n,x,y,r,u,v,t*xtra 
taxtra— .9210340371976183d! 
if (n.eq.l.dO)  go  to  2 

Q 

c  ***  calculate  exponential 
c 

i«int(real(z) ) 

v-exp(z-cmplx(dble( float (i) ) , O.dO) ) 
return 

Q  ##★ 

e  ***  calculate  cosine 
c  *** 

2  x-Z 

y— (O.dO, -l, do ) *z 
r-dabs (y) 
i-idint (r) 
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v»d*xp(r-dbl*( float (i) ) )  *.SdO 
if (r.gt.taxtra)  go  Co  4 
r»d*xp(-2.d0*r) 
go  to  6 

4  r-O.dO 

6  u-v*(l.do+r) *dcoa(x) 
v»v*(r-l.dO)  *d*in(x) 
if (y. It. 0. do)  v— v 
w»caplx(u, v) 
return 


block  data 

comaon  /info/  uc,vc,ac,be,a,b,e,d,co*120,co*30,third 

conpl«x*16  w, nu, nutol3 , nuto23 , nuto43 , nuto2 , nuto4 , nu23xz , 

♦  a(6,2) ,b(6,2) , c(6, 2) , d ( 6 , 2 ) ,ar(2) ,br(2) ,er(2) ,dr(2) , 

+  u(3) ,v(3) ,aaua,bsun, csum,dsun, 

+  » ,  zata , z*tal2 , ztan3  2 , phiof z , tauinv , tl , t2 , t3 , t4 , 

+  al , aip , bi , bip , ail , ailp , ai2 , ai2p , 

>  00*120,00*30, 

+  x,t,t*rnl,t*ra2 

r*al*8  uc(14) ,vc(14) ,ac(5) ,bc(5) , th*ta, third, si, s2,s3, 

+  raalnu , imagnu , magnu 

data  uc/ 

+  . 12 3 00 000 00 00 000 Od  0,  .2083333333333333d  0,  . 0703125000000000d  0, 

4>  .4010416666666667d  0,  . 3342013888888889d  0,  . 0732421875000000d  0, 

+  . 3912 1093 73000000d  0,  . 1846462673 611111d+l ,  . 1023812396430617d+l , 
+  • 1121520996093750d  0,  . 23640S6914062500d+l ,  . 8789123335156243d+l , 

+  .  11207002 616223 00d+2  ,  .  4669384423426247d+l/ 

data  vc/ 

+  . 3 7 30 000 00000 00 00d  0,  . 2916666666666667d  0,  . 1171873000000000d  0, 

+  . 51362500000Q0000d  0,  . 3949632777777778d  0,  . 1023390625000000d  0, 

+  . 10892373 12500000d+l ,  . 2130533854166667d+l,  . 114649643 1327160d+l , 
+  . 14419 355 664062 30d  0,  . 2793920893437500d+l ,  . 9961006673177080d+l , 

+  . 123 8 6687 102 14 12 0d+ 2 ,  . 3073635242854613d+l/ 

data  ac/ 

+  .1000000000000000d+l,  .  1041666666666667d  0,  .  3335034722222222d-l , 

+  . 12822637 4556327 Id  0,  . 2918490264641404d  0/ 

data  be/ 

+  . 1000000000000000d+l,  145333333333333'*'*  0 987413 1944444444d-l , 

+-.1433120539158950d  0, 31722720267841.'.  V 

data  a/ 

+  (-. 444 4 44 4 64 3 9 638 60d-0 2 ,  . 1420578595930520d-13)  , 

+ (- . 1463 7 077 8 9 13 904 9d“02 ,  . 6340517965151910d-15) , 

+(  . 7064169774714094d“03 ,  . 1150347373 14 0440d-13 ) , 

+(  . 6729874534333S3Sd"03 ,  . 4301377701085880d-13 ) , 

+(  . 134002 42 88 054 082d-03 , 12 3 14 13 8 2 3 62 1110d-ll ) , 
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+  (  — ,  5766054 14 88 9 113 2d-04 , 
*(  . 6937359165991589d-03, 
+  (  .  3686593402959254d-03  , 
+(--2698466000992447d-03, 
+(-.3512938924677978d-03, 
+  (-. 102 57 87 2 2 62 3 05 oOd-O 3 , 
+(  ,3103622516889943d-04, 

data  b/ 

+(  . 17 9988722084405 Id -01, 
+  (  .888888903661063513-02, 
+(  . 1625687323118995d-02 , 
+(-.3642847477297854d-03, 
+  (-. 3 02 0603 9 07 8 3 827 6d-03 , 
+(-.58443«7298607676d-04, 
+  (-. 1492829626703820d-02, 
+  (-. 1394063 233151223d-02, 
+  (-. 3 8 2096 14 987 088 59d-0 3 , 
+(  . 16909714262S3441d-03 , 
+(  . 1709979970784747d-03 , 
+  (  . 4163 922 9 4 68 3 448 5d-04 , 

data  c/ 

+(  . 158740 10558 16959d  00, 
+(  ,2519842152827235d-01, 
+(-.3301586923087801d-02, 
+(-. 2356591754094817d-02 , 
+  (-.  8  674  2964.8  6680  654d-04  , 
+(  ,2775244443466439d-03, 
+(-. 21692 1922963 63 lld-02 , 
+( -.3443422 45622 139 3d-03 , 
+(  . 7803 7825 3 22 3 2 62 5d-03 , 
+(  . 3 8 13 588 013 2 150 2 Sd-03 , 
+  (  - . 93 954262994 572 43d-04 , 
+  (-. 15039 12 08622 063 2d-03 , 

data  d/ 

+(  . 7 3 01587 5 29 5208 6 8d-0 2 , 
+(  . 3328274163 3 18927d-02 , 
+  (  - . 2 83 79404 64 10067 3d-0 3 , 
+  ( - . 7 6 15 127 119 2 9 09 3 4d-0 3 , 
+  (-. 23 9007037 2 062 00 2d-0 3 , 
+(  . 423 674 652927213 3 d-04 , 
+  (-. 937 299784 02 5 66 60d-0 3 , 
+(-. 6275701960924115d-03, 
+  (  . 184903 88 6388677 ld-03 , 
+(  . 3795955396627962 d-03 , 
+(  . 1370678212350365d-03 , 
+ (-. 2256873 55 157875 ld-04 , 


.  4537367583220980d-ll) , 
.  18973 33806124 500d-12 )  , 
.3371533582241440d-ll) , 
.  3 102 3 3 252 8113 870d- 10 )  , 
.2412 173 021438 3 0Qd-09 )  , 
.  3 114846206597 Q60d-09)  , 
.1779089985273  050d-07) / 


. 10 0 179 6 15 9 8 715 lOd- 16) , 
.  6331537179006460d-16) , 
.  32100 68 39563 3 560d-15) , 
.  63 1628 2 8 66 674 OlOd- 14 )  , 
.  2654113441221400d-13 )  , 
.  6757485415387840d-12)  , 
.2157 677521550440d-13 )  , 
. 39 299 27 3 63 489 88 Od- 13  )  , 
.  252138C004952490d-ll)  , 
.  9523018697274730d-ll)  , 
. 737096 12 65482 4 OOd- 10) , 
.  25037893 08434 250d-09) / 


.  3756978 60 60 1592 Od- 17)  , 
.  18 5 2 643 77 50 668 9 Od-l 5 )  , 
.  933 1802  48739 1570d.-16)  , 
. 4 468 2 0079 3 2 64 19 Od- 15) , 
.  8205527 812277 57 Od-14 ) , 
.  3393 508378723 270d- 13 )  , 
.  118 877 117 2 507 OOOd- 14 )  , 
.  28585622523 90 12 Od- 13 )  , 
.  103 53 588842 88 OSOd- 12 )  , 
.  3 0584 2 9 6779 264 lOd- 11)  , 
.  11013 067008 14700d- 10)  , 
.  8 64 984890787 19 60d- 10)  / 


. 144 4 64 103 5 6 52 92 Od- 15)  , 
. 6647 2 3 073 0883 52 Od-13) , 
.  115 567 IS 027 14 180d- 13 )  , 
.  4346463579853 540d-13 )  , 
. 12 4 43 4 2 5722 14 57 Od- 11) , 
. 4472204080723  670d- 11)  , 
,  1904932 268 853 110d-ll)  , 
.  3 40255509 23 0573 Od-12)  , 
.  3068433 18 2874 190d- 10)  , 
.  2 415609997 2 46610d-09)  , 
.  78538127087  82 050d-09 )  , 
.  1790520486266400d-07)  / 


data  COS120/ (-.5000000 0000 00 OOOd  00,  .8660254037844386d  00)/ 
data  CO830  /(  .3660254037844386d  00,  . 5000000000000000d  00)/ 
data  third/  .3333333333333333d  00/ 
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CHIVER 

This  program  calculates  Che  partial  wave  series  solution  for  the 
backscattaring  form  function  of  a  tungsten  carbide  sphere.  The  program  was 
adapted  from  one  written  by  R.  C.  Chiver  at  al.  and  referenced  in  the  paper 
given  at  the  beginning  of  the  program.  The  program  results  were  checked 
against  those  of  G.  C.  Gaunaurd  et  al.  in  Ref.  11  of  Chapter  3  in  the  ka 
region  from  1  £  ka  £  70.  The  form  function  of  other  spheres  can  be 
calculated  by  changing  the  material  parameters. 
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CCCCCCCCCCCCCCC 
C  C 
C  CHIVES  C 
C  C 
CCCCCCCCCCCCCCC 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c  c 

C  THIS  PROGRAM  CALCULATES  FORM  FUNCTIONS  USING  PROGRAMS  REFERENCED  IN  C 
C  R.C.  CHI VERS  AND  L.W.  ANSON  "CALCULATIONS  OF  THE  BACHS CATTERING  C 
C  AND  RADIATION  FORCE  FUNCTIONS  OF  SPHERICAL  TARGETS  FOR  USE  IN  ULTRA-  C 
C  SONIC  BEAM  ASSESSMENT. "ULTRASONICS , JAN.  1982 ,PP. 25-34. THE  PROGRAM  C 
C  OUTPUTS  THE  KA  OF  THE  SPHERE  AND  THE  ABSOLUTE  VALUE  OF  THE  FORM  C 
C  FUNCTION  AT  THAT  KA.I  FOUND  THAT  THE  PROGRAM  REQUIRES  HIGH  PRECISION  C 
C  ARITHMETIC  OR  IT  WILL  HAVE  NUMERICAL  DIFFICULTIES. THE  PROGRAM  ALSO  C 
C  BECOMES  EXPENSIVE- I?  TOU  DESIRE  TO  GO  TO  VERY  LARGE  HA’S  (>30) .  C 
C  C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIMENSION  PB(300) ,PCB<300) ,PS8(300) ,KAM(2000) ,FABS(2000) 

REAL  HAM 
COMPLEX  CFAC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  SET  MATERIAL  PARAMETERS  FOR  THE  SPHERE.  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DATA  RHO , VC , VS  /7. 7, 3960. 0,3240.0/ 

RH-1.0/(RHO*2.0) 

RAT1-U83 . 0/VC 
RAT2-1483.0/VS 
WRITE (6,5)  RHO, VC, VS 

5  FORMAT ( 3X , 5H  HHO-,FS.2,4H  C1-,F8.1,4H  C2«,F8.1) 

WnT'TT  t  k  A.  \ 

4  FORMAT ( 11X , 3H  RA,9X,3H  /ff/) 

XMAX-50 

X- 1.0 
1C-0 

10  x-x+o.os 

K-K+l 

T-1.25*X 

XI- RAT1*X 
X2-RAT2*X 
X2S-X2**2 
XSN-SIN(X) 

XCS-COS (X) 

CALL  BESS (X, ?S,T) 

TERM" ( XSN/X) /PB ( 1 ) 

BN1-PB(2) "TERM 
CALL  BESS (XI, PCS, T) 

CBN1-PCB ( 2 ) /PCS ( 1 ) 

CALL  BESS (X2 , PS3, T) 

SNNO— XCS/X 
XSIGN-1.0 
3 BET-0 . 0 
SALP-0.0 
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N-0 

SUM-0 . 0 

CALL  STVAL(X, XI , X2 , ALPN, BETH, CSN1, BH1 , SNN1 ,RH, XSN , X2S ,  XCS ) 

20  C0NTINU2 
CBN1-1.0 

CSN2-PCB (N+3 ) /PCS (N+2 ) 

Q- FLOAT (N+l) 

AHUM- (Q*C3N1 ) - (X1*CBN2 ) 

ADEN* ( (Q-1.0) *C3N1) -(X1*CBN2) 

DNUM-  ( ( (X2S/2 . 0 )  -0*  (Q-l .  0 ) )  *CBN1 )  -  ( 2 . 0*X1*CBN2 ) 

SBN1-1.0 

S8N2-P3B (N+3 ) /PSB (N+2 ) 

BNUM-(2.0*Q*(Q+1) )  *S8N1 

EDEN-( ((2.0*(Q**2) )-(X2S+2.0)  )  *S8N1)+(2.0*SBN2*X2) 
ENUM-2.0*O*(Q+l-0)  *(  ( (1.0-<J)  *SBN1)  +  (X2*SSN2  )  ) 

FN1-X2S*RH*( (ANUM/ADEN) -(BNUM/EDEN) ) / ( (DNUM/ ADEN) 

6- (ENUM/EDEN) ) 

30  CONTINUE 

BN2-PB (N+3 ) *TERM 

<ari-(  (FNl-C)  *BN1)+(X*BN2) 

BN1-BN2 

SNN2- ( ( (2 . 0*0) +1. 0) *SNN1/X) -3NNO 
HN1- ( ( FNl-O) *SNN1 ) + (X*SNN2 ) 

SNNO-SNN1 

SNN1-SNN3 

DEN- (GN1**2 ) + (HN1**2 ) 

ALPN1— (GN1**2)/DEN 
8ETN1—  (<SN1*HN1)/0EN 

ADD-Q* (ALPN+ALPN1+ (2 . 0*ALPN*AL?N1)  +  ( 2 . 0  * BETN*  BETN1 ) ) 
FAC-XSIGN*(2.0*Q-1.0) 

SALP-SALS»+FAC*ALPN 
SBET-SBET+FAC* BETN 
ALPN-ALPN1 
BETN-BETN1 
SUM— SUM+ADO 

IF  ( (AB3( ADD/SUM) ) .LT.1.0E-8)  00  TO  90 

J-0 

N-N+l 

XSIGN— XSIGN 
00  TO  20 
90  J-J+l 

IF  (J.EQ.5)  00  TO  100 
N-N+l 

XSIGN— XSIGN 
GO  TO  20 

100  TR-(-4.0/(X**2) ) *SUM 
FAC— XSIGN*  (2. 0*0+ 1.0) 

SALP-SALP+FAC+ALPN1 
SBET-S8ET+FAC+BETN1 
CFAC-CMPLX ( SBET , -SALP) * ( 2 . 0/X) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 

C  THE  VARIABLE  FINF  IS  THE  ABSOLTE  VALUE  OF  THE  BACKS CATTERING 
C  FORM  FUNCTION. THE  VALUE  OF  (KA)  AND  (FINF)  ARE  OUTPUT  BELOW. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


noon 


non  non 
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riKT-CABS (CPAC) 

KAM(X) -X 

pabs(X)-fin? 

WRITER, 40)  KAM(X) , PASS (X) 
40  FORMAT ( 5X , F10 . 2 , 2F14 . 4 ) 
IF(X.GS.XMAX)  SO  TO  SO 
GO  TO  10 
50  CONTINUE 
STOP 
END 


SUBROUTINE  STVAL ( X , XI , X2  ,  AI*PN ,  SETN ,  CSN1, BN1, SNN1,RH, XSN, X2S 

■  #  aCj) 

TERM»CBN1*X1 

PN— X2S*RH*TERM/  ( (X2S/2 . 0)  -2 .  0*TERM) 

SNN1—  (  (XCS/  (X**2 )  )  +  ( XSN/X)  ) 

GN- ( PN* ( XSN/X) ) + ( SN1 * X) 

HN- ( - ( PN*XCS )  /X)  +  ( X*SNN1 ) 

DEN- (GN**2 ) + (HN**2 ) 

A1PN—  (GN*«*2 )  /DEN 
BETN—  (GN*HN)  /DEN 
RETURN 
END 


SUBROUTINE  BESS  (X,  PA,  T) 

DIMENSION  PA(300) 

L-IPIX(T)+13 

PA(L+2)«0.0 

PA(L+l)-l.0E-3Q 

DO  10  I-1,L 

H-(L+1)-I 

R-PLOAT(M) 

PA(M)-( ( (2  .  Q*R)  + 1. 0) *PA(M+1)/X) -PA(M+2) 

CONTINUE 

RETURN 


